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Abstract 


This report describes the high speed Dynamic Data Acquisition System 
(DDAS) which provides the capability for the simultaneous measurement of 
velocity, density, and total temperature fluctuations. The system of 
hardware and software is described in context of the wind tunnel 
environment . 

The DDAS replaces both a recording mechanism and a separate data 
processing system. The data acquisition and data reduction process has 
been combined within DDAS. DDAS receives input from hot wires and 
anemometers, amplifies and filters the signals with computer controlled 
modules, and converts the analog signals to digital with real-time 
simultaneous digitization followed by digital recording on disk or tape. 
Automatic acquisition (either from a computer link to an existing wind 
tunnel acquisition system, or from data acquisition facilities within 
DDAS) collects necessary calibration and environment data. The generation 
of hot wire sensitivities is done in DDAS, as is the application of 
sensitivities to the hot wire data to generate turbulence quantities. The 
presentation of the raw and processed data, in terms of root mean square 
values of velocity, density and temperature, and the processing of the 
spectral data is accomplished on demand in near -real -time with DDAS. 

This paper describes the interface to DDAS and the internal mechanisms of 
DDAS. A summary of operations relevant to the use of the DDAS is also 
provided. 
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Symbols 


Ai - A g 

E 

e' 


u 



Constants in equation (1) 
mean voltage across wire 

instantaneous voltage across wire (less the mean) 
instrumentation amplifier scalar 


velocity sensitivity 


3 log e 

3 log u p,T 0 ,T w 


density sensitivity 


3 log e 

3 log p u, T 0 ,T w 


temperature sensitivity 


3 log e 

3 log T 0 u,p,T w 


mean total temperature 


mean temperature of heated wire 


u mean velocity 

p mean density 


1. INTRODUCTION 

Recent advancements have been made in hot wire anemometry techniques which 
allow a three wire probe to separate three components of the perturbations 
in the flow field. Velocity, density and total temperature fluctuations 
can be determined as a function of three parallel hot wires, since at 
subsonic and transonic speeds it is generally conceded that the voltage 
measured across a heated wire mounted normal to the flow and operated with 
a constant temperature anemometer is a function of velocity, density and 

l 

total temperature. Under these conditions, a single equation is obtained 
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for the fluctuating voltage across a single wire which is a function of 
the three variables - velocity, density, and total temperature. 
Quantitative measurements for the three fluctuations in the flow variables 
have used probes with three wires mounted normal to the flow and operated 
at three different "overheats". 

The development of a dedicated hardware and software system to support hot 
wire anemometry at NASA Langley Research Center in the Fluid Dynamics 
Branch of the Transonic Aeronautics Division was precipitated by the 
necessity to process simultaneous hot wire data from three wire probes 
more rapidly than previously possible. 

Prior to the development of the DDAS, all data was acquired on FM tape, 
and' all processing was done in an off-line batch mode. This method 
delayed recognition of faulty or incomplete data, and test results were 
often delayed several months. 

During a flow diagnostics test in the 8 Foot Transonic Pressure Tunnel 

(8'TPT) at NASA Langley Research Center in. January of 1988 the DDAS was 
connected in parallel to the existing test instrumentation systems to 
provide an initial test bed for the new system. See Figure 1. The DDAS 
was not designated as a primary data acquisition or reduction system, but 
it soon became apparent that the data logging capabilities would be 
especially helpful in collecting the hot wire calibration data in an 
easily manageable format. The hot wire calibration data and the 
generation of hot wire sensitivities were processed only by the DDAS and 
the calibration data and sensitivities were used both by DDAS and by other 
data processing facilities. As a test of the digitization and recording 
capability, dynamic data was routinely digitized in parallel with the FM 
tape recordings. As soon as adequate calibration data was collected, the 
DDAS processed some of the data, and provided velocity, density, and total 
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temperature turbulence measurements. These results compared favorably 
with subsequent off-line batch data processing. 

A second test was supported to compare hot wire techniques and laser 
velocimetry techniques in the Basic Aerodynamic Research Facility 

(B.A.R.F.) , 3 ’ 5 


The DDAS provides the processes necessary to: 

1) acquire the hot wire calibration data 

2) acquire the dynamic hot wire data 

3) generate hot wire coefficients and sensitivities 

4) compute velocity, density, and total temperature fluctuations 

5) compute other statistical relationships 

6) provide spectral analysis 

7) manage data 

8) produce reports and plots 

2. SYSTEM DESCRIPTION 

DDAS is a system of hardware and software based on systems purchased from 
Data Laboratories, Ltd., Precision Filters, Inc., and Hewlett Packard 
Corporation. Modifications and enhancements to the software and hardware 
have converted a waveform recorder into a hot wire anemometry acquisition 
and processing system specifically tailored to the three wire technique 
that yields separate velocity, density and total temperature components of 
turbulence . 


2.1. HARDWARE 
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The system (fig. 1) is divided into an analog front end, and a computer- 
based processing and display section. The analog front end consists of a 
filter/amplifier subsystem, a high speed digitizer, and a low speed 
digitizer (or a data link to another acquisition computer) . All are fully 
computer controlled. The processing and display subsystem controls the 
analog subsystems, and receives the digitized data, processes the data, 
displays the data, and stores the data in a permanent file. 

2.1.1. Amplifier and Filter Subsystem 

The analog signals from the hot wire anemometers are first routed to the 
Precision Filters, Inc. precision amplifier and filter subsystem. This 
subsystem is currently configured for four channels, providing support 
for only one three wire probe. Each channel successively passes the 
anemometer signal through a pre amplifier, high pass filter, low pass 
filter, and a post amplifier. The full bandwidth capability of each 
channel is .1Hz to 200KHz, but the high pass and low pass filters usually 
provide a narrower bandwidth (1Hz to 5KHz). The high pass filter acts as 
the anti-aliasing filter for the high speed digitizer. 

2.1.2. High Speed Digitzer Subsystem 

A high speed digitizer, called a Multitrap modular waveform recorder by 
Data Laboratories, Ltd., is configured to digitize up to 14 channels of 
fluctuating hot wire data (three channels are required for each 3 -wire 
probe) at rates of up to one million (1M) samples per second (6 channels 
at up to 1M samples/sec, and 8 channels at up up to 256 thousand (256K) 
samples/second. This data is stored temporarily in the Multitrap buffer 
memories (up to 256,000 samples per channel), and then transferred to the 
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HP 9000/330 computer at about 100,000 samples per second - one channel 
memory at a time - via a dedicated 16 bit parallel bus (GPIO) . 

2.1.3. Low speed digitizer 

This subsystem is not currently implemented, and an existing wind tunnel 
data acquisition system provided DDAS the functions of a low speed 
digitizer subsystem. However, the optional low speed digitizer subsystem 
would consist of a multiplexer and digitizer selected for collection of 
mean values, not fluctuating values. This subsystem would be used to 
collect calibration data and tunnel parameter data, which would be logged 
for further processing of the calibration and fluctuating data. 

2.1.4. Computer link to tunnel computer 

The General Purpose Interface Bus (GPIB) , an IEEE488 standard bus, is used 
to receive static data from the existing tunnel data acquisition system 
computer. The tunnel computer transmits a packet of data relevant to the 
tunnel conditions and hot wire calibration data. This link was selected 
because of its availability in both the tunnel computer and in the DDAS 
computer. It provides an 82,000 byte per second transfer rate, which is 
more than adequate to receive as many as four complete ASCII data packets 
per second. The tunnel computer actually sent only one packet per second. 

Use of the existing tunnel data acquisition system to collect the tunnel 
conditions and the additional mean values related to the hot wire 
calibration data eliminates the need for a parallel hardware system (the 
low speed digitizer) , and the need to develop instrument calibration 
software and hardware. It does, however, provide additional work for the 
tunnel computer personnel to configure their acquisition setup to handle 
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the additional channels, and to provide the GPIB software to generate the 
data packets for the DDAS . 

The link is configured with the DDAS end as not system controller, and as 
device 01. This was accomplished by setting switches on the HP 98624A HP- 
IB Interface Card inserted into the computer specifically for the link. 

The select code was set to 8; interrupts are not relevant, since they are 
not used. The wind tunnel computer providing the data packet is 
configured as system controller, and outputs ASCII data packets at a rate 
set by the wind tunnel computer. 

The packet is read into a DDAS packet buffer with one program statement in 
the subroutine Get_packet in module MODUS R2 : ENTER Pkt_sc ; Pkt$ (*) , where 
Pkt_sc is equal to 8, and the Pkt$ array was sized for 47 each 80 
character strings. 

The packet format is shown in Table 1. 

2.1.5. Computer Peripherals 

2. 1.5.1. Disk storage 

75Mb of non removable disk storage is available for programs and data. In 
addition, a 1.2Mb removable disk drive is available for program 
development and hot wire calibration data. 

2. 1.5. 1.1. 55Mb hard disk 

The computer then transfers the data to a 55Mb hard disc at about the same 
100,000 samples per second - one channel at a time. 
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2. 1.5. 1.2. 20Mb hard disk 

Programs and support software is stored on the 20Mb hard disk. 

2. 1.5. 2. Tape storage 

Once the data disc is full, the data is copied to a 67Mb tape cartridge 
for permanent storage. 

2. 1.5.3. Display 

A color CRT is the system console and data display. 

2. 1.5.4. Plotter 

An 8 pen autoload flatbed plotter is available for plot generation, and is 
used to display dynamic data and hot wire calibration data. 

2. 1.5. 5. Printer 

A dot matrix printer is available for data display. It can produce screen 
dumps , but is used primarily to generate a record of the hot wire 
calibration data, and, as data processing is accomplished, the results of 
the processing are printed. 


2.2. SOFTWARE 


2.2.1. Software environment 



Page 9 


All DDAS programs operate under a BASIC operating system, in an 
interpretive BASIC language. Several compiled subroutines are a part of 
the ACQUIRE software system to enhance computational speed in some parts 
of the software. 

2.2.2. Baseline software 

The ACQUIRE^ software system, provided by Data Laboratories, Ltd. is the 
basis for the Dynamic Data Acquisition System (DDAS) used for the 
acquisition of hot wire anemometry data. 

Since the ACQUIRE package is an off-the-shelf product, no attempt to 
describe its full capability will be made. 

ACQUIRE has been slightly modified in only one area - the addition of a 
sequence number to a dynamic data disk file was inhibited if it was not 
necessary to discriminate between two files with the same name. See 
Section 2 . 2 . 6 . 2 . 1 . 2 . f or further discussion on the requirement for the 
modification. 

2.2.3. Hot wire application software 

Major additions were made to ACQUIRE in the form of two sub-programs: 
M0DUSR1 and M0DUSR2. These user written modules are configured according 
to guidelines provided by ACQUIRE, so that they will be automatically 
included in ACQUIRE. Appendix A contains the full program listings of 
these two modules . 

The functions implemented in both M0DUSR1 and M0DUSR2 are listed in 
Appendix B. 
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2.2.4. Configuration files 

Configuration files used by ACQUIRE for the 8'TPT test include acquisition 
setup parameters, hardware configuration parameters, and default display 
parameters. They are set, saved and stored by a variety of ACQUIRE 
functions . 

2.2.5. Sequence program files 

The sequence program functions of ACQUIRE provide a mechanism for 
specifying a series of functions to be accomplished. Both the 
initialization and acquisition sequence program files used to tailor the 
DDAS for support of the 8'TPT test are listed in Appendix C. 

Although the configuration files and sequence files are not the easiest to 
configure, the result is a system that is literally a "turnkey” system. 
Turn on the hardware, allow the hardware and software to be configured, 
and press a button to simultaneously acquire calibration and dynamic data. 

2.2.6. Hot wire data acquisition 

The ACQUIRE software system was augmented to support the specific 
requirements of hot wire anemometry systems currently in use at NASA 
Langley Research Center in the Fluid Dynamics Branch of the Transonic 
Aeronautics Division. Software design and implementation followed both 
form and style of the supplied ACQUIRE software, maintaining the 
appearance of a seamless environment within ACQUIRE. This feature 
resulted in a software system for an instrument that has evolved from a 
waveform recorder to a hot wire anemometry system designed to acquire and 
process both mean hot wire calibration data and fluctuating hot wire data. 
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The end result has been an easily used flow diagnostics instrument that 
minimizes the researchers workload. 

With hot wires, there are most often two concurrent tasks: calibration of 
the wires, and acquisition of the dynamic, or fluctuating data. This 
system is designed to calibrate and process three wire probe data. The 
current implementation supports two three wire probes, and acquires, 
processes, and stores them separately. 

Once the instrument is configured, each data point is acquired with the 
push of a single button -- one button recording. Data processing does 
require a few more button sequences, but only because the researcher 
provides more direction in the data reduction process. 

For the 8'TPT flow diagnostics test in January of 1988, the hot wires had 
not been previously calibrated, so concurrent calibration data and dynamic 
data acquisition was necessary. 

2. 2. 6.1. Calibration 

Calibration of three wire probes require a complete data system to acquire 
mean (static) conditions of both the operational environment and of the 
mean hot wire values. 

It is assumed that the hot wires are sensitive to velocity, density, and 
i 

total temperature. To determine what the sensitivity is to each variable, 
the mean voltage output of each hot wire must be measured at each 
combination of velocity, density, and temperature. The tunnel run 
schedule was configured to assure that adequate data jpoints are taken to 
provide a realistic profile of sensitivities. 
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The run schedule was also selected to expose the hot wire probe to the 
highest dynamic pressures first, so that if a wire is going to break, then 
the least amount of tunnel time will be lost. 

Each of the three parallel hot wires on a probe are operated at different 

overheats, to encourage a wide separation of sensitivity between each hot 
1 

wire . 

The calibration data consists of mean values, which do not require the 
high sampling rates normally invoked to digitize the dynamic data, so the 
calibration data is not acquired through the high speed digitizers , An 
existing wind tunnel data acquisition system collected the data, and then 
transferred it through a GPIB link to the DDAS computer. Several data 
packets of data are averaged and then stored in a formatted record on 
disk. 

2. 2. 6. 1.1. Observation files 

This calibration data is stored in a file called an observation file. 

Other related mean data is also stored in the observation file: 
tunnel conditions, 
test identification parameters, 

auxiliary data (such as RMS microphone readings, and 

amplifier gain settings), and 
simple calculated data (such as density, velocity, 

static temperature, the logs 
of a variety of data, and the 
products of the logs of a 
variety of data. 


! 
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The observation "file' 1 is really several files. The first is a file 
containing the data received from a packet sent from the MODCOMP via the 
HP-IB link. (Actually, several packets are averaged, and the average is 
stores as an observation in the first observation file.) The second and 
third files each contain data related only to a specific three wire probe. 
These files also contain tunnel condition and test identification data as 
well, but only the tunnel computer data and simple, computed values 
related to a specific probe is contained in these files. 

The internal format of the observation files was carefully selected to 
conform to the format specifications of an existing statistics package. 

The Basic Data Statistics package from HP has historically been the 
statistics package used to reduce the hot wire calibration data, so the 
file format was made compatible with that package. 


2. 2. 6. 2. Dynamic data 

The acquisition of the dynamic hot wire data is entirely accomplished by 
the off-the-shelf ACQUIRE software. ACQUIRE has all the mechanisms 
necessary to configure the actual data acquisition hardware and the 
capability to manage the data once it has been digitized and buffered by 
the high speed digitizer hardware, which includes the transfer from the 
buffer to computer memory, and the transfer of the data to disk. These 
functional modules are "strung together" in a sequence program mechanism, 
which is also an inherent part of ACQUIRE. 


2. 2. 6. 2.1. Fluctuating data 
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Fluctuating data is either digitized data from the hot wire anemometers, 
or it is calculated data from the process of computing velocity, density, 
and temperature fluctuations, which is discussed later. In either case, 
the result is a time varying array of samples of data. 


2. 2. 6. 2. 1.1. Fluctuating data volume 

The sheer volume of fluctuating data is worth note: since each channel can 
handle 256K samples (512K bytes) at a time, and there are 3 channels per 
probe, and the same amount of results exist, 256K X 2 X 3 X 2 = 31457K 
bytes per observation. For 117 observations (8' Transonic Pressure Tunnel 
Test 934), 368M bytes of data storage becomes necessary. 

2. 2. 6. 2. 1.2. Fluctuating data disk file naming convention 

DDAS collects and generates multiple channels of fluctuating data files 
for each test condition or observation. These files are related to a 
specific record in an observation file, which contains non- fluctuating 
scalar data related to the observation. The relationship of the 
fluctuating data file names to the test condition and to the observation 
file and record within the observation file is specifically defined by 
convention. The use of a naming convention allows data reduction programs 
to associate all data files necessary for data reduction and for naming 
resultant fluctuating data files. Table 2 details the naming format. 

The dynamic data samples digitized by the MULTITRAP digitizer hardware are 
stored by the waveform recorder function of ACQUIRE on a channel -per- file 
basis. Fluctuating data names are generated whenever a hot wire 
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calibration observation is logged, so that a naming convention is followed 
- should a request to digitize hot wire fluctuating data be processed. 

2.2.7. Coefficient calculations 

The process of providing the sensitivities necessary to convert three hot 
wire data arrays into velocity, density and temperature fluctuation arrays 
first requires that the calibration data be processed by multiple linear 
regression techniques to produce a set of coefficients for each of the 
three hot wires. Since the relationship of the performance of the hot 
wire is highly nonlinear in relationship to the velocity, density, and 
temperature, up to 10 coefficients are required (Eight are in use, as 

shown ^ : 

log E - Aj+ A 2 log u + A 3 log p + A 4 log T 0 

o 

+A S log u log p + A 6 log u log T 0 
+ A 7 log p log T 0 

+ A 8 log u log p log T 0 (1) 

Since velocity, density, and temperature are all known for each 
observation (as collected by the DDAS form the tunnel data acquisition 
system - in the form of P^, P g , and T T ) , the most direct solution is 

through multiple linear regression. 

Whenever requested by the operator, a Multiple Linear Regression routine 
(which is a specifically modified version of the Hewlett Packard routine 
MLR which was purchased as part of a statistics package) is invoked, which 
calculates coefficients for each hot wire on each probe. These 
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coefficients can then be stored in a coefficient disk file related to each 
probe. (This internal MLR routine is not currently implemented.) 

Alternatively, coefficients calculated in a separate multiple linear 
regression package may be read from a disk file generated by that package, 
or the coefficients may be manually entered through the keyboard. 

2. 2. 7.1. Sensitivity calculations 

The coefficients, which represent the hot wire relationship to velocity, 
density, and temperature, are combined with specific test conditions, 
which have been stored in an observation record of the observation file. 

s u - A 2 + A s log p + A 6 log T 0 
+ A g log p log T 0 
S p - A 3 + A 6 log u + A 7 log T 0 
+ A 8 log u log T 0 
S T 0 “ + A 6 log u + A 7 log p 

+ Ag log u log p 

New ACQUIRE functions created in M0DUSR2 allow the appropriate calibration 
file to be specified, and the beginning and ending observations and 
beginning and ending probes to be selected for the computations. For each 
observation and each probe, the log values of velocity, density and 
temperature are retrieved from the appropriate record in the observation 
file. Once the computation is completed for each probe, the resultant 
senstivitys are inserted into existing, but as yet unused variables in the 
previously recorded observation. 


(2a) 

(2b) 

(2c) 
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2. 2. 7. 2. Calculating velocity, density, and temperature fluctuations 

Once the sensitivities have been calculated, the operator may request that 
the dynamic data for a given set of observations and probes be processed 
in a way that yields dynamic waveforms representing fluctuating velocity, 
density and temperature (instead of 3 fluctuating voltages) and with 
turbulence figures and other statistical performance characteristics. 

The equation that defines the relationship of voltages to turbulence 
parameters is\ 
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By rearranging again, which involves inverting the sensitivity matrix, 
solve for the three unknowns: 



The computation of the instantaneous velocity, density, and temperature is 
accomplished as shown: 


for each probe and each observation to be processed, 
for each of the three hot wires 

the mean hot wire voltage (E) is retrieved 

the gain (G) for the fluctuating hot wire voltage is retrieved 
the three sensitivities (u,p,T 0 ) are retrieved and placed in the 

sensitivity matrix 

the dynamic data file is retrieved from disk, and placed in memory 
the sensitivity matrix is inverted 
for each of the instantaneous samples 
for each of the three hot wires 
compute : 


e ' 1 
E G, 


w 


and store in the independent variable matrix 
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matrix multiply the inverted sensitivity array by the 
independent variable array, and place the instantaneous 


T t 

o 

— \ 


turbulence ratios (tt, , and in memory 

u p i 0 


u # o' io 

compute the RMS values of , and re- 


store the RMS values of velocity, density, and temperature 

for each observation into existing, but as yet 


u' 

U 


2 - !•■ 


P ’ T 0 

unused variables in the previously recorded observation 


store the fluctuating velocity, density, and temperature waveforms in 
disk files for later spectral investigations, utilizing existing 
functions of ACQUIRE 


2.2.8. Precision filter system control 

The computer control of the filters and amplifiers allows adaptive 
processing of various hot wire signals, which are dependent upon a variety 
of operational parameters. The software interface allows full control of 
each functional module within the Precision Filter system - including the 
calibration module, and also allows the interrogation of all status and 
condition data - including the calibration module. The modules are 
connected to provide a full calibration sequence, and to allow full 
operator control, semiautomatic or automatic operation. (This feature not 
implemented as of April 88). 
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2.2.9. Utility Functions 

Other utility functions were implemented to enhance operational 
characteristics of the system. The ability to eject plots, and the ability 
to plot "special” hot wire calibration data, are "plot utilitys" . The 
ability to list categories of files on the disk, to purge extraneous files 
- or groups of files, the ability to copy or move files - or groups of 
files - to another disk (or tape) are "file utilitys". 

Function LOGFILLE TO PC was used to transfer observation files containing 
logged and computed data to another system. A GPIB bus connects DDAS to 
the PC, where a GPIB card (National Instruments or HP) is installed. 
Appendix D contains the PC BASIC program used with the HP GPIB card to 
receive and store the data on the PC disk. 

3. OPERATION 

Operation of DDAS begins prior to the tunnel operation. Configuration of 
both the ACQUIRE software and the high speed digitizer is accomplished 
from within the software. Configuration files of various types are 
generated by the software and are recallable either at power on time or 
from within a sequence program. 

3.1. System setup - hardware 

The initial configuration of hardware is accomplished only once. The 
assignment of device addresses is as follows: 

GPIB devices: 

printer 701 

20Mb program disk 703,0 
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67Mb data tape 

703,1 

plotter 

705 

high speed digitizer 

708 

(control link) 


computer link 

801 

55Mb data disk 

1400 


GflO (16 bit parallel) devices: 

high speed digitizer 12 
(data link) 

3.2. System setup - software 

The ACQUIRE operating system was configured to operate within the memory 
constraints of 4Mbytes, and to be configured for 14 channels of 
digitizers, and 14 channels of data memory. Refer to the ACQUIRE 
operations manual for further details. 

3.3. ACQUIRE installation 

Upon receipt of the software, the installation procedure defined by the 
manufacturer allows the software to be configured for existing hardware, 
including amount of computer memory, number of digitizers in the digitizer 
chassis, and the maximum number of channels in the computer memory at any 
one time. 


3.4. System variables 
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A file structure is maintained in the ACQUIRE software for containing a 
wide variety of currently selected operating parameters. This mechanism 
allows the operator to interactively select preferred operational 
conditions, and then store the "sysvars" on disk for later retrieval. 

These parameters include, but are not limited to, display format, memory 
length for each channel, waveform file names, waveform channel selection, 
system variables file name, binary switch name, plot file name, and 
sequence file name. To recall a specific set of system variables 
automatically at power on time, the system variables are stored in a file 
called "AUTO VARS" . 

3.5. Binary switches - digitizer configuration - MULTITRAP 

The high speed digitizer is configured utilizing an interactive session to 
select sampling rates, gains, trigger modes, data block size, etc., and 
then the configuration of the binary switches within the digitizer are 
saved in a binary switch configuration file. To recall a specific 
configuration for the digitizer automatically at power on time, the binary 
switches are stored in a file called "AUTOSW" . 

3.6. Plot setup 

The format of a plot is defined interactively and may then be saved on a 
plot file. The actual data is not saved in the file. Once the waveform 
channel in memory is selected, the position, scaling, and labeling of the 
axis is defined, and waveform labeling is determined. Once all channels 
are positioned and defined, the plot title is defined, and the plot 
configuration is saved to disk. 


3.7. Sequence program - initialization 
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An initialization sequence program is interactively generated, which will 
determine a sequence of functions to be performed to set up DDAS to a 
configuration relevant to a specific wind tunnel test. See Appendix C for 
a listing of the initialization sequence used for the 8'TPT test. To 
recall the initialization sequence program at power on time, the sequence 
is stored in a file called " AUTOS EQ " . 

3.8. Sequence program - acquisition 

A run sequence program is interactively generated, which will determine a 
sequence of functions to be performed to: 
log calibration data 
digitize and store fluctuating data 

plot hot wire mean voltages vs. mass flow, pu, (see Fig. 2) 
print a report displaying many parameters of the current observation 
whenever the operator presses a single button. See Appendix C for a 
listing of the run sequence used for the 8'TPT test. The run sequence 
program is automatically loaded by the initialization sequence program, so 
that once all configurations are defined, powering on the system, and 
pressing a button is all that is necessary to simultaneously acquire both 
calibration and fluctuating hot wire data. 

3.9. File transfer to PC 

The probe log data files - one or both - can be transferred to a PC via a 
dedicated GPIB cabled between DDAS and a PC. The PC BASIC program 
,, XFR.HP ,f (see Appendix D) should be started first, and then, before 
providing the requested file name, invoke the DDAS function LOGFILE TO PC. 
Refer to* the relevant function sheet in Appendix B for details on proper 
configuration prior to starting the transfer. When the file name is then 
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entered into the PC, which defines where the data is to be stored, the 
transfer will begin. 

4. SYSTEM STRENGTHS 

4.1. ACQUIRE 

ACQUIRE, in combination with the hardware is a very versatile waveform 
recorder: 

It controls all the hardware associated with the system. 

It manages hardware and software configuration - via files. 

It manages process, or "sequence" files. 

It manages dynamic data files and internal arrays of data. 

It provides a choice of operator dialogue techniques, including: 
cursor, menu, and command line entry. 

It provides data display management. 

It provides the waveform plotting capabilities. 

A Digital Signal Processing package is included which provides: 

Fast Fourier Transforms 
filters 

power spectrum 
transfer functions 

The acquisition of the dynamic data, and the storing of the dynamic data 
is a very significant strength of ACQUIRE. But most importantly, the 
internal design allowed application routines to be written into ACQUIRE, 
which produces a set of software that appears to the user to be a single 
entity, without seams, and fully integrated. 
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4.2. Data logging 

The internal log file format allows direct access by a commercially 
available statistics package, which includes a multiple linear regression 
analysis capability necessary to generate coefficients used in creating 
hot wire sensitivities. 

4.3. Computer link 

4.3.1. PDAS to tunnel computer 

A software/hardware link is currently used with the MODCOMP data 
acquisition computer to receive mean data values, but a self-contained, 
accurate and reliable static data acquisition subsystem could be 
integrated, making the DDAS self-contained. The use of an existing data 
acquisition system for the collection of mean values transfers the 
instrument calibration requirements for those values to another system. 


4.3.2. DDAS to PC 

The LOGFILE TO PC function to transmit the logged and computed parameters 
to another system, where the data is reformatted and imported to a 
spreadsheet program (Lotus Symphony) for further analysis and data 
presentation. 


5. LIMITATIONS 


5.1. Uncalibrated wires 
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Hot wire calibration currently consumes the major portion of the tunnel 
operation time. Although not a limitation of the DDAS , the process of 
calibrating hot wire probes relative to temperature, density and pressure 
is currently the most expensive part of the three wire technique. 

Pre-calibrated wires would allow real time processing of the voltages from 
the three wires into the velocity, density, and temperature components of 
turbulence. The facility would be much less expensive to construct and 
operate than the wind tunnel to be supported, since the size could be much 
smaller, and the tolerable turbulence levels could be higher, since only 
the mean values of velocity, density, and temperature are used in 
determining hot wire sensitivities. 

For wind tunnels not capable of independently controlling velocity, 
density (or total pressure) and temperature, the three wire technique 
requires that the wires be pre-calibrated, since the sensitivities could 
not be properly determined in such a wind tunnel. 

Although a hot wire calibration tunnel has been partially constructed, it 
is not yet operational due to manpower and funding constraints. A data 
logging program module developed for DDAS is available as a module for 
eventual integration with an instrumentation system expressly for the hot 
wire calibration facility. 

5.2. Data storage 

The acquisition of 2 . 5 seconds (50KHz bandwidth) of fluctuating data 
representing a single 3-wire probe hot wire output requires the rapid 
digitization, processing, display and storage, of 1.5Mb. of data. 150Mb of 
data could easily be collected in 8 hours of transonic wind tunnel 
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testing. The hardware originally purchased with ACQUIRE can adequately 
digitize 14 hot wire channels of the dynamic components. Modifications 
have been made to rapidly transfer the digitized data to the computer, 
adequate hardware exists to transfer the data to permanent storage. But 
only 55Mb of conventional disk space is available for data. A Write Once, 
Read Many (WORM) laser disk drive ($14K) would dramatically improve the 
storage capability, since WORM drives can typically store 600-800Mb of 
data per disc. 

5.3. Compute speed 

5.3.1. Hardware - central processing unit (CPU) 

The real limitation of this system was - and is - in the processing speed 
of the CPU. The original HP 9000/310 CPU was about as fast as an IBM 
PC/AT, and often took minutes to perform a simple evaluation of a few 
thousand points of dynamic data from a single channel. The upgrade to an 
HP 9000/330 (for $13K) in the beginning of 1988 improved the processing 
performance somewhat, but the array processing problem is still not being 
met head on. 

5.3.2. Data structure 

For each computation the array structure requires an indexing algorithm to 
access each sample. Through the indexing mechanism, and by representing 
sample values in an integer format, at least a four- fold savings of 
computer memory and disk space is realized. But the saving of space 
(memory) has become an unnecessary and unacceptable tradeoff. All 
generated data had been simply rescaled (as ratios) existing data, so the 
linear coefficients were easily determined for the new integer data 
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arrays. However, the solution for three unknowns in three simultaneous 
equations does not allow for a simple determination for the linear 
coefficients to scale the integer data. To compute instantaneous 
turbulence fluctuation, each instantaneous voltage must be translated to 
floating point by applying first order coefficients. Then the floating 
point computations (a floating point matrix multiply operation is in 
itself not a fast operation) are accomplished for each instant in time. 

But the answers are in a floating point format and no linear coefficients 
have been determined to convert the floating point answers to a range of 
integer values. Therefore, the hot wire computations are accomplished 
twice - once to determine linear coefficients, and once to store the 
instantaneous turbulence fluctuation in an integer format. Both the 
integer format and the index algorithm produce excessively slow computing 
processes . 

5.3.3. Operating system 

The existing BASIC operating system is a single user, single task 
executive. It cannot support high speed communication via Ethernet. It 
cannot support concurrent operations; program development, data 
acquisition, data processing, and data communications cannot all be 
executing concurrently. 

5.3.4. Programming language 

The BASIC language is an interpreter, rather than a compiler, which trades 
off execution speed for ease of program development and maintenance. 
Although the ACQUIRE software takes advantage of some compiled and 
assembled subroutines - for speed - all of the application software is 
still interpreted. 
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6. RECOMMENDATIONS 

Several solutions exist - they all require much larger investments of time 
and money than a mere doubling of financial resources. 

6.1. Improvement options 

6.1.1. Array processor hardware 

A dedicated array processor is available from Analogic Corp ($27K) which 
is designed to interface to both the HP 9000/330 and the HP software. 
Modifications to the ACQUIRE Digital Signal Processing software (DSP) , or 
development of a user-provided DSP routine to replace the ACQUIRE DSP 
software (3 man-months, est.) would be required 

6.1.2. Faster CPU 

A larger HP 9000/350 CPU ($29K) would quadruple the processing speed, and 
still be able to run the existing ACQUIRE software. 

A combination of the Analogic array processor and the HP 9000/350 would 
provide the best performance possible - without abandoning the ACQUIRE 
software . 

6.1.3. Utilize UNIX operating system 

Provide the multitasking, multiuser environment necessary to support 
concurrent operations, Ethernet (TCP/IP) communications, a choice of 
programming languages - including interpretive and compiled, and a wider 
marketplace for software and hardware solutions like nine track magnetic 
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tape support, laser printer support, graphics and statistics support, and 
data management support. 

6.1.4. Abandon ACQUIRE 

Abandon the ACQUIRE software system, and actively search for a software 
system that operates in the UNIX environment, has the potential for 
supporting the high speed digitizer, can acquire, process, display, and 
save both the mean data and fluctuating data more rapidly than ACQUIRE. 

6.1.5. Remote processing 

Processing the data elsewhere: ACD, MODCOMP or other larger computer 
resource. The solution is suggested by the existence of other 
computational resources that may be made available, including the tunnel 
computer, and would provide parallel processing of the DDAS data once the 
hot wire sensitivities are available, and once the instantaneous data has 
been acquired and transferred to the other resource. This approach 
assumes that a viable communications link like Ethernet is available. At 
NASA LaRC, this capability is called LaRCNET. Although this link is 
proposed for the East Area of LaRC - where 8'TPT is located, its presence 
is still about 2 years distant. LaRCNET also implies - by its very 
existence - that the ACD computational resources will be in great demand. 
The MODCOMP connection, however, proposes a much closer solution. 

Although not an array processor machine, and not yet capable of 
communicating via LaRCNET, the access via a local Ethernet (to eventually 
be a part of LaRCNET) is scheduled for the forth quarter of 1988. 


6.2. Preferred solution 
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The preferred solution is: abandon ACQUIRE for a UNIX compatible set of 

software, translate existing hot wire programs to the UNIX environment, 
purchase new statistical software, and purchase an array processor and a 
faster CPU (the HP 9000/350). About 1 man-month would be required for 
conversion of the hot wire software, and about 3 man-months would be 
required to integrate all the various software and hardware modules. This 
solution minimizes the engineering integration risks attendant in any 
system of this complexity. 
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Figure 1. DDAS System Block Diagram 
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Figure 2. Plot: hotwire voltage vs. mass flow 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


3FT TEST 934 - DATA REDUCTION - 
TEST 934 

RUN 17 LOG PILE 934REDUCE 

POINT S 

TUNNEL 

CONDITIONS 

.9001725 
9.79489233333 
709.55 
465 . 4 

80 . 3028133333 
058.406598442 
.018230240804 


MEAN C HU 1) 
MEANCHU2 ) 
MEANCHU3) 


Mach 

Reynolds No. 

Pt 

Ps 

Tt 

Ue 1 oc i t y 
Dens i t y 
LOGCRhoU) 


OBSERVATION « II 


LOCAL 

WALL 

PROBE 

CONDITIONS 


708.603433333 

477.^245 

540.302813333 

832 72 01491334 

.0195583700825 

1.18876033942 

4.82403933333 

4.67863933333 

5.40758666667 


LOCAL 

FLOOR 

PROBE 

CONDITIONS 


710.075833333 

475 . 0303 

540.302813333 

839.235011329 

.0185023169039 

1.19112144477 

5.75614 

4 . 93762oo6667 

5.3251 


SCU) (HUD 
SCRho) (HU1 ) 
S( To ) (HUD 
S(U) (HU2) 
S(Rho) ( HU2 ) 
S( To ) (HU2) 
S(U) (HU3 ) 
S ( Rho ) ( HU3 ) 
S( To ) ( HU3 ) 


. 0801095340403 
,246631035249 
-.371606324654 
. 0815796200448 
.22379374977 
-.790702032582 
. 0757205995021 
.231703194152 
-.523546487225 


. 0711717036965 
. 151745706072 
- .217170651926 
. 0114987894803 
. 179632925078 
-.724845767689 
. 0299907458573 
. 1945071901 67 
-.427212036868 


u ‘/U ( rms ) 
p ' /P ( rms ) 
to'/To (rms) 


.0132928317592 

.00400499662361 

.000224062285135 


-9.99999999999E*6 

-9.99999999999E*6 

-9.99999999999£*6 


R(RhoU) 

R ( UT0 ) 
R(RhoTO) 


-.997767275712 

.872062903427 

-.07577096438* 


-9.99999999999E*o 

-9.99999999999E** 

-9.99999999999E*6 


M* /M 
P’/P 


. 0 0930 0624276 07 -9 . 99999999999E*o 

. 007255242 05695 -9 . 99999999999E*6 


Figure 3 . Observation Report 
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POINT 3, 
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Figure 4. Plot: Waveforms - Floor strut 
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Figure 5. Plot: Waveforms * Wall strut 
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Pkt$( 1) <stx> 

(ignored) 


( 2) 

TEST 


( 3) 

RUN 


( 4) 

POINT 


( 5) 

TIME 


( 6) 

DATE 


( 7) 

Ps 

tunnel static pressure (psf) 

( 8) 

Pt 

tunnel total pressure (psf) 

( 9) 

Tt 

Tunnel total temperature (deg F) 

(10) 

Mach number 


(11) 

Reynolds number (per chord foot) 

(12) 

PtSl 

Strut 1 (Wall) total pressure (psf) 

(13) 

PsSl 

static pressure (psf) 

(14) 

TtSl 

total temperature (deg F) 

(15) 

PtS2 

Strut 2 (Floor) total pressure (psf) 

(16) 

PsS2 

static pressure (psf) 

(17) 

TtS2 

total temperature (deg F) 

(18) 

PtS3 

Strut 3 (Unused) 

(19) 

PsS3 


(20) 

TtS3 


(21) 

P1HW1 

Strut 1 Hot wire 1 mean voltage 

(22) 

P1HW2 

2 mean voltage 

(23) 

P1HW3 

3 mean voltage 

(24) 

P1HW4 

4 mean voltage 

(25) 

P2HW1 

Strut 2 Hot wire 1 mean voltage 

(26) 

P2HW2 

2 mean voltage 

(27) 

P2HW3 

3 mean voltage 

(28) 

P3HW1 

Strut 3 Hot wire 1 mean voltage 

(29) 

P4HW1 

Strut 4 Hot wire 1 mean voltage 

(30) 

P5HW1 

Strut 5 Hot wire 1 mean voltage 

(31) 

Kulitel 

Microphone RMS voltage 

(32) 

2 


(33) 

3 


(34) 

4 


(35) 

5 


(36) 

6 


(37) 

HW1-4GAIN 

Gain code representing instrument gain 

(38) 

HW5GAIN 


(39) 

HW6GAIN 

Wires 5, 6, and 7 are on probe 2, 

(40) 

HW7GAIN 

wires 1, 2 and 3 

(41) 

KulitelGAIN 

Gain code representing instrument gain 

(42) 

2 


(43) 

3 


(44) 

4 


(45) 

5 


(46) 

6 


(47) 

<ETX> 

(ignored) 


Table 1. Fluctuating Data File Name Format 
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character 

R 

P 

V 

D 

T 

character 

A 

B 

character 

rr 

character 


Fluctuating Data File Name Format 
1 

Real time digitized data 
Playback (FM tape) digitized data 

Velocity ratio - computed data 
Density ratio - computed data 
Temperature ratio - computed data 

2 

"A" 3 -wire probe (wall strut - test 934) 
"B" 3 -wire probe (floor strut - test 934) 


34 

run number (00-99) 
56 


pp point number (00-99) 


character 


78 


cc channel number (00-07) 


1 

character 90 

ss sequence number - if any assigned 


example: RA171203 


real time, probe A, 
run 17, point 12, 
no sequence number 
assigned 


The ACQUIRE software module M0DUSR2 has been written to utilize the 
fluctuating data file naming conventions described above . 

The ACQUIRE software module MODGEN (provided by Data Laboratories) was 
modified to not automatically add a sequence number in column 9 if not 
necessary to differentiate between two files with the same name (in 
columns one through eight) . 


Table 2. Fluctuating Data File Name Format 


Page 40 


APPENDICES 
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APPENDIX A. Program Listings 

This appendix contains the DDAS program listings for all the 
application specific code to acquire, process, display, store 
and transmit the hot wire data. 

The program is separated into two modules: M0DUSR1 provides 
general utility functions. M0DUSR2 provides all hot wire 
specific functions . 
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26042 Modusrl : SUB Modusrl ( INTEGER Routine , Code , OPTIONAL REAL Rvar,Rvar$) 

26043 ! filename M0DUSR1 

26044 ! issue 1 

26045 ! mod 0 

26046 ! date 01 Oct 1987 

26047 ! mod 1 

26048 ! date 22 Dec 1987 

26049 ! programmer S. CLUKEY, Vigyan Research Assoc. 

26050 ! 

26051 ! This program becomes a part of "ACQUIRE", and provides additional 

26052 ! utility functions. As the need for additional functions increases, 

26053 ! so will the functions implemented in this program. 

26054 ! 

26055 OPTION BASE 1 

26056 COM /Arr/ INTEGER Arrowvar(*) 

26057 COM /Cross/ Crval(*) , INTEGER Crvar(*) ,Crtab(*) , Crmptr 

26058 COM /Curs/ Caddrx(*) , INTEGER Cincr(*) , Cpos (*) , Scursor (*) , Cactive , Cpixmin 
(*) ,Cpixmax(*) .Cflag 

26059 COM /Error/ INTEGER Errf.Errtype 

26060 COM /Inp/ Inpstr$ , INTEGER Ilinex, Iliney , Ilinefd, Insertf , Inptype , Inpstr_p 
os,Inpx,Inpy,Inpfd 

26061 COM /Input/ Tinput$ , INTEGER Kposx,Kposy,Quitcode 

26062 COM /Keys/ Keylab$(*) , INTEGER Keymap(*) .Keymenu, Key incr.Okey type, Okeymen 
u 

26063 COM /Localvar/ Lvar$ 

26064 COM /Mem/ Bsv$ (*), Sw$(*) , Sw(*) , INTEGER X(*) , Isv(*) ,Memlenb(*) ,Memlenu(*) 
,Memstartb(*) ,Memstartu(*) , Nummem , Maxnmem , Memmaxl , Tnmem 

26065 COM /Menut/ Menulab$ ,Mvar$(*) ,Mlit$ ,Mvar(*) , INTEGER Mivar<*) ,Menutab(*) , 
Keycode(*) , Sysmod(*) ,Nummi terns .Nummods .Numkeys .Menuptr .Xindex, Yindex 

26066 COM /Param/ Mval$ ,Mval(*) , INTEGER Mstack(*) ,Mvallist(*) .Mstackptr .Mvalpt 
r .Mvalstrptr 

26067 COM /Rcl/ Rclstr$ , INTEGER Rclstr_ptr(*) .Rclptr .Rclnum 

26068 COM /Screen/ INTEGER Garray(*) ,Ctextx(*) .Ctexty .Ctextn, Ctextw, Ocpos (*) ,0 
smp tr , Os type , Update type 

26069 COM /Scrtab/ Cmd_exec$ , INTEGER Morefl,Smtab(*) , Srnptr , Smnum, Smvptr 

26070 COM /Scrvars/ INTEGER Crtvar(*) .Pwidth, Endline, Nextl, Topline 

26071 COM /State/ Status$ , INTEGER Statx, Staty , Statfd, Statusp , Conffl 

26072 COM /Sysvar/ INTEGER Stype,Schg,Sysrec,Sysinit,Seqrunfl,Sysflags,Prdev,G 
rtype 

26073 COM /Trvars/ Trval(*) , INTEGER Trmem(*) ,Trmembit(*) ,Trcrt(*) ,Tryytr(*) ,Tr 
active , Tron , Troverlay , Trf lags , Tr label , S trace , Numtr 

26074 '. 

26075 ! 

26076 ! 

26077 ! 

26078 COM /Plotl/ INTEGER Titlexcoor , Titleycoor , Titlesize .Ntrace , T_chan(*) , Lty 
pe(*) ,Secondc(*) , Titlep , Titlepc , Ploth, Plname$ , Plotstring$ 

26079 COM /Plot2/ REAL Xstart(*) ,Xend(*) ,Ystart(*) ,Yend(*) ,Xorigin(*) , Yorigin( 
*) ,Xmin(*) ,Xmax(*) ,Ymin(*) ,Ymax(*) ,Xtic(*) ,Ytic(*) 

26080 COM /Plot3/ INTEGER Xlabelp(*) ,Xlabelpc(*) , Ylabelp(*) ,Ylabelpc(*) .Commen 
tp(16,8) , Commentpc(16 , 8) , Comments ize(*) ,Npoint(*) , Labelsize (*) 

26081 COM /Plot4/ INTEGER Comcount(*) .Tabcount (*) , Eventbit(*) , Tabvalp(*) , Tabva 
lpc(*) .REAL Tabvalx(*) ,Comxcoor(*) ,Comycoor(*) 

26082 COM /Plot5/ INTEGER Created, Noofpen, Plotdev, Plgrid(*) , Tlabp(*) ,Tlabpc(*) 
,Cur_trace,Plzref(*) ,Xtype(*) ,Ytype(*) 

26083 ! 

26084 ! 

26085 ! 

26086 Usercom: ! 

26087 COM /Usrl/ Cat_array$ (800) [80] ,Fil_nam$[10] , Fil_grp$ ( 10] ,Sym_tbl$(4) [1] ,Nu 
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m_obs_plotted 

26088 COM /Usrl/ From_disk$ { 10 ] , To_disk$ [ 10 ] , File_grp$ [ 10] , File_nani$ [ 10 ] 

26089 COM /Usr2/ Sfn$,REAL Hwsens(*) , Sensinv(*) ,Mean(*) ,Mean_parara(*) , Enorm(*) , S 
tddev_param(*) , Max_param(*) ,Min_paraa(*) ,Vo_param(*) , Vs_param(*) 

26090 COM /Usr2/ Gain_code _14(*) ,Gain_code_57(*) 

26091 COM /Usr2/ C_names$(*) , P_c(*) 

26092 COM /Usr2/ Log_fn$ ,Data_set_title$ , Logged_var_name$(*) , Obs_rec(*) , ?kt_sc , P 
kt$(*) , Pkt_avg(*) ,Num_avgs ,Max_vars ,Rcvd_vars 

26093 COM /Usr2/ Subfile__names$(*) , Subfile_chartst(*) , Initial_obs , Ending_obs , Num 
_obs_recd , Nua_obs_printed , Max_obs_rec , Tag_pkt$ (*) 

26094 COM /Usr2/ Initial_probe , Ending_probe 

26095 COM /Usr2/ A_fn$ , A_set_title$ , A_var_name$(*) , A_rec(*) 

26096 COM /Usr2/ A_subfile_names$(*) , A_sub_chartst(*) 

26097 COM /Usr2/ B_fn$ ,B_set_title$ ,B_var_name$(*) ,B_rec(*) 

26098 COM /Usr2/ B_subfile_names$(*) ,B_sub_chartst(*) 

26099 COM /Usr2/ C_fn$ ,C_set_title$ ,C_var_name$(*) ,C_rec(*) 

26100 COM /Usr2/ C_subfile_names$(*) ,C_sub_chartst(*) 

26101 COM /Usr2/ Hv_rms(*) 

26102 INTEGER V(3) 

26103 INTEGER I 

26104 DIM V$ [ 30] 

26105 SELECT Routine 

26106 CASE 1 ! init pass ! 

26107 RESTORE Menulist 

26108 LOOP 

26109 READ V$ 

26110 EXIT IF V$-"***" 

26111 READ V(*) 

26112 CALL Chkmitem(V(*) ,V$) 

26113 END LOOP 

26114 CASE 2 ! init pass 2 

26115 RESTORE Keylist 

26116 LOOP 

26117 Errfl-0 


26118 Err type- 1 

26119 READ V(l) 

26120 EXIT IF V(1)<0 

26121 READ V(2) 

26122 CALL Chkkey<V(l) ,V(2)) 

26123 END LOOP 

26124 ! 

26125 Menulist: ! 

26126 ! LABEL, Function, Flagl , Flag2 

26127 DATA "UTILITY" ,4000 , 0, 6 

26128 DATA "PLOT UTILITYS" , 4001 , 0 , 6 

26129 DATA "PLOT EJECT" , 4002 , 0 , 262 

26130 DATA "TAG PLOT" , 4003 , 0 , 262 

26131 DATA "TAG PICTURE" , 4004 , 0 , 262 

26132 DATA "PLOT LOG_E" , 4005 , 0 , 262 

26133 DATA "PICTURE LOG_E" , 4006 , 0 , 262 

26134 DATA "FILE UTILITYS" , 4020, 0 , 6 

26135 DATA "CAT GROUP" ,4021 , 8704, 260 

26136 DATA "PURGE" , 4022 , 0 , 6 

26137 DATA "PURGE GROUP" ,4023 , 8704, 390 

26138 DATA "PURGE FILE" ,4024 , 8704 , 390 

26139 DATA "FILE COPY" ,4025 , 0 , 6 

26140 DATA "FROM DISK" , 4026 , 8704 , 22 

26141 DATA "TO DISK" , 4027 , 8704 , 22 

26142 DATA "FILE GROUP" ,4028 , 8704, 22 

26143 DATA "COPY FILES" , 4029 , 8704 , 394 


ORIGINAL PAGE IS 
OF POOR QUALITY 



26144 

26145 

26146 

26147 

26148 

26149 

26150 

26151 

26152 

26153 

26154 

26155 

26156 

26157 

26158 

26159 

26160 
26161 
26162 

26163 

26164 

26165 

26166 

26167 

26168 

26169 

26170 

26171 

26172 

26173 

26174 

26175 

26176 

26177 

26178 

26179 

26180 
26181 
26182 

26183 

26184 

26185 

26186 

26187 

26188 

26189 

26190 

26191 

26192 

26193 

26194 

26195 

26196 

26197 

26198 

26199 

26200 
26201 
26202 

• 26203 
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DATA "MOVE FILES" , 4030 , 3704 , 394 
DATA "***" 

j 

Keylist: ! 

DATA 0,4000 
DATA 100,4001 
DATA 110,4002 
DATA 120,4003 
DATA 130,4004 
DATA 140,4005 
DATA 150,4006 
DATA 200,4020 
DATA 210,4021 
DATA 220,4022 
DATA 221,4023 
DATA 222,4024 
DATA 240,4025 
DATA 241,4026 
DATA 242,4027 
DATA 243,4028 
DATA 245,4030 
DATA 247,4029 
DATA -1,-1 

i 

CASE 3 !RUN TIME INITIALIZATION 
Fil_nam$-" " 

Fil_grp$-"" 

RESTORE Symbols 
READ Sym_tbl$(*) 

Symbols: ! 

DATA "+" 

DATA "X" 

DATA "o" 

DATA 

i 

CASE 4 ! Power on initializations: 

Num_obs_plotted-0 

i 

CASE ELSE 
Usercode: ! 

SELECT Code 

CASE 4002 '.PLOT EJECT 

SELECT Routine 
CASE 31 

OUTPUT Plotdev; "PG" 

Num_obs_plotted-0 ! RESET OBSERVATIONS POINTER for Function 4005 
Num_obs_printed-0 ! RESET OBSERVATIONS PRINTED for Function 4107 
END SELECT 

CASE 4003,4004 !PL0T TAG PICTURE TAG 

SELECT Routine 
CASE 31 

Errf-0 

IF Errf-0 AND Code-4003 THEN PLOTTER IS Plotdev, "HPGL" 

IF Errf-0 THEN 
OFF TIMEOUT 7 
Yminoff-0 
Ymoff-1 
X_gdu_min-0 
X_gdu_max-100*RAT 10 
Y_gdu_min-0 
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26204 

26205 

26206 

26207 

26208 

26209 

26210 
26211 

max* Ymo f f +Ym ino f f 


Y_gdu_max-100 

DEG 

PEN 1 

IF Titlesize<-.8 THEN Titlesize-4 
CSIZE . 65*Titlesize*Ymoff , . 65*RATIO*. 5 
LORG 3 
LDIR 0 

VI EWPORT X_gdu_min , X_gdu_max , Y_gdu_min*Ymo f f+Ymino f f , Y_gdu 


26212 

26213 

26214 

26215 

26216 

26217 

26218 

26219 

26220 
26221 
26222 

26223 

26224 

26225 

26226 

26227 

26228 

26229 

26230 

26231 

26232 

26233 

26234 

26235 

26236 

26237 ! 

26238 ! 


WINDOW X_gdu_min , X_gdu_max , Y_gdu_min , Y_gdu_max 
MOVE X_gdu_max- ( . 20*X_gdu_max) , Y_gdu_max- 10 
FOR 1-1 TO 10 
SELECT I 

CASE 1,2,5 !TEST, RUN, DATE only 
) ##### LABEL Logged_var_name$ ( I ) &" ”&Tag_pkt$(I) 

END SELECT 
NEXT I 

FOR Probe-Initial_probe TO Ending_probe 
SELECT Probe 
CASE 1 

PEN 1 

LABEL "WALL PROBE" 

CASE 2 
PEN 2 

LABEL "FLOOR PROBE" 

END SELECT 
NEXT Probe 
END IF 
PENUP 

VI EWPORT X_gdu_min , X_gdu_max , Y_gdu_min , Y_gdu_max 
WINDOW X_gdu_min , X_gdu_max , Y_gdu_min , Y_gdu_max 
MOVE X_gdu_max , Y_gdu_raax 
PLOTTER IS CRT, "INTERNAL" 

END SELECT 


26239 Plot log e : ! 

26240 CASE 4005,4006 !PLOT LOG(E) vs L0G(Rho*U) 


26241 

26242 

26243 

26244 

26245 

26246 

26247 

26248 

26249 

26250 

26251 

26252 

26253 

26254 

26255 

26256 

26257 

26258 

26259 

26260 
26261 
26262 


SELECT Routine 
CASE 31 

IF Num_obs_plotted-0 THEN !Get the axis plotted 

IF Code-4005 THEN OUTPUT Plotdev; "PG” ! Ej ect old plot 
Mvar$(3)-"RhoU" 

CALL Routine(22 , 618) !PLOT NAME-RhoU 

CALL Routine (3 1,6 16) ! LOAD PLOT 

IF Code-4005 THEN 

CALL Routine(31, 615) ! PLOT PICTURE 

ELSE 

CALL Routine(31 , 612) ! REDRAW PICTURE 
END IF 
END IF 

IF Code-4005 THEN 

PLOTTER IS Plotdev, "HPGL" 

OFF TIMEOUT 7 

ELSE 

PLOTTER IS CRT," INTERNAL" 

END IF 

Yminoff-0 

Ymoff-1 

IF TitlesizeO. 8 THEN Titlesize-4 
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26263 

26264 

26265 ! 

26266 

26267 

26268 

26269 

26270 

26271 

26272 

26273 

26274 

26275 

26276 

26277 


CSIZE . 65*Titlesize*Ymoff , . 65*RATIO*.5 

LORG 5 ! Center the symbol 

FOR Probe-Initial_probe TO Ending_probe 
PEN Probe 

FOR This_obs— Initial_obs TO Ending_obs 
SELECT Probe 
CASE 1 

Num_wires-4 
CASE 2 

Num_wires-3 
CASE ELSE 

Num_wires-1 
END SELECT 

VIEWPORT Xstart(l) ,Xend(l) , Ystart(l)*Ymoff+Yminoff ,Yend 


( 1 ) *Ymof f+Yminof f 


26278 

26279 

26280 
26281 
26282 

26283 

26284 

26285 

26286 

26287 

26288 

26289 

26290 

26291 

26292 

26293 

26294 

26295 

26296 

26297 

26298 

26299 

26300 

26301 

26302 

26303 

26304 

26305 

26306 

26307 

26308 

26309 

26310 

26311 

26312 

26313 

26314 

26315 

26316 

26317 

26318 

26319 

26320 

26321 


WINDOW Xmin(l) ,Xmax(l) ,Ymin(l) ,Ymax(l) 

!###### FOR Wire-1 TO Nunt_wires 
FOR Wire-1 TO 1 
SELECT Probe 
CASE 1 

Pressure-A_rec (4 , This_obs ) 

X_val-MAX( - E6 , A_rec ( 13 , This_obs ) ) 

Y_va 1-MAX ( - E6 , A_rec ( 21+Wire , This_obs ) ) 
CASE 2 

Pressure-B_rec (4 , This_obs ) 

X_va 1-MAX ( - E6 , B_rec ( 13 , This_obs ) ) 
Y_val-MAX( - E6 , B_rec ( 2 1+Wire , This_obs ) ) 
END SELECT 
MOVE X_val,Y_val 
SELECT Pressure 
CASE 700. TO 720. 

LABEL Sym_tbl$(l) 

CASE 850. TO 880. 

LABEL Sym_tbl$(2) 

CASE 1400. TO 1500. 

LABEL Sym_tbl$(3) 

CASE 1700. TO 1800. 

LABEL Sym_tbl$(4) 

END SELECT 
NEXT Wire 
NEXT This_obs 
IF Probe-1 THEN 
MOVE 1. , .6639 
DRAW 1.54, .77832 
LINE TYPE 4 
PEN Probe+3 
MOVE 1. , .67889 
DRAW 1.58, .72076 
LINE TYPE 1 
END IF 

NEXT Probe 

• PLOTTER IS CRT, "INTERNAL" 

END SELECT 

CASE 4021 ! CAT A SELECT GROUP OF FILES 

SELECT Routine 
CASE 21 

Mvar $ ( 3 ) -Fi l_grp$ 

CASE 22 

Fil_grp$-Mvar$ ( 3 ) 
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26322 

26323 

26324 

26325 

26326 

26327 

26328 

26329 

26330 

26331 

26332 

26333 

26334 

26335 

26336 

26337 

26338 

26339 

26340 

26341 

26342 

26343 

26344 

26345 

26346 

26347 

26348 

26349 

26350 

26351 

26352 

26353 

26354 

26355 

26356 

26357 

26358 

26359 

26360 

26361 

26362 

26363 

26364 

26365 

26366 

26367 

26368 

26369 

26370 

26371 

26372 

26373 

26374 

26375 

26376 

26377 

26378 

26379 

26380 

26381 


CASE 31 

CAT; SELECT Fil_grp$ 

END SELECT 

CASE 4023 ! PURGE A SELECT GROUP OF FILES 

SELECT Routine 
CASE 21 

Mvar$ ( 3 ) -Fil_grp$ 

CASE 22 

Fil_grp$-Mvar$ ( 3 ) 

CASE 31 

CAT TO Cat_array$(*) ; SELECT Fil_grp$,NO HEADER, COUNT Num_in_grp 

Fil_grp$-"" 

FOR 1-1 TO Num_in_grp 

Fil nam$-Cat_array$(I) [ 1 ; 10 ] 

PURGE Fil_nam$ 

NEXT I 
END SELECT 

CASE 4024 ! PURGE A FILE 

SELECT Routine 
CASE 21 

Mvar $ ( 3 ) -Fi l_nam$ 

CASE 22 

Fil nam$“Mvar$ ( 3 ) 

CASE 31 

IF LEN ( F i l_nam$ ) >0 THEN PURGE Fil_nam$ 

CALL Pline(0 , "FILE "&Fil_nam$&" HAS BEEN PURGED") 

END SELECT 

CASE 4026 'FROM DISK 

SELECT Routine 
CASE 21 

Mvar$ ( 3 ) -From_disk$ 

CASE 22 

From_disk$-Mvar$ ( 3 ) 

END SELECT 

CASE 4027 !TO DISK 

SELECT Routine 
CASE 21 

Mvar$ ( 3 ) -To_disk$ 

CASE 22 

To_disk$-Mvar$(3) 

END SELECT 

CASE 4028 .'FILE GROUP 

SELECT Routine 
CASE 21 

Mvar$ ( 3 ) -File_grp$ 

CASE 22 

File_grp$-Mvar$ ( 3 ) 

END SELECT 

CASE 4029,4030 !COPY OR MOVE A SELECT GROUP OF FILES 

!( Copy leaves the files on both from and to disks.) 

! ( Move leaves the files only on the to disk, by purging 
! the files successfully copied from the 'from' disk.) 
SELECT Routine 
CASE 21 

Mvar$ ( 3 )-File_grp$ 

CASE 22 

File_grp$-Mvar$(3) 

CASE 31 

CALL Pline(0, "Selecting files. Please wait....’’) 

CAT From_disk$ TO Cat_array$(*) ; SELECT File_grp$,NO HEADER, COUN 
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T Num_in_grp 

26382 ! 

26383 

26384 

26385 

26386 

26387 

26388 

26389 

26390 

26391 

26392 

26393 

26394 

26395 

26396 

26397 ! 

26398 

26399 

26400 


PRINTER IS PRT; WIDTH 108 

PRINT CHR$ ( 12 ) ! Form Feed 

Num_in_grp ; " files have been selected for copying. 


" FROM" ;TAB(30) ; "TO" 
From_disk$ ;TAB(30) ;To_disk$ 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

FOR 1-1 TO Num_in_grp 

PRINT TAB (20) ; Cat_array$(I) [1,10] 

NEXT I 
PRINT 
PRINT 
PRINT 

PRINTER IS CRT 

Numbe r_c op ie d-0 
PRINTER IS PRT 

PRINT "The following files "&CHR$( 27 )&"&dD"&" 


HAVE BEEN COPIED 


&CHR$(27)&"&d<a"&"to "&To disk$ 


26401 

26402 

26403 

26404 

26405 

26406 

26407 

HAS BEEN COPIED FROM 

26408 

26409 

26410 

26411 

N COPIED FROM DISK 

26412 


26413 

26414 
ied. 

26415 

26416 

26417 
ied. 

26418 

26419 

26420 

26421 

26422 

"&CHR$(128)&" 

26423 ! 

26424 ! 

26425 

26426 

26427 

26428 
ONALS" 

26429 

26430 

26431 

26432 

26433 


PRINT 

PRINTER IS CRT 
FOR 1-1 TO Num_in_grp 
ON ERROR GOTO 26414 
File_nam$-Cat_array$(I) [ 1 ; 10 ] 

COPY File_nam$&From_disk$ TO File_nam$&To_disk$ 

CALL Pline(0,"FILE"&CHR$(129)&" "&File_nam$&" "&CHR$(128)&" 
DISK "&From_disk$&" TO DISK "&To_disk$) 

PRINTER IS PRT; WIDTH 108 
Numb e r_c op i e d-Numb e r_c op i e d+ 1 
Cat_array$ (Number_copied) [1,10] -File_nam$ 

PRINT "FILE"&CHR$(129)&" "&File_nam$&" "&CHR$(128)&"HAS BEE 
&From_disk$&" TO DISK "&To_disk$ 

PRINTER IS CRT 
GOTO 26420 

PRINT TAB (10) ; "File "&File_name$&" was NOT successfully cop 


PRINT ERRM$ 
PRINTER IS PRT 
PRINT TAB(IO) ; ' 

PRINT ERRM$ 
PRINTER IS CRT 
NEXT I 

ON ERROR CALL Error 
CALL Pline(0 , " 


File "&File_name$&" was NOT successfully cop 


"&CHR$(129)&" FILE COPYS COMPLETE 


") 


IF Code-4030 THEN ! If a MOVE, then purge the origional file 
PRINTER IS PRT; WIDTH 108 
PRINT 

PRINT Number_copied; " FILES HAVE BEEN COPIED; PURGING ORIGI 

ON ERROR GOTO 26432 
FOR 1-1 TO Number_copied 

PURGE Cat_array$ ( I ) [ 1 , 10 ] &From_disk$ 

NEXT I 

ON ERROR CALL Error 
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26434 

26435 

26436 

26437 

26438 

26439 

26440 

26441 

26442 

26443 


PRINT " MOVE FILES function is complete" 
PRINT 

PRINTER IS CRT 
END IF 
END SELECT 
END SELECT 
END SELECT 
SUBEND 

j 
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24000 Modus r 2 : SUB Modusr2( INTEGER Routine , Code , OPTIONAL REAL R 
var ,Rvar$) 


24002 

24004 

24006 

24007 

24008 

24012 

24013 

24014 

24016 

24017 

24018 

24019 

24020 

24021 

24022 

24023 

24024 

24025 ! ! 

24026 

24027 

24028 

24030 

24031 ! ! 

24032 

24033 

24034 

24035 !! 

24036 

24037 

24038 

24039 

24040 

24041 

24042 

24043 

24044 

24045 

24046 

24047 

24048 

24049 

24050 

24051 

24052 

24053 

24054 

24055 

24056 

24057 

24058 

24059 

24061 

24062 

24063 

24064 

24065 

24066 

24067 


filename MODUSR2 
issue 1 
mod 1 

date 28 Sept 1987 
mod 6 

date 8 Apr 1988 

programmer S. CLUKEY, Vigyan Research Assoc. 

This routine logs mean value data received from another CPU 

via the HP- IB bus. This bus is SC-8, and the primary address is 01. 

This end is NOT System Controller! 

The variables are logged in a disc file that contains a header 
record, a variable names record, and 100 observation records. 

It is an ASCII file. 


This routine also computes coefficients of calibration for multiple 
3 wire probes thru the use of multiple linear regression. 
Alternatively, the coefficients can be entered either by reading 
a coefficient file, or by manually keying in the coefficients - 

With the calibration coefficients, this routine can generate 
sensitivities for the designated probe for each observation 
previously logged. 


This routine can then apply the sensitivities to dynamic HOTWIRE data! 
which is automatically loaded: 


in trace 

1 HOTWIRE 1 

2 HOTWIRE 2 

3 HOTWIRE 3 


and returns: 

fluctuating velocity (u) in trace 4 
fluctuating density (p) in trace 5 
fluctuating temperature (To) in trace 6 


These calculations are performed only between cursor positions 
-if - both cursors of trace 1 are active; or from the 
only cursor to the end of the memory. 

Otherwise, sensitivity coefficients are applied to the entire sample. 


then 


[where XXX is 
the hotwire data disk 


Traces 1, 2, and 3 are automatically loaded from disk files that were 
previously recorded using a naming convention defined here: 

If the log filename is of the format "XXX" 

the test number] 

file names would have the format: 
"MPXXYYZZnn” 

[where M is the mode as follows: 

R - Real time digitization 
P - Playback digitization 
P is the probe selection as follows: 

A - probe A 
B - probe B 
XX is the RUN number 
YY is the POINT number 

of the data most recently logged, 
and therefore most likely to define 
related mean conditions, and will 
contain calculated sensitivities, 
gains of the fluctuating voltages, etc 
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ZZ is the digitizer channel number (this 
naming convention assumes a 1-to-l 
relationship : 

probe "A" 

chan 1 — hotwire 1 
2 2 

3 3 

probe "B" 

5 1 

6 2 

7 3 

nn is the "serial number" automatically 

applied (unfortunately) by ACQUIRE, 
and must be ignored. 

NOTE: the hotwire data disk file name is defined by ACQUIRE function 

"FILENAME". 


Using the sensitivities, mean voltages, and gains from the "A" or "B" 
probe file, 

24092 ! traces 4, 5, and 6 are scaled and created to represent the ratios: 

24093 ! trace 4 — velocity fluctuations / mean velocity 

24094 ! trace 5 - density fluctuations / mean density 

24095 ! trace 6 - temperature fluctuations / mean temperature 

24096 ! 

24097 ! The mean value is removed from these ratios, and an rms value of 

24098 ! the three ratios is stored as variables 47, 48, and 49 of the 

24099 ! appropriate probe file - "A" or "B". 

24100 ! 

24101 ! ! ! Additionally, the ratios above can be retrieved and used to compute 

24102 ! correlations between velocity, density, and temperature fluctuatio 
ns, 

24103 ! and then, using these correlations, go on and compute massflow and 

24104 ! pressure fluctuations. These results are stored much as the ratio 

s 

24105 ! above are stored. 

24106 ! 

24107 OPTION BASE 1 

24137 COM /Arr/ INTEGER Arrowvar(*) 

24138 COM /Cross/ Crval(*) , INTEGER Crvar(*) ,Crtab(*) .Crmptr 

24139 COM /Curs/ Caddrx(*) , INTEGER Cincr (*) , Cpos (*) , Scursor(*) , Cactive , Cpixmin 
(*) ,Cpixmax(*) ,Cflag 

24140 COM /Error/ INTEGER Errf, Err type 

24141 COM /Genvar/ Filename$ (*) , Discdev$ , INTEGER Discmap(*) .Recordnum 

24143 COM /Inp/ Inpstr$ , INTEGER Ilinex, Iliney, Ilinefd, Insertf , Inptype , Inpstr_p 
os , Inpx , Inpy , Inpf d 

24144 COM /Input/ Tinput$ , INTEGER Kposx.Kposy.Quitcode 

24145 COM /Keys/ Key lab$ (*), INTEGER Keymap(*> .Keymenu.Keyincr , Okeytype , Okeymen 
u 

24146 COM /Localvar/ Lvar$ 

24147 COM /Mem/ Bsw$(*) , Sw$(*) , Sw(*) , INTEGER X(*) , Isw(*) ,Memlenb(*) ,Memlenu(*) 
,Memstartb(*) ,Memstartu(*) ,Nummem,Maxnmem,Memmaxl .Tnmern 

24148 COM /Menut/ Menulab$ ,Mvar$(*) ,Mlit$ ,Mvar(*) , INTEGER Mivar(*) ,Menutab(*) , 
Keycode(*) , Sysmod(*) , Nurami terns , Nummods , Numkeys .Menuptr .Xindex, Yindex 

24149 COM /Param/ Mval$ ,Mval(*) , INTEGER Mstack(*) ,Mvallist(*) .Mstackptr .Mvalpt 
r.Mvalstrptr 

24150 COM /Rcl/ Rclstr$ , INTEGER Rclstr_ptr(*) .Rclptr.Rclnum 

24151 COM /Screen/ INTEGER Garray(*) ,Ctextx(*) ,Ctexty ,Ctextn,Ctextw,0cpos(*) ,0 
smptr.Ostype .Update type 


24070 

24071 

24072 

24073 

24074 

24075 

24076 

24077 

24078 

24079 

24080 

24081 

24082 

24083 

24084 

24085 

24086 

24090 

24091 ! ! 
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24152 COM /Scrtab/ Cmd_exec$ , INTEGER Morefl , Smtab(*) , Smptr , Smnum , Smvp tr 

24153 COM /Scrvars/ INTEGER Crtvar(*) , Pwidth , Endline , Nextl , Topline 

24154 COM /StaCe/ Status$ , INTEGER Statx, Staty , Statfd, Statusp , Conf fl 

24155 COM /Sysvar/ INTEGER Stype , Schg, Sysrec , Sysinit , Seqrunfl , Sysflags , Prdev, G 
rtype 

24156 COM /Trvars/ Trval (*), INTEGER Trmem(*) ,Trmembit(*) ,Trcrt(*) ,Tryytr(*) ,Tr 
active , Tron , Troverlay , Trf lags , Tr label , S trace , Numtr 

24157 ! 

24158 ! 

24159 ! 

24160 Usercom:COM /Usr2/ Sfn$ [ 20 ] , REAL Hwsens(3 , 3) , Sensinv(3 , 3) ,Mean(3) ,Mean_par 
am(3) , Enorm(3) , Stddev_param(3) ,Max_param(6) ,Min_param(6) , Vo_param(6) ,Vs_param(6) 

24161 COM /Usr2/ Gain_code_14(17) ,Gain_code_57(8) 

24162 COM /Usr2/ C_names$ (10) [ 10 ] , P_c(3 , 10) ! P_c is the coefficient file 

24163 COM /Usr2/ Log_fn$[20] ,Data_set_title$[80] ,Logged_var_name$(50) [10] ,Obs_re 
c(50, 300) , Pkt_sc , Pkt$(47) [80] , Pkt_avg(45) ,Num_avgs ,Max_vars ,Rcvd_vars 

24164 COM /Usr2/ Subfile_names$(20) [10] , Subfile_chartst(20) , Initial_obs , Ending_o 
bs ,Num_obs_recd,Num_obs_printed,Max_obs_rec ,Tag_pkt$(10) [ 80] 

24165 COM /Usr2/ Initial_probe , Ending_probe 

24166 ! 

24167 ! 

24168 ! The A, B, and C files below contain the "calculated" variables 

24169 ! related to: A) wall strut, B) floor strut, and C) Kulites 

24170 ! (and 'other' 'big end' wires) 

24171 ! 

24172 COM /Usr2/ A_fn$ [ 20 ] , A_set_title$ [80] , A_var_name$ (50) [ 10] , A_rec(50 , 300) 

24173 COM /Usr2/ A_subfile_names$(20) [10] ,A_sub_chartst(20) 

24174 ! 

24175 COM /Usr2/ B_fn$ [20] , B_set_title$ [80] ,B_var_name$(50) [ 10] , B_rec(50 , 300) 

24176 COM /Usr2/ B_subfile_names$(20) [ 10] ,B_sub_chartst(20) 

24177 ! 

24178 COM /Usr2/ C_fn$[20] ,C_set_title$[80] ,C_var_name$(50) [10] ,C_rec(50, 300) 

24179 COM /Usr2/ C_subfile_names$(20) [10] ,C_sub_chartst(20) 

24180 ! 

24181 COM /Usr2/ Hw_nns(3) 

24182 ! 

24183 INTEGER Xch(6) ,Xpt(6) ,Nch(6) ,Npt(6) ,Xchmax,Xptmax, Ymax,Xchmin,Xptmin, Yrain, 
J ,Wtr,Wmem(6) ,Xsch(3) ,Xspt(3) , Ip, Wire 

24184 INTEGER Routinl 

24185 REAL Temp,Wmark(6 , 2) , Vo (6) , Vs (6) , Sum_param(3) ,Sumsq_param(3) ,Vdp(3) ,Wire_g 
ain(4) 

24186 INTEGER V(3) 

24187 DIM V$ [ 30 ] , File_comments$ [ 80 ] 

24188 Routinl-Routine 

24189 SELECT Routine 

24190 CASE 1 ! init pass ! 

24191 RESTORE Menulist 

24192 LOOP 

24193 READ V$ 

24194 EXIT IF V$-"***" 

24195 READ V(*) 

24196 CALL Chkmitem(V(*) ,V$) 

24197 END LOOP 

24198 CASE 2 ! init pass 2 

24199 RESTORE Keylist 

24200 LOOP 

24201 Errfl-0 

24202 Err type-1 

24203 READ V(l) 

24204 EXIT IF V(1)<0 
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24205 READ V(2) 

24206 CALL Chkkey(V(l) ,V(2)) 

24207 END LOOP 

24208 ! 

24209 Menulist: ! 

24210 ! LABEL, Function, Flagl,Flag2 

24211 DATA "HOTWIRE MENU" ,4100, 0, 6 

24212 DATA "COEF FILENAME" ,4101 , 8704, 22 

24213 DATA "LOAD COEFS" ,4102,0,262 

24214 DATA "ENTER COEFS" ,4103 , 0 , 388 

24215 DATA "CALC VEL etc" ,4104,0, 262 

24216 DATA "LOG FILENAME" ,4105 , 8704, 22 

24217 DATA "LOAD LOGFILE" ,4106,0,262 

24218 DATA "PRNT LOGFILE" ,4107,0,262 

24219 DATA "LOG DATA POINT" ,4108 , 0, 390 

24220 DATA "LOG DATA" ,4109,0,6 

24221 DATA "SAMPLES TO AVG" ,4110, 24608 , 22 

24222 DATA "COMPUTE SENS 4111 , 0, 390 

24223 DATA "INITIAL OBS" ,4112 , 24608 , 22 

24224 DATA "ENDING OBS" ,4113 , 24608 , 22 

24225 DATA "HOTWIRE CALC" ,4120,0, 6 

24226 DATA "GET COEF" ,4122,0,6 

24227 DATA "COMPUTE COEFS" ,4123 , 0, 390 

24228 DATA "STORE COEFS" ,4124,0,262 

24229 DATA "CODE TO GAINS" ,4125 , 0, 388 

24230 DATA "FILE TRANSFERS" ,4126,0,6 

24231 DATA "LOGFILE TO PC" ,4127, 0, 388 

24232 DATA "Remake Probe" ,4128,0,388 

24233 DATA "SELECTOR" ,4129 , 0 , 6 

24234 DATA "INITIAL PROBE" ,4130,24608,22 

24235 DATA "ENDING PROBE" ,4131,24608 , 22 

24236 DATA "COMPUTE R etc" ,4132,0,388 

24237 DATA "EDIT COEFS" ,4133 ,0, 388 

24238 DATA "***" 

24239 ! 

24240 Keylist: ! 


24241 

DATA 0,4100 

! HOTWIRE MENU 

24242 

DATA 300,4129 

! SELECTOR 

24243 

DATA 320,4112 

! INITIAL OBS 

24244 

DATA 330,4113 

! ENDING OBS 

24245 

DATA 340,4130 

[INITIAL PROBE 

24246 

DATA 350,4131 

[ENDING PROBE 

24247 

DATA 370,4110 

[SAMPLES TO AVG 

24248 

DATA 400,4120 

[HOTWIRE CALC 

24249 

DATA 420,4122 

[GET COEFS 

24250 

DATA 421,4101 

[COEF FILENAME 

24251 

DATA 422,4102 

[LOAD COEFS 

24252 

DATA 423,4103 

[ENTER COEFS 

24253 

DATA 424,4133 

! EDIT COEFS 

24254 

DATA 425,4123 

[COMPUTE COEFS 

24255 

DATA 427,4124 

[STORE COEFS 

24256 

DATA 430,4111 

[COMPUTE SENS 

24257 

DATA 440,4104 

[CALC VEL etc 

24258 

DATA 450,4132 

[COMPUTE E etc 

24259 

DATA 460,4125 

'.XLATE GAIN CODES 

24260 

DATA 470,4128 

[Remake Probe Data 

24261 

DATA 500,4109 

[LOG DATA 

24262 

DATA 510,4105 

[LOG FILENAME 

24263 

DATA 520,4106 

[LOAD LOGFILE 

24264 

DATA 530,4107 

! PRNT LOGFILE 
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24265 

DATA 

570.4108 

! LOG DATA POINT 

24266 

DATA 

600,4126 

! FILE TRANSFERS 

24267 

DATA 

610,4127 

! LOG FILE TO PC 

24268 

DATA 

-1,-1 


24269 

! 



24270 

CASE 3 

! RUN TIME 

VARIABLE INITIALIZATION 

24271 

RESTORE Coef names 

24272 

READ 

C_names$(*) 


24273 

Coef_names: ! 

THESE ARE THE ORDER IN WHICH THE COEFFICIENTS ARE PROC 

ESSED 




24274 

DATA 

"CONSTANT” 


24275 

DATA 

"L(U) " 


24276 

DATA 

"L(Rho) " 


24277 

DATA 

"L(T0) " 


24278 

DATA 

"L(R)L(T0) " 


24279 

DATA 

”L(R)L(U) ” 


24280 

DATA 

"L(U)L(T0) " 


24281 

DATA 

"LULRLTO" 


24282 

DATA 

" (L(U) )**2" 


24283 

DATA 

"unused" 



24284 ! 

24285 RESTORE Code_14 

24286 READ Gain_code_14(*) ! GAIN CODE CONVERSION TO GAINS FOR HOTWIRES 1-4 

24287 Code_14: ! 

24288 DATA 0 ! CODE 0 off; not defined 

24289 DATA .25 ! 1 

24290 DATA .5 ! 2 

24291 DATA 1. 

24292 DATA 1 . 99 

24293 DATA 3.98 

24294 DATA 7.84 

24295 DATA 15.8 

24296 DATA 31.6 

24297 DATA 63. 

24298 DATA 125. 

24299 DATA 251. 

24300 DATA 501. 

24301 DATA 1000. 

24302 DATA 1995. 

24303 DATA 3981. 

24304 DATA 7943. 

24305 RESTORE Code_57 

24306 READ Gain_code_57(*) !GAIN CODE CONVERSION TO GAINS FOR HOTWIRES 5-7 

24307 Code_57 : » 

24308 DATA 0 ! CODE 0 

24309 DATA 1. ! CODE 1 

24310 DATA 2. 

24311 DATA 5. 

24312 DATA 10. 

24313 DATA 20. 

24314 DATA 50. 

24315 DATA 100. 

24316 Max_obs_rec-300 

24317 Max_vars-50 

24318 Rcvd_vars-45 

24319 CASE 4 ! POWER ON VARIABLE INITIALIZATION 

24320 ! 

24321 Sfn$-"C2" ! COEFFICIENT FILENAME 

24322 Initial_probe-l 

' 24323 Ending_probe-2 
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24324 MAT P_c- (0) 

24325 Num_avgs-5 

24326 Log_fn$-"LQGFILE" 

24327 lnitial_obs-0 

24328 Ending_obs-0 

24329 Num_obs_printed-0 

24330 MAT Pkt$- ("0") 

24331 CASE ELSE 

24332 Usercode:! 

24333 SELECT Code 

24334 ! 

24335 ! 

24336 

24337 

24338 

24339 

24340 

24341 

24342 

24343 ! 

24344 ! 

24345 

24346 

24347 

24348 

24349 

24350 

24351 

24352 

24353 
*Log(T0) 

24354 


CASE 4101 

SELECT Rout ini 
CASE 21 

Mvar$(3)-Sfn$ 
CASE 22 

Sfn$-Mvar$(3) 
END SELECT 


! enter HW COEFFICIENT FILE NAME 


'.get old data, set up parameters 
Ipresent setting 
'.store new data 


! COMPUTE SENSitivities 


*Log(Rho) +A7 


CASE 4111 

SELECT Routine 
CASE 31 

FOR Probe-Initial_probe TO Ending_probe 
FOR R-Initial_obs TO Endlng_obs 
FOR W-l TO 3 

SELECT Probe 
CASE 1 

! S(U) -A2 +A6 

+A8 *Log(Rho)Log(T0)+A9*2.*Log(U) 

A_rec ( 29+V, R) -P_c ( W , 2 ) +P_c (W , 6 ) *A_rec ( 11 , R) +P_c 
( W , 7 ) *A_rec ( 12 , R) +P_c (W , 8 ) *A_r ec ( 14 , R)+P_c (W , 9 ) *2 . *A_r ec ( 10 , R) 

24355 ! S(Rho) -A3 +A5 *Log(T0) +A6 

*Log(U) +A8 *Log(U)Log(T0) 

24356 A_rec ( 3 3+W , R) -P_c ( W , 3 ) +P_c ( W , 5 ) *A_rec ( 12 , R) +P_c 
( W , 6 ) *A_rec ( 10 , R) +P_c ( W , 8 ) *A_r ec ( 16 , R) 

24357 ! S(T0) -A4 +A5 *Log(Rho) +A7 

*Log(U) +A8 *Log(U)Log(Rho) 

24358 A_rec(37+W,R)-P_c(W,4)+P_c(W, 5)*A_rec(ll ,R)+P__c 
( W , 7 ) *A_rec ( 10 , R) +P_c ( W , 8 ) *A_rec ( 1 5 , R) 

24359 CASE 2 

24360 ! S(U) -A2 +A6 *Log(Rho) +A7 

*Log(T0) +A8 *Log(Rho)Log(T0)+A9*2.*Log(U) 

24361 B_rec ( 29+W , R)-P_c ( W , 2 ) +P_c (W , 6 ) *B_rec ( 11 , R) +P_c 
( W , 7 ) *B_rec ( 12 , R) +P_c ( W , 8 ) *B_rec ( 14 , R) +P_c ( W , 9 ) *2 . *B_rec ( 10 , R ) 

24362 ! S(Rho) -A3 +A5 *Log(T0) +A6 

*Log(U) +A8 *Log(U)Log(T0) 

24363 B_rec ( 33+W , R)-P_c (W , 3 ) +P_c ( W , 5 ) *B_rec ( 12 , R) +P_c 
( W , 6 ) *B_rec ( 10 , R) +P_c ( W , 8 ) *B_rec ( 16 , R) 

24364 ! S(T0) -A4 +A5 *Log(Rho) +A7 

*Log(U) +A8 *Log(U)Log(Rho) 

24365 B_rec ( 37+W , R)-P_c ( W , 4) +p_c (W , 5 ) *B_rec ( 11 , R) +P_c 
( W , 7 ) *B_rec ( 10 , R) +P_c ( W , 8 ) *B_rec ( 15 , R) 

24366 END SELECT 

24367 NEXT W 

24368 NEXT R 

24369 NEXT Probe 

24370 GOSUB Log_vars '. Update the disk files 

24371 END SELECT 
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24372 

24373 

24374 

24375 

24376 

24377 

24378 

24379 

24380 

24381 

24382 

24383 

24384 

24385 

24386 

24387 

24388 

24389 

24390 

24391 

24392 

24393 

24394 

24395 

24396 

24397 

24398 

24399 

24400 

24401 

24402 

24403 

24404 

24405 

24406 

24407 

24408 

24409 

24410 

24411 

24412 

24413 

24414 


CASE 4123 ! COMPUTE COEFficienCs 

SELECT Routinl 
CASE 31 

CALL Mir ! Perform Multiple Linear Regression 

! for the probes selected 

END SELECT 

CASE 4124 ! STORE COEFficient file 

SELECT Routinl 
CASE 31 
Errf-0 

CREATE BDAT Sfn$, 1,256 
ASSIGN @Disk TO Sfn$ 

IF Errf-0 THEN OUTPUT @Disk;P_c(*) 

ASSIGN @Disk TO * 

END SELECT 

CASE 4102 ! LOAD COEFficient file 

SELECT Routinl 
CASE 31 
Errf-0 

ASSIGN @Disk TO Sfn$ 

IF Errf-0 THEN ENTER (§Disk; P_c(*) 

ASSIGN @Disk TO * 

PRINTER IS PRT 
PRINT CHR$(12) 

PRINT "THESE ARE THE COEFFICIENTS FROM COEFFICIENT FILE ";Sfn$ 

PRINT 

PRINT 

PRINT 

PRINT TAB(5) ; "WIRE 1" ;TAB(25) ; "WIRE 2" ;TAB(45) ; "WIRE 3" 

PRINT 

FOR 1-1 TO 10 

PRINT P_c (1,1); TAB ( 20 ) ; P_c (2,1) ;TAB(40) ; P_c(3 , 1) 

NEXT I 

PRINTER IS CRT 
END SELECT 

CASE 4103,4133 ! ENTER COEF file; EDIT COEF file 

SELECT Routinl 
CASE 31 

IF Code-4133 THEN GOTO Edit_coefs 
FOR 1-1 TO 3 
BEEP 

FOR J-l TO 10 

CALL Pline(0, "Enter the "&C_names$(J)&" (A"&VAL$ (J)&" ) 


coefficient for wire "&VAL$(I)&": 




24415 ON ERROR RECOVER 24418 

24416 Temp$-FNInput$ ( 3 ) 

24417 P_c(I,J)-VAL(Temp$) 

24418 ON ERROR CALL Error 

24419 NEXT J 

24420 NEXT I 

24421 Edit_coefs: ! 

24422 REPEAT 

24423 CALL Pline(2 , "wirel wire2 wire3") 

24424 FOR Coef_num-l TO 10 

24425 CALL Pline(2+Coef_num, VAL$ (P_c(l , Coef_num) )&" 

" &VAL$ ( P_c ( 2 , Coe f_num) ) &" " &VAL$ ( P_c ( 3 , Coe f~num) ) ) 

24426 NEXT Coef_num 

24427 CALL Pline(17 , "ARE THESE COEFFICIENTS CORRECT FOR PROBE "&V 


AL$(Initial_probe)&"? ") 

24428 Ans$-FNInput$ ( 2 , "NO" ) 
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24429 

") 


CALL Pline(17 , " 


24430 

24431 

24432 

24433 

ts need changing: 

24434 

24435 

24436 

is co be entered: 


SELECT Ans$ 

CASE "YES" , "YE" , "Y" 

CASE ELSE 

CALL Pline(15 , "Enter the wire # (1-3) whose coeffficien 

") 

I-VAL(FNInput$(3 , "1") ) 

IF Id OR I>3 THEN GOTO 24433 

CALL Pline ( 15 , "Enter the coefficient number (A#) whcih 

") 


24437 

24438 

24439 


coefficient for wire "&VAL$(I)&": 


J-VAL ( FNInput$ ( 3 , " 1 n ) ) 

IF Jd OR J>10 THEN GOTO 24436 

CALL Pline(15, "Enter the "&C_names$(J)&" (A"&VAL$(J)&") 


’) 


24440 

24441 

24442 

24443 

24444 

))&" 

24445 

24446 

24447 

24448 

24449 ! 

24450 ! 

24451 ! 

24452 ! 

24453 

24454 

24455 

24456 

24457 

24458 

24459 

24460 

24461 

24462 

24463 

24464 

24465 
s_obs)+l) 


ON ERROR RECOVER 24443 
Temp$-FNInput$ ( 2 ) 

P_c ( I , J ) -VAL(Temp$ ) 

ON ERROR CALL Error 

CALL Pline ( 2+J , VAL$ (P_c ( 1 , J ) ) &" "&VAL$ ( P_c ( 2 , J 

"&VAL$ ( P_c ( 3 , J ) ) ) 

END SELECT 

UNTIL Ans$-"Y" OR Ans$-"YES" OR Ans$-"YE'' 
Schg-BINIOR(Schg,2 A 6+2 A 7) ! CLEAR SCREEN+CLEAR GRAPHICS 
END SELECT 


CASE 4125 1XLATE GAIN CODES TO GAINS 

SELECT Routinl 
CASE 31 

FOR This_obs-Initial_obs TO Ending_obs 

FOR Probe-Initial_probe TO Ending_probe 
! GET MEAN VALUE & GAIN FROM THE LOG FILE 
SELECT Probe 
CASE 1 

FOR Wire-1 TO 4 

! #### SELECT Obs_rec(36,This_obs) 

SELECT 1 
CASE 1 TO 16 

!##### Wire_gain(Wire)-Gain_code_14(Obs_rec(36 ,Thi 


24466 

24467 

24468 

24469 

24470 

24471 

24472 

24473 

24474 

24475 

24476 

24477 

24478 

24479 

e,This_obs)+l) 


Wire_gain(Wire)-Obs_rec(36 ,This_obs) 
A_rec(41+Wire ,This_obs)-Wire_gain(Wire) 
A_var_name $ ( 4 1+W i r e ) - " GAIN " &VAL$ (Wire) 

CASE ELSE 

Obs_rec ( 36 , This_obs ) -0 ! DEFAULT GAIN 
A_var_name $ ( 4 1+Wi r e ) - " GAIN " &VAL$ (Wire) 

END SELECT 
NEXT Wire 
CASE 2 

FOR Wire-1 TO 3 

!#### SELECT Obs_rec(36+Wire ,This_obs) 

SELECT 1 
CASE 1 TO 7 

!##### Wire_gain(Wire)-Gain_code_57(Obs_rec(36+Wir 


24480 

24481 


Wire_gain(Wire)-Obs_rec(36+Wire ,This_obs) 
B_rec (41+Wire , This_obs ) -Wire_gain(Wire ) 
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24482 

24483 

24484 

24485 

24486 

24487 

24488 

24489 

24490 

24491 

24492 

24493 ! 

24494 ! 

24495 

24496 

24497 

24498 

24499 
vation 


B_var_name$ (41+Wire )— "GAIN "&VAL$ (Wire) 
CASE ELSE 

B_rec (41+Wire ,This_obs)-0 ! Default gain - 0 
B_var_name$( 41+Wire) -"GAIN "&VAL$(Wire) 

END SELECT 
NEXT Wire 
END SELECT 
NEXT Probe 
NEXT This_obs 
GOSUB Log_vars 
END SELECT 


CASE 4128 ! REMAKE PROBE DATA 

SELECT Routinl 
CASE 31 

FOR This_obs-Initial_obs TO Ending_obs 

CALL Pline(Staty, "Recalculating Tunnel parameters for obser 
"&VAL$(This obs) ) 


24500 

24501 

24502 

24503 

24504 !#### 


GOSUB Tun_vars 

FOR Probe-Initial_probe TO Ending_probe 
SELECT Probe 
CASE 1 

CALL Pline(Staty, "Recalculating probe A for observa 


tion "&VAL$(This_obs)&" 

24505 

24506 

24507 !#### 

tion "&VAL$(This_obs)&" 


") 

GOSUB Wall_vars 
CASE 2 

CALL Pline(Staty, "Recalculating probe B for observa 
") 


24508 

24509 

24510 

24511 


GOSUB Floor_vars 
END SELECT 
NEXT Probe 
CALL Pline(Staty , " 


") 


24512 

24513 

24514 

24515 

24516 

24517 

24518 

24519 

24520 

24521 

24522 

24523 

24524 

24525 

24526 

24527 

24528 

24529 

24530 

24531 

24532 

24533 

24534 

24535 

24536 

24537 


NEXT This_obs 
GOSUB Log_vars 
END SELECT 


CASE 4132 ! COMPUTE R [CORRELATIONS] etc 

SELECT Routinl 
CASE 31 

! CLEAR ALL DATA TRANSFERS TO START 

Mivar(14)-0 
Mvar(2)-0 

CALL Routine(22 , 239) ! DATA FILE MAP 

! LOAD ONLY CHANNELS 2,3,4 

FOR 1-2 TO 4 
Mvar(2)-L 
Mivar(14)-I 

CALL Routine(22 , 239) ! DATA FILE MAP 

NEXT I 

FOR Probe-Ini tial_probe TO Ending_probe 

FOR This_obs-Initial_obs TO Ending_obs 

DISP "PROBE " &VAL$ ( Probe )&"; OBS "&VAL$(This_obs)&" ; 
! SPECIFY FILENAMES: 

FOR 1-2 TO 4 

Filename$(I)-" 

SELECT I 
CASE 2 


tt . 
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24538 

24539 

24540 

24541 

24542 

24543 

24544 

24545 

24546 

24547 
tc 

24548 

24549 

24550 
! Run 

24551 

24552 

24553 

24554 

24555 

24556 

! Point 


Filename$ (I) [1, 1 ] -"V" 

Filename$ ( I) [7 , 8]-" 04 11 
CASE 3 

Filename$(I) [1, 1]-"D" 

Filename$(I) [7,8]-"05" 

CASE 4 

Filename$(I) [1, 1]-"T" 

Filename$(I) [ 7 , 8 ]-"06" 

END SELECT 

Filename$(I) [ 2 , 2] -CHR$ (64+Probe) ! "A" for Probe 1, e 

SELECT INT(0bs_rec(2 , This_obs) ) 

CASE 0 TO 9 

Filename$ (I) (3 ,4] -"0"&VAL$ (0bs_rec ( 2 , This_obs ) ) 
CASE 10 TO 99 

Filename$(I) [ 3 , 4]-VAL$ (0bs_rec(2 , This_obs) ) 

END SELECT 

SELECT INT(0bs_rec(3 ,This_obs) ) 

CASE 0 TO 9 

Filename$ ( I ) [ 5 , 6 ] -"0"&VAL$ (Obs_rec ( 3 , This_obs ) ) 


24557 CASE 10 TO 99 

24558 Filename$ (I) [ 5 , 6 ] -VAL$ (0bs_rec ( 3 , This_obs ) ) 

24559 END SELECT 

24566 Mivar(14)-I 

24567 CALL Routine(21 , 232) .'DISPLAY THE FILENAME 

24568 NEXT I 

24569 ! LOAD THE THREE CHANNELS INTO MEMORYS 2,3,4 

24570 ! (2-V 3-D 4-T) : 

24571 CALL Routine(31 , 235) ! LOAD DATA 

24572 ! 

24573 FOR Xp-1 TO 3 

24574 DISP "CORRELATION "6.VAL$ (Xp)&" 

24575 Mivar(14)-1 

24576 SELECT Xp 

24577 CASE 1 

24578 ! COPY MEMORY 3 TO MEMORY 1, leaving D in 1 and V in 2: 

24579 Mvar(2)-3 

24580 CALL Routine(22 , 280) IMEMORY COPY V TO 1 

24581 CASE 2 

24582 ! COPY MEMORY 4 TO MEMORY 1, leaving T in 1 and V in 2: 

24583 Mvar(2)-4 

24584 CALL Routine (22 , 280) ! MEMORY COPY T TO 1 

24585 CASE 3 

24586 ! COPY MEMORY 3 TO MEMORY 1, and 

24587 ! COPY MEMORY 4 TO MEMORY 2, leaving D in 1 and T in 2 : 

24588 Mvar(2)-3 

24589 CALL Routine(22 , 280) ! MEMORY COPY D TO 1 

24590 Mivar(14)-2 

24591 Mvar(2)-4 

24592 CALL Routine(22 , 280) IMEMORY COPY T TO 2 

24593 END SELECT 

24594 ! MULTIPLY MEMORY 1 BY MEMORY 2: 

24595 Mivar(14)-0 

24596 CALL Routine(31 , 533) '.MEMORY MULTIPLY 

24597 ! FIND MEAN OF MEMORY 1 

24598 Mivar(14)-1 

24599 CALL Routine (31, 511) ! MEAN 

24600 Temp_mean-Mvar(3) 
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24601 ! STORE MEAN 

24602 

24603 

24604 

24605 

24606 

bs) ! (d'/D)*(v'/V) 

24607 

24608 

bs) ! (v'/V)*(t'/T) 

24609 

24610 

bs) ! (d'/D)*(t'/T) 

24611 

24612 

24613 

24614 

24615 

24616 

bs) ! (d'/D)*(v'/V) 

24617 

24618 

bs)!(v'/V)*(t'/T) 

24619 

24620 

bs) ! (d'/D)*(t'/T) 


INTO DATA BASE AS D*V, V*T, or D*T CORRELATION 
SELECT Probe 
CASE 1 

SELECT Xp 
CASE 1 ! D*V 

Temp_rms-A_rec(48 ,This_obs)*A_rec(47 ,This_o 
CASE 2 ! V*T 

Temp_rms-A_rec(47 ,This_obs)*A_rec(49 ,This_o 
CASE 3 ! D*T 

Temp_rms-A_rec(48 ,This_obs)*A_rec(49 ,This_o 
END SELECT 

A_rec(25+Xp , This_obs )-Temp_mean/Temp_rms 
CASE 2 

SELECT Xp 
CASE 1!D*V 

Temp_rms-B_rec(48 ,This_obs)*B_rec(47 ,This_o 
CASE 2 ! V*T 

Temp_rms-B_rec(47 ,This_obs)*B_rec(49 , This_o 
CASE 3!D*T 

Temp_nns— B_rec (48 , This_obs)*B_rec(49 , This_o 


24621 

24622 

24623 

24624 ! 

24625 

24626 

24627 

24628 

24629 

24630 

24631 

24632 

24633 

24634 

24635 

24636 

24637 

24638 

24639 Etc4132 : ! 

24640 

24641 

24642 

24643 

24644 

24645 

24646 

24647 

24648 

24649 

24650 

24651 

24652 

24653 

24654 


END SELECT 

B_rec ( 25+Xp , This_obs ) -Temp_mean/Temp_ms 
END SELECT 

SELECT Probe 
CASE 1 

A_var_name $(26)-"R( RhoU ) " 

A_var_name$ ( 27 ) -”R(UT0 ) ■ 

A_var_name $ ( 2 8 ) - " R ( RhoTO ) " 

CASE 2 

B_var_name$(26)-"R(R h ° u ) " 

B_var_name $ ( 2 7 ) R (UTO ) - 
B_var_name $ ( 2 8 ) - " R ( RhoTO ) ” 

END SELECT 
GOSUB Log_vars 
NEXT Xp 
NEXT This_obs 
NEXT Probe 

COMPUTE H'/M, P'/P 

FOR Probe-Ini tial_probe TO Endin gp robe 

FOR This_obs-Initial_obs TO Ending_obs 
Mach-Obs_rec(9 ,This_obs) 

M2-Mach*Mach 
SELECT Probe 
CASE 1 

U-A_rec ( 47 , This_obs ) 

Rho-A_rec (48 , This_obs ) 

T0-A_rec (49 , This_obs ) 

R_urho-A_rec ( 26 , This_obs ) 

R_crho-A_rec ( 28 , This_obs ) 

R ut-A_rec ( 27 , This_obs ) 

CASE 2 

U-B_rec(47 ,This_obs) 

Rho-B rec(48 ,This_obs) 
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24655 

24656 

24657 

24658 

24659 

24660 

24661 

24662 

24663 

24664 

24665 

24666 

24667 

24668 

24669 

24670 

24671 

24672 

24673 

24674 

24675 

24676 

24677 

24678 

24679 

24680 

24681 

24682 

24683 

24684 

24685 

24686 

24687 ! 

24688 ! 

24689 

24690 ! 

24691 

24692 

24693 

24694 

24695 

24696 

24697 

ac a time 


T0-B_rec(49 ,This_obs) 

R_urho-B_rec ( 26 , This_obs ) 

R_trho-B_rec(28 ,This_obs) 

R_ut-B_rec ( 27 , This_obs ) 

END SELECT 

IMASSFL is M'/M 

Mass f 1-SQR (U*U+2*R_urho*U*Rho+Rho*Rho ) 

! PRESS is P'/P 

Press-Rho*Rho+M2*M2*( 1 . 44*T0*T0+ . 16*U*U) 
Press-Press+M2*(2 . 4*R_trho*Rho*T0 - . 8*R_urho*U*Rho) 
Press-Press+M2*M2*( - . 96)*R_ut*U*TO 
Press-SQR( Press) 

SELECT Probe 
CASE 1 

A_rec(41,This_obs)-Massfl 
A_rec(46 ,This_obs)-Press 
CASE 2 

B_rec(41,This_obs)-Massfl 
B_rec(46 ,This_obs)-Press 
END SELECT 
GOSUB Log_vars 
NEXT This_obs 
SELECT Probe 
CASE 1 

A_var_name$(41)--M'/M" 

A_var_name$ ( 46 ) - " P ' /P " 

CASE 2 

B_var_name$ (41 )-"M ' /M" 

B_var_name$ (46)-"P'/P" 

END SELECT 
NEXT Probe 
END SELECT 


! CALC VEL etc 

CASE 4104 ! Perform translation from voltages to 

! velocity, density, temperature 

SELECT Routinl 

CASE 31 .'perform action 

Idiag-0 

FOR Probe-Ini tial_probe TO Ending_probe 

FOR This_obs-Initial_obs TO Ending_obs! Do one observation 


24698! 

24699! 

24700! 

24701! 

24702! 

24703! 

24704! 

24705! 

24706 

24707 

24708 
24709! 
24710! 

24711 

24712 
LE 


FOR ALL SAMPLES IN THE MEMORIES - AS DEFINED BY THE CURSORS OF 
OF TRACE 1, CALCULATE THE EQUIVALENT INSTANTANEOUS : 

VELOCITY, DENSITY, TEMPERATURE 


Mivar(14) — 1 
Mvar(2)-0 

CALL Routine(21,239) !Turn off all data files (maps) 


FOR Wire-1 TO 3 

SELECT Probe! GET MEAN VALUE & GAIN FROM THE LOG FI 
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24713 


CASE 1 

24714 


Mean (Wire ) -A_r ec ( 17+Wire , This_obs ) 

24715 


Wire gain(Wire)-A rec (41+Wire , This obs) 

24716 


CASE 2 

24717 


Mean(Wire)-B_rec(17+Wire , Inieial_obs) 

24718 


Wire gain(Wire)-B rec (41+Wire .This obs) 

24719 


END SELECT 

24720 


FOR Sens-1 TO 3! GET THE SENSITVITYS FROM THE LOG F 

ILE 

24721 


SELECT Probe 

24722 


CASE 1 

24723 


Hwsens(Wire , Sens)-A rec(25+(Sens*4)+Wire , Th 

is obs) 

24724 


CASE 2 

24725 


Hwsens(Wire , Sens)-B rec(25+(Sens*4)+Wire , Th 

is obs) 

24726 


CASE 3 

24727 


Hwsens(Wire , Sens)-C rec(25+(Sens*4)+Wire ,Th 

is_obs) 

24728 


END SELECT 

24729 


NEXT Sens 

24730* 

24731* 


SET UP DATA FILE MAP FOR TRACES 1,2,3 ONLY 

24732 


Mivar(14)-Wire 

24733 


Mvar(2)-1 

24734 


CALL Routine(22 , 239) !TURN ON daca files l,2,or3 

247-35 


Mivar ( 14 ) -3+Wire 

24736 


Mvar ( 2 ) -0 

24737 


CALL Routine (22, 239 )! TURN OFF data files 4,5,or6 

24738! 

24739 


Temp$[l]-"R" 

24740 


Temp$ [ 2 ] -CHR$ ( 64+Probe ) 

24741 


SELECT I NT (Obs rec (2 ,This_obs) ) 

24742 


CASE 0 TO 9 

24743 


Temp$ [ 3 , 4 ] -" 0 " &VAL$ ( INT (Obs rec(2,This obs)))'.R 

UN 

24744 


CASE 10 TO 99 

24745 


Temp$ [3,4] -VAL$ ( INT ( Obs_rec ( 2 , This_obs ) ) ) ! RUN 

24746 


END SELECT 

24747 


SELECT INT(0bs_rec(3 , This_obs) ) 

24748 


CASE 0 TO 9 

24749 


Temp$ [5,6] -"0"&VAL$ ( INT (Obs rec(3,This obs)))!P 

OINT 

24750 


CASE 10 TO 99 

24751 


Temp$ [5,6] -VAL$ ( INT (Obs_rec ( 3 , This_obs ) ) ) ! POINT 

24752 


END SELECT 

24753 


SELECT Probe 

24754 


CASE 1 

24755 


Temp$ [7,8] -"0"&VAL$ (Wire) 

24756 


CASE 2 

24757 


Temp$ [7 , 8 ]-"0"&VAL$ (Wire+4) 

24758 


END SELECT 

24759 


Temp$ [ 9 , 9 ]-"l" 

24760 

DISP Temp$ 


24761 

WAIT 1 


24762 


Filename$(Wire)-Temp$ 

24763 


Mivar (14) -Wire 

24764 


CALL Routine (21 , 232) ! DISPLAY THE FILENAME 

24765! 

24766! 
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24767 

! TRANSFER CHANNEL CHARACTERISTICS 

24768 

! TO COMPUTED CHANNELS 

24769 

FOR 1-1 TO 15 

24770 

Sw(Wire+3 , I)-Sw(Wire ,1) 

24771 

NEXT I 

24772 

NEXT Wire 

24773! 


24774 

CALL Routine (43, 235) ! LOAD DATA 

24775 

CALL Routine(31, 235) 

24776 

FOR 1-1 TO 3 ! For each channel 

24777 

! CLEAR MEMORY TAGS FOR "REMOVE MEAN" FUNCTION 

24778 

FOR J-0 TO 2 

24779 

Sw$(I , 13+J)-"" 

24780 

NEXT J 

24781 

Mivar(14)-I 

24782 

CALL Routine (31, 524) ! REMOVE MEAN 

24783 

CALL Routine(31,513) ! FIND RMS 

24784 

Hw rms(I)-Mvar(3) 

24785 

NEXT I 

24786! 


24787! 

TRACE 1, 2, AND 3 POINT TO THE MEMORIES 

24788! 

CONTAINING THE FLUCTUATING COMPONENT ONLY 

24789! 


24790 ! 


24791! 

SET UP DATA FILE MAP FOR RESULTS: 4,5,6 ONLY 

24792 

FOR Results-4 TO 6 

24793 

Mivar ( 14) -Resul ts - 3 

24794 

Mvar(2)**0 

24795 

CALL Routine (22, 239 )! TURN OFF data files 1,2,3 

24796 

Mivar (14) -Results 

24797 

Mvar(2)-1 

24798 

CALL Routine (22, 239) !TURN ON data files 4,5,6 

24799! 


24800 

SELECT Results 

24801 

CASE 4 ! Velocity 

24802 

Temp$[l]-"V" 

24803 

CASE 5 ! Density 

24804 

Temp$[l]-"D" 

24805 

CASE 6 ! Temperature 

24806 

Temp$[l]-"T" 

24807 

END SELECT 

24808 

Temp$ [ 2 ] -CHR$ ( 64+Probe ) 

24809 

SELECT INT(0bs_rec(2 ,This_obs) ) 

24810 

CASE 0 TO 9 

24811 

Temp$ [ 3 , 4 ] -" 0 " &VAL$ ( INT ( Obs_rec ( 2 , This_obs ) ) ) ! R 

UN 


24812 

CASE 10 TO 99 

24813 

Temp$ {3,4] -VAL$ ( INT ( Obs_rec ( 2 , This_obs ) ) ) ! RUN 

24814 

END SELECT 

24815 

SELECT INT(Obs_rec(3 ,This_obs) ) 

24816 

CASE 0 TO 9 

24817 

Temp$ [ 5 , 6 ] - " 0 " &VAL$ ( INT ( Ob s_rec ( 3 , This_ob s ) ) ) ! P 

OINT 


24818 

CASE 10 TO 99 

24819 

Temp$ [5,6] -VAL$ ( INT (Ob s_rec ( 3 , This_obs ) ) ) ! POINT 

24820 

END SELECT 

24821 

SELECT Probe 

24822 

CASE 1 

24823 

Temp$ [ 7 , 8 ] - " 0 " &VAL$ (Results) 

24824 

CASE 2 
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24825 

24826 

24827 

24828 

24829 

24830 

24831 

24832 

24833 

24834 
24835! 
24836! 
24837! 

24838 

24839 

24840 

24841 

24842 

24843 

24844 

24845 

24846 

24847 
24848! 
24849! 
24850! 

24851 

24852 

24853 

24854 

24855 


i 

i 

t 


Temp$ [7,3 ] -"0 "&VAL$ (Results+4) 

END SELECT 

Filename$ (Results)-Temp$ 

Mivar(14)-Results 

CALL RouCine(21 , 232) ! DISPLAY THE FILENAME 
NEXT Results 

FILENAMES SET UP TO STORE DATA FILES 
MAT Sensinv- INV(Hwsens) 


DETERMINE SCALING (from a sampling of the data) 
Wtr-1 

Numb_samp-Sw(Wtr ,2) 

SELECT Numb_samp 
CASE <100 

Numb_sub-Numb_s amp 
CASE 100 TO 10000 
Numb_sub— 100 
CASE >10000 

Numb_sub-Numb_samp/100 
END SELECT 

DETERMINE Mean, Standard deviation, Min, Max, range, Offset, Scale 
of velocity, density, and temperature. 


Value-0 

FOR Wtr-1 TO 6! pick up trace number 

Wmem(Wtr)-Trmem(Wtr) !pick up the memory number 
IF (BIT(Cactive,Wtr*2-l) ) THEN ! ACTIVE 1st CURSOR 
Wmark(Wtr,l)-Caddrx(Wtr*2-l) !get the position o 


f first cursor 


24856 
CURSOR 

24857 

24858 

24859 

e last point in memory 

24860 

24861 

24862 
R 

24863 

of second cursor 


IF (BIT(Cactive , Wtr*2) ) THEN ! and ACTIVE 2nd 

Wmark(Wtr , 2)-Caddrx(Wtr*2) 

ELSE ! and IN-ACTIVE 2nd CURSOR 

Wmark(Wtr , 2)-Sw(Wtr , 2) !no second cursor; us 

END IF 

ELSE ! IN- ACTIVE 1st CURSOR 

IF (BIT(Cactive ,Wtr*2) ) THEN ! ACTIVE 2nd CURSO 

Wmark(Wtr , l)-Caddrx(Wtr*2) ! get the position 


24864 

24865 

24866 

24867 

24868 

24869 

position is <- 

24870 

24871 

24872 

24873 

24874 


ELSE 

Wmark(Wtr, l)-0 
END IF 

Wmark(Wtr , 2 )-Sw(Vtr , 2 ) 

END IF 

IF Wmark(Wtr,l)>Wmark(Wtr,2) THEN iassure starting 

Temp-Wmark(Wtr,l) ! to ending position 
Wmark(Wtr,l)-Wmark(Wtr,2) . 

Wmark(Wtr , 2) -Temp 
END IF 

Xch(Wtr)-( (FNMems(l .Wmem(Wtr) )+Wmark(Wtr , 1) ) DIV 10 


24)+l ! starting row address 

24875 Xpt(Wtr)-((FNMems(l,Wmem(Wtr))+Wmark(Wtr,l)) MOD 10 

24)+l! starting column address 

' 24876 Nch(Wtr)-(Wmark(Wtr , 2) -Wmark(Wtr , 1) ) DIV 1024!numbe 
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r of rows 



24877 


Npt (Wtr)-(Wmark(Wtr , 2) -Wmark(Wtr , 1) ) MOD 1024 ! numbe 

r of columns 

in last row 

24878 


Vo(Wtr)-Sw(Trmem(Wtr) ,22)/100!y OFFSET offset 

24879 


Vs (Wtr)-Sw(Trmem(Wtr) , 21) ! y GAIN (VOLTS) sensitivi 

ty 



24880 


NEXT Wtr 

24881 


Strtpt-MAX(Wmark(l , 1) ,Wmark(2,l) ,Wmark(3,l)) 

24882 


Stoppt-MIN(Wmark(l , 2) ,Wmark(2,2) ,Wmark(3,2)) 

24883 ! 



24884 ! FOR 

THE 

SAMPLE POINTS: 

24885 


FOR I-Strtpt TO Stoppt STEP Numb samp/Numb sub 

24886 


FOR Wire-1 TO 3 

24887 


IF FNMems(l,Tnnem(Wire) )+I<-Sw(Wire ,2) THEN 

24888 


Xsch(Wire)-((FNMems(l,Trmem(Wire))+I) DIV 1 

024)+l 



24889 


Xspt(Wire)-((FNMems(l,Trmem(Wire))+I) MOD 1 

024)+l 



24890 


Temp-X(Xsch(Wire) .Xspt(Wire)) 

24891 


Value-( (Temp/65536) -Vo(Wire) )*Vs(Wire) 

24892 


Enorm(Wire)-Value/(Mean(Wire)*Wire_gain(Wir 

e)) 



24893 


IF Idiag-9 AND Wire-1 AND I<900 THEN 

24894 


PRINTER IS PRT 

24895 


PRINT I ; "th Sample - (Wire" ; Wire - X: 

" ; Temp ; " X(volts): 

" ; Value ; " Enorm: " ; Enorm(Wire) 

24896 


PRINTER IS CRT 

24897 


END IF 

24898 


END IF 

24899 


NEXT Wire 

24900 


MAT Vdp- Sensinv*Enorm 

24901 


FOR Ip-1 TO 3 

24902 


SELECT I 

24903 


CASE Strtpt 

24904 


Max_param( Ip+3 ) -Vdp ( Ip ) 

24905 


Min_param( Ip+3 ) -Vdp ( Ip ) 

24906 


Sum_param ( Ip ) -Vdp ( Ip ) 

24907 


Sumsq_param( Ip ) -Vdp ( Ip ) *Vdp ( Ip ) 

24908 


CASE ELSE 

24909 


IF Vdp(Ip)>Max_param(Ip+3) THEN Max_param(I 

p+3)-Vdp(Ip) 



24910 


IF Vdp ( Ip ) <Min_param( Ip+3 ) THEN Min_param(I 

p+3)-Vdp(Ip) 



24911 


Sum_param( Ip ) -Sum_param(Ip ) +Vdp ( Ip ) 

24912 


Sumsq_param(Ip)-Sumsq_param(Ip)+(Vdp(Ip)*Vd 

p(Ip)) 



24913 


IF Idiag AND Ip-3 THEN 

24914 


PRINTER IS PRT 

24915 


PRINT "TEST VALUES n ;Vdp(*) 

24916 


PRINT "MAX ” ; Max param(*) 

24917 


PRINT "MIN " ;Min param(*) 

24918 


PRINTER IS CRT 

24919 


END IF 

24920 


END SELECT 

24921 


NEXT Ip 

24922 


NEXT I 

24923! 



24924! CALCULATE THE SLOPE (Vs_paraa) AND INTERCEPT (Vo_param) 

24925 


FOR Ip-1 TO 3 

24926 


Mean param( Ip ) -Sum_param ( Ip ) /Numb_sub 
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24927 IF ( ( Sumsq_param( Ip) - (Numb_sub*Mean_param( Ip ) *Mean 
par am (Ip) ) )/ (Mean_param(Ip) -1) )>-0 THEN 

24928 Stddev_param(Ip)-SQR( (Sumsq_param( Ip) - (Numb_sub 
*Mean_param(Ip)*Mean_param(Ip) ) )/(Mean_param(Ip) - 1) ) 


24929 


ELSE 

24930 


Stddev param(Ip)-0 

24931 


END IF 

24932 


Vo_param( Ip+3 )-(Max_param( Ip+3 )+Min_par am ( Ip+3 ) )/2 . 

24933 


IF Max_param( Ip+3 ) -Vo_param( Ip+3) OO THEN 

24934 


Vs_param(Ip+3)-10*( (Max_param(Ip+3) -Vo_param( Ip 

+3))/32767) 


24935 


ELSE 

24936 


Vs_param(Ip+3)-10 . 

24937 


END IF 

24938! 

STORE OFFSETS 

AND SLOPES IN THE MEMORIES POINTED TO BY TRACES 4,5,6 

24939 


Sw(Ip+3 , 28)-Vo_param(Ip+3) 

24940 


Sw(Ip+3 , 27)-Vs_param(Ip+3) 

24941! 

IDENTIFY DISPLAY AS USER UNITS 

24942 


Isw(Ip+3 , 13)-1 

24943 


SELECT Ip 

24944 


CASE 1 

24945 


Sw$(Ip+3.H)-"u , /U" 

24946 


CASE 2 

24947 


Sw$(Ip+3 , ll)-"p'/T" 

24948 


CASE 3 

24949 


Sw$(Ip+3, ll)-"t'/T" 

24950 


END SELECT 

24951 


NEXT Ip 

24952! 



24953! 

COMPDT 

E u, p. To 

24954! 



24955 


FOR I-Strtpt TO Stoppt 

24956 


FOR Wire-1 TO 3 

24957 


Temp-X(Xch(Wire) .Xpt(Wire)) 

24958 


Value— ( (Temp/65536 ) -Vo (Wire ) ) *Vs ( Wire ) 

24959 


Enorm(Wire)-Value/(Mean(Wire)*Wire gain(Wire) ) 

24960 


Xpt(Wire)-Xpt(Wire)+l 

24961 


IF Xpt(Wire)>1024 THEN 

24962 


Xpt(Wire)-l 

24963 


Xch(Wire)-Xch(Wire)+l 

24964 


END IF 

24965 


NEXT Wire 

24966! 



24967 


MAT Vdp- Sensinv*Enorm 

24968 


IF Idiag THEN 

24969 


PRINTER IS PRT; WIDTH 134 

24970 


IF I-S Crept THEN 

24971 


PRINT "MAX_P: " ;Max_param(*) 

24972 


PRINT "MIN_P: " ;Min_param(*) 

24973 


PRINT 

24974 


PRINT "SLOPE: " ; Vs_param(*) 

24975 


PRINT "INTERCEPT: " ; Vo_param(*) 

24976 


PRINT 

24977 


PRINT "Means: " ;Mean(*) 

24978 


PRINT "Gains: " ; Wire_gain(*) 

24979 


PRINT 

24980 


END IF 

24981 


PRINTER IS CRT 

24982 


END IF 

24983 


FOR Ip-4 TO 6! PUT u, p. To IN MEMORIES UNDER TRAC 
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ES 4, 5, 6 


24984 

) 

Value-INT( (Vdp(Ip-3) -Vo_param(Ip) )/Vs_param( Ip) 

24985 

IF Idiag AND Ip-6 THEN 

24986 

IF I MOD 100-0 AND I<1 500 THEN 

24987 

PRINTER IS PRT; WIDTH 134 

24988 

PRINT "SAMPLE: ";I; "TRACE: ";Ip,"Vdp: " 

; Vdp(Ip-3) , " Value: 

Value 

24989 

PRINTER IS CRT 

24990 

END IF 

24991 

END IF 

24992 

SELECT Value! Limit the range to the 16 bit int 

eger 


24993 

CASE -32768 TO 32767 

24994 

X(Xch(Ip) ,Xpt( Ip)) -Value 

24995 

CASE >32767 

24996 

X(Xch(Ip) ,Xpt(Ip) )-32767 

24997 

CASE <32768 

24998 

X(Xch(Ip),Xpt(Ip))— 32768 

24999 

END SELECT 

25000 

Xpt(Ip)-Xpt(Ip)+l 

25001 

IF Xpt(Ip)>1024 THEN 

25002 

Xpt(Ip)-l 

25003 

Xch(Ip)-Xch(Ip)+l 

25004 

END IF 

25005 

NEXT Ip 

25006 ! 

DISP I 

25007 

NEXT I 

25008 

FOR Ip-1 TO 3! For each ratio (u'/U,p ' /? , tO'/TO) 

25009 

FOR 1-4 TO 6! CLEAR MEMORY TAGS FOR "REMOVE MEAN” F 

UNCTION 


25010 

FOR J-0 TO 2 

25011 

Sw$(I , 13+J)-"" 

25012 

NEXT J 

25013 

NEXT I 

25014 

Mivar(14)-Ip+3 

25015 

CALL Routine (31, 524) ! REMOVE MEAN 

25016 

CALL Routine(31 , 513) ! FIND RMS 

25017 

Value-Mvar ( 3 ) 

25018 

SELECT Probe 

25019 

CASE 1 

25020 

A rec(46+Ip , This obs)-Value 

25021 

CASE 2 

25022 

B rec(46+Ip,This_obs) -Value 

25023 

END SELECT 

25024 

NEXT Ip 

25025! 


25026! 


25027 

CALL Routine (43, 234) ! STORE FLUCTUATING DATA 

25028 

CALL Routine (31, 234) 

25029! 


25030! 


25031 ! i i i i • i 


25032 

RESTORE A_names 

25033 

READ A_var_name$(*) 

25034 

RESTORE B_names 

25035 

READ B_var_name$(*) 

25036 !!!!!!! 


25037 

GOSUB Log_vars ! Save the computed values on disc 

25038 ! 
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25039 

PRINTER IS PRT 

25040 

PRINT CHR$(12) ! FORM FEED 

25041 

PRINT Data_set_title$ ;TAB(40) ; "PROBE " ; Probe ; TAB (60) ; "0 

BSERVATION * ";This obs 


25042 

PRINT 

25043 

PRINT "TEST " ;0bs rec(l,This obs) 

25044 

PRINT "RUN " ;0bs rec(2,This obs) 

25045 

PRINT "POINT " ;0bs rec(3,This obs) 

25046 

PRINT "P_TOTAL ";0bs rec(7,This obs) 

25047 

PRINT "P STATIC ";0bs rec(6,This obs) 

25048 

PRINT "T TOTAL ";0bs rec(8,This obs) 

25049 

PRINT "MACH # " ;0bs rec(9,This obs) 

25050 

PRINT 

25051 

rt 

PRINT TAB(20) ; "WIRE 1" ;TAB(40) ; "WIRE 2" ; TAB(60) ; "WIRE 3 

25052 

PRINT "Vjnean" ;TAB(20) ;Mean(l) ;TAB(40) ;Mean(2) ;TAB(60) ; 

Mean(3) 


25053 

PRINT "V_rms" ; TAB (20) ;Hw_rms(l) ;TAB(40) ;Hw_rms(2) ;TAB(6 

0) ;Hw rms(3) 


25054 

PRINT "GAIN"; TAB (20) ; Wire_gain(l) ;TAB(40) ; Wire_gain(2) ; 

TAB(60) ;Wire gain(3) 


25055 

PRINT "SENSITIVITY"; TAB (20) ;Hwsens(l , 1) ;TAB(40) ;Hwsens( 

2,1) ; TAB (60) ;Hvsens(3,l) 


25056 

PRINT TAB(20) ;Hwsens(l,2) ;TAB(40) ;Hwsens(2,2) ;TAB(60) ;H 

wsens(3 , 2) 


25057 

PRINT TAB(20) ;Hwsens(l, 3) ;TAB(40) ;Hwsens(2,3) ;TAB(60) ;H 

wsens(3 , 3) 


25058 

SELECT Probe 

25059 

CASE 1 

25060 

Vel_f-A_rec(47 ,This_obs) 

25061 

Dens_f-A_rec(48 ,This_obs) 

25062 

Temp_f-A - rec (49 , This_obs ) 

25063 

CASE 2 

25064 

Vel_f-B_rec(47 ,This_obs) 

25065 

Dens_f-B_rec ( 48 , This_obs ) 

25066 

Temp f-B rec(49 ,This_obs) 

25067 

END SELECT 

25068 

PRINT 

25069 

PRINT "u'/U_rms" ;TAB(20) ; Vel_f 

25070 

PRINT "p'/T rms" ;TAB(20) ; Dens f 

25071 

PRINT "t'/T rms" ; TAB ( 20 ); Temp f 

25072 

PRINT 

25073 

PRINTER IS CRT 


25074 NEXT This_obs 

25075 ! 

25076 ! 

25077 ! 

25078 NEXT Probe 

25079 ! 

25080 END SELECT 

25081 Schg-BINI0R(Schg, 16384) !SET BIT 14 TO SAY WAVEFORMS CHANGED 

25082 ! 

25083 ! 

25084 ! 

25085 Log file menu: ! 

25086 CASE 4105 ! enter LOG FILE NAME 

25087 SELECT Rout ini 

25088 CASE 21 !get old data, set up parameters 

25089 Mvar$(3)-Log_fn$ Ipresent setting 

25090 CASE 22 ! store new data 
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! LOAD LOG FILE 


25091 Log_fn$-Mvar$(3) 

25092 END SELECT 

25093 CASE 4106 

25094 SELECT Routinl 

25095 CASE 31 

25096 ON ERROR GOTO No_file 

25097 ASSIGN @Disk TO Log_fn$ 

25098 ON ERROR CALL Error 

25099 ENTER @Disk, 1 ; Data_set_title$ ,Max_obs_rec ,Max_vars , Logged_var_n 
ame$(*) ,Numsubfile,Subfile_names$(*) , Subfile_chartst(*) 

25100 Num_obs_recd-Subfile_chartst(l) 

25101 ! 

25102 REDIM Obs_rec(Max_vars , Max_obs_rec ) ,A_rec(Max_vars ,Max_obs_rec) 

, B_rec(Max_vars ,Max_obs_rec) 

25103 ! 

25104 ENTER @DIsk,2 

25105 ENTER @Disk;Obs_rec(*) 

25106 ! Set up WALL strut file 

25107 ASSIGN @Diska TO Log_fn$&"A" 

25108 ENTER @Diska, 1 ;A_set_title$ .Dummy , Dummy, A_var_name$(*) , Dummy ,A_ 
subfile_names$(*) , A_sub_chartst(*) 

25109 ENTER @Diska,2 

25110 ENTER @Diska; A_rec(*) 

25111 ! Set up FLOOR strut file 

25112 ASSIGN @Diskb TO Log_fn$&"B" 

25113 ENTER @Diskb, 1 ;B_set_title$ , Dummy .Dummy, B_var_name$ (*) , Dummy , B_ 
subfile_names$(*) ,B_sub_chartst(*) 

25114 ENTER @Diskb,2 

25115 ENTER @Diskb;B_rec(*) 

25116 ! Set up OTHER strut file 

25117 ! ! ! ! ASSIGN @Diskc TO Log_fn$&"C" 

25118 !!!! ENTER @Diskc , 1 ;C_set_title$ , Dummy , Dummy , C_var_name$ (*) , Dummy , C 

_sub f i 1 e_names$ ( * ) , C_s ub_c harts t (* ) 

25119 !!!! ENTER (§Diskc,2 

25120 !*!! ENTER @Diskc ; C_rec(*) 

25121 GOTO Got_it 

25122 No_file: !NO FILE. !OK, SO CREATE THE FILE 

25123 ON ERROR CALL Error 


25124 

25125 

25126 

25127 

25128 

25129 

25130 

25131 !!!! 
0 

25132 ! ! ! ! 
•25133 

25134 

25135 

25136 

25137 

25138 

25139 

25140 

25141 

25142 

25143 

25144 


Errf— O 

CREATE BOAT Log_fn$ , INT(8*Max_vars*Max_obs_rec/1280)+2 , 1280 
ASSIGN @Disk TO Log_fn$ 

CREATE BOAT Log_fn$&"A" , INT(8*Max_vars*Max_obs_rec/1280)+2 , 1280 
ASSIGN @Diska TO Log_fn$&"A" 

CREATE BOAT Log_fn$&"B" . INT(8*Max_vars*Max_obs_rec/1280)+2 , 1280 
ASSIGN @Diskb TO Log_fn$&"B" 

CREATE BDAT Log_fn$&"C" , INT(8*Max_vars*Max_obs_rec/1280)+2 , 128 

ASSIGN @Diskc TO Log_fn$&"C" 

File_conunents$-" " 

DISP "ENTER THE LOG FILE DATA SET TITLE - Return IF NONE"; 

INPUT " " , File_comments$ 

IF LEN(File_comments$)>0 THEN 

Data_set_title$-File_comments$ 

A_set_title$-"WALL STRUT : : : "&File_comments$ 

B_set_title$— "FLOOR STRUT : : : "&File_comments$ 

C_set_title§— "KULITE, MISC ; : : "&File_conments$ 

ELSE 

Data_set_title$-"" 

A_set_title$-"" 

B set title$-"" 


Page 70 


25145 

25146 

25147 

25148 

25149 

25150 

25151 

25152 

25153 

25154 

25155 

25156 

25157 

25158 

25159 

25160 

25161 

25162 

25163 

25164 

25165 

25166 

25167 

25168 

25169 

25170 

25171 

25172 

25173 

25174 

25175 

25176 

25177 

25178 

25179 

25180 

25181 

25182 

25183 

25184 

25185 

25186 

25187 

25188 

25189 

25190 

25191 

25192 

25193 

25194 

25195 

25196 

25197 

25198 

25199 

25200 

25201 

25202 

25203 

25204 


C_set_citle$-" " 

END IF 

Num_obs_recd-0 

Numsub file-0 

MAT Subf ile_names$— ("") 

MAT A_subf ile_names$— ("") 
MAT B_subf ile_names$— ("") 
MAT C_subf ile_names$- ("") 
MAT Subfile_chartst- (0) 

MAT A_sub_chartst- (0) 

MAT B_sub_chartsc- (0) 

MAT C_sub_chartst- (0) 

MAT 0bs_rec- (-9999999.99999) 
MAT A_rec- (-9999999.99999) 
MAT B_rec- (-9999999.99999) 
MAT C_rec- (-9999999.99999) 

! ~ 

! 


RESTORE Var_names ! NOTE: 

! IN FUTURE LINKS, THESE NAMES WILL PRECEDE 
! THE VALUES IN THE DATA PACKET RECEIVED. 
READ Logged_var_name$ ( *) 

Var names:! THE VARIABLE NAMES THAT EACH OBSERVATION RECORD CONTAINS: 


DATA 

"TEST" ! 

1 

DATA 

"RUN" ! 

2 

DATA 

"POINT" ! 

3 

DATA 

"TIME" ! 

4 

DATA 

"DATE" ! 

5 

DATA 

" Ps " ! 

6 

DATA 

"Pt" ! 

7 

DATA 

"Tt" ! 

8 

DATA 

"MACH" ! 

9 

DATA 

"REYNO" ! 10 

DATA 

"PtSl" 


DATA 

"PsSl" 


DATA 

"TtSl" 


DATA 

"PtS2" 


DATA 

"PsS2" 


DATA 

"TtS2" 


DATA 

"PtS3" 


DATA 

”PsS3" 


DATA 

"TtS3" 


DATA 

"P1HW1" 


DATA 

"P1HW2" 


DATA 

"P1HW3" 


DATA 

"P1HW4" 


DATA 

"P2HW1" 


DATA 

"P2HW2" 


DATA 

"P2HW3" 


DATA 

"P3HW1" 


DATA 

"P4HW1" 


DATA 

"P5HW1" 


DATA 

"KULITE1" 

! 30 

DATA 

"KULITE2" 

! 31 

DATA 

"KULITE3" 

! 32 

DATA 

"KULITE4" 

! 33 

DATA 

"KULITE5" 

! 34 

DATA 

"KULITE6" 

! 35 

DATA 

"HW1-4 GAIN" 

DATA 

"P2HW1GAIN 

If 
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25205 DATA "P2HW2GAIN" 

25206 DATA "P2HW3GAIN" 

25207 DATA "K1 GAIN" 

25208 DATA "K2 GAIN" 

25209 DATA "K3 GAIN" 

25210 DATA "K4 GAIN" 

25211 DATA "K5 GAIN" 

25212 DATA "K6 GAIN" 

25213 DATA "" 

25214 DATA "" 

25215 DATA "" 

25216 DATA "" 

25217 DATA "" 

25218 ! 

25219 ! 

25220 Subfile_chartst(l)-Nvun_obs_recd 

25221 OUTPUT @Disk, 1 ; Data_set_title$ ,Max_obs_rec ,Max_vars , Logged_var 
name$(*) .Numsubf lie , Subfile_names$(*) , Subf ile_chartst(*) 

25222 STATUS @Disk, 3 ;Norecs .Nobpr 

25223 CONTROL @Disk, 7; Norecs.Nobpr 

25224 ENTER @Dlsk,2 

25225 OUTPUT @Disk;0bs_rec(*) 

25226 ! 


25227 RESTORE A_names ! NOTE: 

25228 ! THESE ARE THE VARIABLE NAMES FOR THE 

25229 ! WALL STRUT DATA FILE "A” 

25230 READ A_var_name $ ( * ) 

25231 A_names:! THE VARIABLE NAMES THAT EACH "A" DATA RECORD CONTAINS: 


25232 

DATA 

"TEST" ! 1 

25233 

DATA 

"RUN" ! 2 

25234 

DATA 

"POINT" ! 3 

25235 

DATA 

"PtSl" 

25236 

DATA 

"PsSl" 

25237 

DATA 

"TtSl" 

25238 

DATA 

"VELOCITY" 

25239 

DATA 

"DENSITY" 

25240 

DATA 

"Ts" 

25241 

DATA 

"L(U) " 

25242 

DATA 

"L(Rho) " 

25243 

DATA 

"L(T0) " 

25244 

DATA 

"L(RhoU) " 

25245 

DATA 

"L(R)L(T0) " 

25246 

DATA 

"L(R)L(U) " 

25247 

DATA 

"L(U)L(T0) " 

25248 

DATA 

"LULRLTO" 

25249 

DATA 

"P1HW1" 

25250 

DATA 

"P1HW2" 

25251 

DATA 

"P1HW3" 

25252 

DATA 

"P1HW4" 

25253 

DATA 

"LOG(PlHWl) ’ 

25254 

DATA 

"L0G(P1HW2) ’ 

25255 

DATA 

"L0G(P1HW3) * 

25256 

DATA 

"L0G(P1HW4) " 

25257 

DATA 

"R(RhoU) " 

25258 

DATA 

"R(UTO) " 

25259 

DATA 

"R(RhoTO) " 

25260 

DATA 

n n 

25261 

DATA 

"S(U)1" 

25262 

DATA 

"S(U)2” 

25263 

DATA 

"S(U) 3" 


! 18 


! 22 


! 26 


! 30 


PRECEDING PAGE BLANK NOT FILMED 
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25264 

DATA "" 


25265 

DATA "S(Rho) 1" 

! 34 

25266 

DATA "S(Rho)2" 


25267 

DATA "S(Rho) 3" 


25268 

DATA "" 


25269 

DATA "S(T0)1" 

! 38 

25270 

DATA "S(T0)2" 


25271 

DATA "S(T0)3" 


25272 

DATA "M'/M" 

! 41 

25273 

DATA "GAIN 1" 

! 42 

25274 

DATA "GAIN 2" 


25275 

DATA "GAIN 3" 


25276 

DATA "GAIN 4" 


25277 

DATA "P'/P" 

! 46 

25278 

DATA "u'/U" 

! 47 

25279 

DATA "p'/P" 


25280 

DATA "to '/To" 


25281 

DATA "" 


25282 ! 



25283 ! 



25284 

OUTPUT @Diska , 1 ; A 

set title$,Max obs rec.Max vars.A var name$(* 

) .Numsubfile , A 

subfile_names$(*) , Subfile_chartst(*) 

25285 

STATUS @Diska, 3 ;Norecs,Nobpr 

25286 

CONTROL @Diska, 7;Norecs , Nobpr 

25287 

ENTER @Diska,2 


25288 

OUTPUT @Diska;A_rec(*) 

25289 ! 



25290 

RESTORE B_names 

! NOTE : 

25291 


! THESE ARE THE VARIABLE NAMES FOR THE 

25292 


! FLOOR STRUT DATA FILE "B" 

25293 

READ B var name$(- 

*) 

25294 B names:! 

THE VARIABLE NAMES 

THAT EACH "B" DATA RECORD CONTAINS: 

25295 

DATA "TEST" ! 1 


25296 

DATA "RUN" ! 2 


25297 

DATA "POINT" ! 3 


25298 

DATA "PtS2" 


25299 

DATA "PsS2" 


25300 

DATA "TtS2" 


25301 

DATA "VELOCITY" 


25302 

DATA "DENSITY" 


25303 

DATA "Ts" 


25304 

DATA "L(U) " 


25305 

DATA "L(Rho) " 


25306 

DATA "L(T0) " 


25307 

DATA "L(RhoU) " 


25308 

DATA "L(R)L(T0) " 


25309 

DATA ”L(R)L(U) " 


25310 

DATA "L(U)L(T0) " 


25311 

DATA "LULRLT0" 


25312 

DATA "P2HW1" 

! 18 

25313 

DATA "P2HW2" 


25314 

DATA "P2HW3" 


25315 

DATA "" 


25316 

DATA "L0G(P2HW1) " 

>22 

25317 

DATA "LOG(P2HW2) " 


25318 

DATA "LOG(P2HW3) " 


25319 

DATA "" 


25320 

DATA "R(RhoU) " 

>26 

25321 

DATA "R(UT0) " 


25322 

DATA "R(RhoTO) " 
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25323 

DATA "" 


25324 

DATA "S(U)1" 

! 30 

25325 

DATA "S(U)2" 


25326 

DATA ”S(U)3" 


25327 

DATA "" 


25328 

DATA "S (Rho) 1" 

! 34 

25329 

DATA "S(Rho)2" 


25330 

DATA "S (Rho) 3" 


25331 

DATA "" 


25332 

DATA "S (TO) 1" 

! 38 

25333 

DATA "S(T0)2" 


25334 

DATA "S(T0)3" 


25335 

DATA "M'/M" 

! 41 

25336 

DATA "GAIN 1" 

! 42 

25337 

DATA "GAIN 2" 


25338 

DATA "GAIN 3" 


25339 

DATA "" 


25340 

DATA "P'/P" 

! 46 

25341 

DATA "u'/U" 

! 47 

25342 

DATA "p'/P" 


25343 

DATA "to '/To" 


25344 

DATA "" 



25345 ! 

25346 ! 

25347 


OUTPUT @Diskb , 1 ;B_set_title$ ,Max obs rec.Max vars , B var name 

1 <4 r* i r- - ~ — — — 


) .Numsubfile ,B_subfile names$(*) , Subfile chartst(*) 

^ c ^ / o ~ _ _ N ' 


25348 

25349 

25350 

25351 


STATUS @Diskb , 3 ; Norecs , Nobpr 
CONTROL @Diskb , 7 ; Norecs , Nobpr 
ENTER @Diskb,2 


25352 

i 

25353 

! RESTORE C names 

25354 


25355 

25356 

! READ C_var_name$(*) 

25357 

C_names : ! THE VARIABLE NAMES 

25358 

DATA "TEST" ! 1 

25359 

DATA "RUN" ! 2 

25360 

DATA "POINT" ! 3 

25361 

DATA "PtS3" 

25362 

DATA "PsS3" 

25363 

DATA "TtS3" 

25364 

DATA "VELOCITY" 

25365 

DATA "DENSITY" 

25366 

DATA "Ts" 

25367 

DATA "L(U) " 

25368 

DATA "L(Rho) " 

25369 

DATA "L(T0) " 

25370 

DATA "L(RhoU) " 

25371 

DATA "L(R)L(T0) " 

25372 

DATA "L(R)L(U) " 

25373 

DATA "L(U)L(T0) " 

25374 

DATA "LULRLT0" 

25375 

DATA "P3HW1" 

25376 

DATA "P4HW1" 

25377 

DATA "P5HW1" 

25378 

DATA "" 

25379 

DATA "L0G(P3HW1) " 

25380 

DATA "L0G(P4HW1) " 

25381 

DATA "L0G(P5HW1) " 


NOTE: 

THESE ARE THE VARIABLE NAMES FOR THE 
'OTHER' STRUT DATA FILE "C" 

iT EACH "C" DATA RECORD CONTAINS: 


! 22 
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25382 

DATA "" 


25383 

DATA "" 


25384 

DATA "" 


25385 

DATA "" 


25386 

DATA " n 


25387 

DATA "S(U)1" ! 30 


25388 

DATA "S (U) 2" 


25389 

DATA "S(U)3" 


25390 

DATA "" 


25391 

DATA "S (Rho) 1" ! 34 


25392 

DATA "S(Rho)2" 


25393 

DATA "S(Rho)3" 


25394 

DATA "" 


25395 

DATA "S (TO) 1" ! 38 


25396 

DATA "S(T0)2" 


25397 

DATA "S(T0)3" 


25398 

DATA " " 


25399 

DATA ■■ 


25400 

DATA "" 


25401 

DATA "" 


25402 

DATA "■ 


25403 

DATA "" 


25404 

DATA "u'/U" 


25405 

DATA "p'/P" 


25406 

DATA "to '/To" 


25407 

DATA "" 


25408 

j 


25409 



25410 

!!!! OUTPUT @Diskc,l;C set title$,Max obs 

_rec , Max_vars , C_var_name$ ( 

*) .Numsubfile ,C subfile names$(*) .Subfile chartst(*) 


25411 

! ! ! ! STATUS @Diskc , 3 ; Norecs , Nobpr 


25412 

!!!! CONTROL @Diskc, 7; Norecs , Nobpr 


25413 

! ! ! ! ENTER @Diskc , 2 


25414 

!!!! OUTPUT @Diskc ; C_rec(*) 


25415 

j 


25416 

Got it: ! 


25417 

ASSIGN @Disk TO * 


25418 

ASSIGN (§Diska TO * 


25419 

ASSIGN (§Diskb TO * 


25420 

! ! ! ! ASSIGN @Diskc TO * 


25421 

Num_obs_plotted-0 


25422 

Num_ob s_p r i n t e d-0 


25423 

Initial_obs-Num_obs_recd+l 


25424 

Ending obs-Initial obs 


25425 

END SELECT 


25426 

CASE 4107 ! PRINT LOG 

FILE INFO 

25427 

SELECT Routinl 


25428 

CASE 31 


25429 

Errf-0 


25430 

j 


25431 

Errorf-0 


25432 

ASSIGN @Disk TO Log_fn$ 


25433 

IF Errorf-0 THEN ENTER @Disk, 1 ; Data_set_title$ ,Max_obs_rec ,Max_ 


vars , Logged_var_name$(*) .Numsubfile , Subfile_names$(*) , Subfile_chartst(*) 

25434 Num_obs_recd-Subfile_chartst(l) 

25435 ! 

25436 REDIM Obs_rec(Max_vars ,Max_obs_rec) , A_rec(Max_vars ,Max_obs_rec) 
, B_rec (Max_vars , Max_obs_rec ) 

25437 ! 

25438 ENTER @Disk,2 
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25439 IF Errf-0 AND Num_obs_recd>-L AND Num_obs_recd<-Max_obs_rec THE 
N ENTER @Disk;0bs_rec(*) 

25440 ASSIGN @Disk TO * 

25441 ! 

25442 Errorf-0 

25443 ASSIGN @Diska TO Log_fn$&"A" 

25444 IF Errorf-0 THEN ENTER @Diska, 1; A_sec_title$ .Dummy .Dummy ,A_var_ 
name$(*) , Dummy ,A_subfile_names$(*) , A_sub_chartst(*) 

25445 ENTER (§Diska,2 

25446 IF Errf-0 AND Num_obs_recd>— I AND Num_obs_recd<-Max_obs_rec THE 
N ENTER @Diska;A_rec(*) 

25447 ASSIGN @Disk TO * 

25448 ! 

25449 Errorf-0 

25450 ASSIGN @Diskb TO Log_fn$6<"B" 

25451 IF Errorf-0 THEN ENTER @Diskb,l;B_set_title$, Dummy, Dummy, B_var_ 
name$(*) , Dummy, B_subfile_names$ (*) ,B_sub_chartst(*) 

25452 ENTER @Diskb,2 

25453 IF Errf-0 AND Num_obs_recd>— 1 AND Num_obs_recd<-Max_obs_rec THE 
N ENTER @Diskb;B_rec(*) 

25454 ASSIGN @Disk TO * 

25455 ! 

25456 Errorf-0 

25457 ! !!! ASSIGN (§Diskc TO Log_fn$6."C" 

25458 !!!! IF Errorf-0 THEN ENTER @Diskc,l;C_set_title$, Dummy, Dummy, C_var_n 

ame$(*) , Dummy ,C_subfile_names$(*) ,C_sub_chartst(*) 

25459 !!!! ENTER @Diskc, 2 

25460 ! ! ! ! IF Errf-0 AND Num_obs_recd>-l AND Ntun_obs_recd<-Max_obs_rec THEN 

ENTER @Diskc;C_rec(*) 

25461 ! ! ! ! ASSIGN (§Disk TO * 

25462 ! 

25463 GOSUB PrinC_log_data 

25464 GOTO End_4107 

25465 ! 

25466 ! 


25467 Print_log_data: '.Subroutine to print a formatted report. 

25468 ! 

25469 ! 


25470 

25471 

25472 

25473 

25474 

25475 
ted 


IF Num_obs_printed— 0 THEN Num_obs_printed— Initial_obs - 1 
WHILE Num_obs_printed<Ending_obs 

Num_obs_printed-Num_obs_printed+l 

PRINTER IS PRT; WIDTH 108 

PRINT CHR$(12) ! Form feed 

PRINT Data_set_title$;TAB(60) ; "OBSERVATION # " ;Num_obs_prin 


25476 

25477 

25478 

25479 

FILE " ; Log_fn$ 

25480 

L";TAB(69); "LOCAL" 

25481 

25482 

25483 

25484 

25485 

IONS"&Nu_$ ;TAB(85) 

25486 

25487 


IF Num_obs_printed>0 AND Num_obs_printed<-Max_obs_rec THEN 
PRINT 

PRINT "TEST " ;Obs_rec(l,Num_obs_printed) 

PRINT "RUN " ;0bs_rec(2 , Num_obs_printed) ; TAB (18) ; "LOG 

PRINT "POINT " ;0bs_rec(3,Num_obs_printed) ;TAB(46) ; "LOCA 

PRINT TAB (46) ; "WALL" ; TAB (69) ; "FLOOR" 

PRINT TAB(23) ; "TUNNEL" ; TAB (46) ; "PROBE" ; TAB (69) ; "PROBE" 
U_$-CHR$(27)&"&dD" (Underline 
Nu_$-CHR$ (27)6." 6.d@ " ! No unde r 1 ine 

PRINT TAB(23) ;U_$6."CONDITIONS"6^u_$ ;TAB(54) ;U_$6."CONDIT 
; U_$6<" CONDITIONS ■ &Nu_$ 

PRINT 

PRINT "Mach" ;TAB(23) ;0bs_rec(9 ,Num_obs_printed) 
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PRINT "Reynolds No . " ; TAB(23 ) ; Obs_rec ( 10 , Num_obs_pr inted 


25488 

) 

25489 PRINT "Pt" ;TAB(23) ;0bs_rec(7 , Num_obs_orinted) ;TAB(46) ;A 
_rec(4,Num_obs_printed) ;TAB(69) ;B_rec(4,Num_obs_printed) 

25490 PRINT "Ps" ; TAB (23) ;0bs_rec(6 ,Num_obs_pr inted) ; TAB (46) ; A 
_rec(5 ,Num_obs_printed) ; TAB (69) ;B_rec(5 ,Num_obs_pr inted) 

25491 . PRINT "Tt" ;TAB(23) ;0bs_rec(8 , Num_obs_printed) ;TAB(46) ;A 

_rec ( 6 , Num_obs_printed) ; TAB ( 69 ) ; B_r ec ( 6 , Num_obs_pr inted) 

25492 PRINT "Velocity" ;TAB(23) ;0bs_rec(46,Num_obs_printed) ;TA 
B(46) ;A_rec(7 ,Num_obs_printed) ;TAB(69) ; B_rec(7 ,Num_obs_pr inted) 

25493 PRINT "Density" ; TAB(23) ; 0bs_rec (47 ,Num_obs_pr inted) ; TAB 
(46) ;A_rec(8,Num_obs_printed) ;TAB(69) ;B_rec(8 ,Num_obs_printed) 

25494 PRINT "LOG(RhoU) " ;TAB(46) ; A_rec(13 , Num_obs_printed) ; TAB 
(69) ;B_rec(13 ,Num_obs_printed) 

25495 PRINT 

25496 PRINT "MEAN(HWl) " ; TAB (46) ;A_rec (18 ,Num_obs_pr inted) ; TAB 
(69) ;B_rec(18 ,Num_obs_printed) 

25497 PRINT "MEAN(HW2) " ;TAB(46) ; A_rec(19 , Num_obs_printed) ;TAB 
( 69 ) ; B_r ec ( 19 , Num_obs_pr inted) 

25498 PRINT "MEAN(HW3) " ;TAB(46) ;A_rec (20, Num_obs_pr inted) ;TAB 
(69) ; B_rec( 20 ,Num_obs_pr inted) 

25499 PRINT 

25500 PRINT "S(U) (HW1) " ;TAB(46) ;A_rec(30 ,Num_obs_printed) ;T 

AB(69) ;B_rec( 30, Num_obs_pr inted) 

25501 “ PRINT "S(Rho) (HW1) " ;TAB(46) ;A_rec( 34, Num_obs_pr inted) ;T 

AB(69) ;B_rec(34,Num_obs_printed) 

25502 PRINT "S(To) (HW1) " ;TAB(46) ; A_rec(38 ,Num_obs__printed) ;T 

AB(69) ; B_rec(38 ,Num_obs_printed) 

25503 PRINT "S(U) (HW2) " ;TAB(46) ;A_rec( 31. Num_obs_pr inted) ;T 

AB(69) ;B_rec (31, Nuu»_obs_pr inted) 

25504 PRINT "S(Rho) (HW2) " ; TAB (46) ; A_rec(35 , Num_obs_pr inted) ; T 
AB(69) ;B_rec(35 ,Num_obs_pr inted) 

25505 PRINT "S(To) (HW2) " ;TAB(46) ;A_rec(39 ,Nvun_obs_printed) ;T 
AB(69) ; B_rec( 39 ,Num_obs_pr inted) 

25506 PRINT "S(U) (HW3) " ; TAB (46) ; A_rec(32 ,Num_obs_printed) ;T 

AB(69) ; B_rec (32 , Nua_obs_pr inted) 

25507 PRINT "S(Rho) (HW3) " ;TAB(46) ;A_rec(36 ,Num_obs_printed) ;T 
AB(69) ;B_rec(36,Num_obs_printed) 

25508 PRINT "S(To) (HW3) " ;TAB(46) ;A_rec (40, Ntun_obs_pr inted) ;T 
AB(69) ;B_rec( 40, Num_obs_pr inted) 

25509 PRINT 

25510 PRINT "u'/U (rms) " ; TAB(46) ; A_rec(47 ,Nura_obs_printed) ; TA 
B(69) ; B_rec (47 ,Num_obs_pr inted) 

25511 PRINT "p'/P (rms) " ;TAB(46) ;A_rec (48 ,Num_obs_pr inted) ;TA 
B(69) ; B_rec (48 ,Num_obs_pr inted) 

25512 PRINT "to'/To (rms) " ;TAB(46) ; A_rec(49 ,Num_obs_printed) ; 
TAB(69) ;B_rec (49, Num_obs_pr inted) 

25513 PRINT 

25514 PRINT "R(RhoU) " ;TAB(46) ;A_rec(26 ,Num_obs_printed) ;TAB(6 

9) ;B_rec(26,Num_obs_printed) 

25515 PRINT "R(UT0) " ;TAB(46) ; A_rec(27 , Num_obs_printed) ;TAB(69 

) ; B_rec ( 27 , Num_obs_printed) 

25516 PRINT "R(RhoTO)" ;TAB(46) ;A_rec(28 ,Num_obs_printed) ;TAB( 

69) ;B_rec(28 ,Num_obs_printed) 

25517 PRINT 

25518 PRINT "M'/K" ;TAB(46) ; A_rec (41, Num_obs_pr inted) ;TAB(69) ; 

B_rec(41 ,Num_obs_printed) 

25519 PRINT "P • /?" ;TAB(46) ; A_rec(46 ,Num_obs_pr inted) ;TAB(69) ; 
B_rec(46 , Nua_obs_printed) 

25520 ! 
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25521 

25522 

25523 

25524 

25525 

25526 

25527 

25528 

25529 

25530 

25531 

25532 

25533 

25534 

25535 

25536 

25537 

25538 

25539 

25540 

25541 

25542 

25543 

25544 

25545 

25546 

25547 

25548 

25549 

25550 

25551 

25552 

25553 

25554 

25555 

25556 

25557 

25558 

25559 

25560 

25561 

25562 

25563 

25564 

25565 

25566 

25567 

25568 

25569 

25570 

25571 

25572 

25573 

25574 

25575 

25576 

25577 

25578 

25579 

25580 


END IF 

PRINTER IS CRT 
END WHILE 

RETURN 


j 

i 

End_4107 : END SELECT 

i 

j 

CASE 4108 

SELECT Rout ini 
CASE 31 
Load_data: ! 

! THIS ROUTINE IS THE ACQUISITION AND DATA LOGGING PORTION OF THE HOTWIRE 
! DYNAMIC DATA SYSTEM. 


STEVE CLUKEY, VIGYAN RESEARCH ASSOCIATES, SEPTEMBER, 1987 

DATA IS RECEIVED FROM ANOTHER COMPUTER (MODCOMP) VIA HP- IB IN ASCII. 
UPON REQUEST OF THE OPERATOR, THIS DATA IS AVERAGED, AND STORED (LOGGED) 
ON DISK AS AN OBSERVATION IN AN X-Y ARRAY. SOME VALUES ARE CALCULATED 
AND STORED ALONG WITH THE DATA RECEIVED FROM THE OTHER COMPUTER. 

THESE DATA OBSERVATIONS WILL BE USED TO COMPUTE HOTWIRE SENSITIVITIES, 
AND THEREFORE THE FORMAT OF THE DISK FILE SHOULD BE COMPATIBLE WITH 
THE COMPUTATIONAL PROGRAMS. 


FOR EACH OBSERVATION, "Max_vars" VARIABLES ARE STORED: 


! THE PROCESS IS: 

! 1 . CREATE AN OBSERVATION FILE, OR DECLARE AN EXISTING FILE, AND SET UP 
! TO APPEND (WITH ADDITIONAL RECORDS) ADDITIONAL OBSERVATIONS. 

! 2 . WHEN REQUESTED BY THE OPERATOR, DATA IS READ FROM THE OTHER COMPUTER 
! VIA HP- IB. BOTH THE DATA HOTWIRE DATA AND THE TUNNEL PARAMETER DATA 
! ARE AVERAGED, AND STORED IN THE NEXT AVAILABLE RECORD. 

! 3 . WHEN REQUESTED BY THE OPERATOR, THE FILE IS CLOSED. 

! AT SOME LATER TIME, THE FILE WILL BE PROCESSED TO PROVIDE A HOTWIRE 
! SENSITIVITY FILE. 


CALL Pline(0, "OBSERVATION LOGGING STARTED") 

WAIT 1 

Pkt_sc-8 

MAT Pkt avg- (0) 

1-0 

FOR 1-1 TO Num_avgs 
GOSUB Enable_link 
GOSUB Get_packet 

IF 1-1 THEN IDon't average test, run, point. . . 
FOR J-l TO 3 

Pkt__avg(J)-VAL(Pkt$(J+l) )*Num_avgs 
NEXT J 

FOR J-4 TO 5 

Pkt_avg(J)-0! THIS WAS DATE AND TIME 
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25581 NEXT J 

25582 FOR J-l TO 10! SAVE THE FIRST 10 FOR "TAG PLOT , PICTURE" 

25583 Tag_pkt$(J)-Pkt$(J+l) 

25584 NEXT J 

25585 END IF 

25586 FOR J-6 TO Rcvd_vars 

25587 ON ERROR GOTO 25589! Ignore error if not a good VAL 

25588 Pkt_avg ( J ) -Pkt_avg ( J ) +VAL ( Pkt$ ( J+l ) ) 

25589 ON ERROR CALL Error 

25590 NEXT J 

25591 NEXT I 

25592 ! ! ! ! ! GOSUB Disable_link 

25593 MAT Pkt_avg- Pkt_avg/(Num_avgs) ! FIND THE AVERAGE 

25594 ! 

25595 ! Get ready to log data - by opening files 

25596 ASSIGN @Disk TO Log_fn$ 

25597 ASSIGN @Diska TO Log_fn$&"A" 

25598 ASSIGN @Diskb TO Log_fn$&"B" 

25599 ! ! ! ! ASSIGN @Diskc TO Log_fn$&"C" 

25600 ! 

25601 ! Find out where we are - retrieve number of observations recorded 

25602 ! 

25603 ENTER @Disk, 1 ;Data_set_title$ ,Max_obs_rec ,Max_vars , Logged_var_n 
ame$(*) .Numsubf ile , Subfile_names$(*) , Subfile_chartst(*) 

25604 Num_obs_recd-Subf ile_chartst(l) 

25605 This_obs-MAX(Num_obs_recd+l , 1) ! bump the observations pointer 

25606 IF This_obs<-Max_obs_rec THEN ! move the data to the observatio 

n record 

25607 FOR 1-1 TO Rcvd_vars 

25608 Obs_rec(I ,This_obs)-Pkt_avg(I) 

25609 NEXT I 

25610 Num_obs_recd-This_obs ! prepare to save the observation poin 
ter 

25611 Ending_obs-Num_obs_recd 

25612 Subfile_chartst(l)-Num_obs_recd 

25613 ! 

25614 FOR 1-1 TO Nummem 

25615! SET UP DATA FILE MAP FOR TRACES 1 THRU NUMMEM 

25616 Mivar(14)-I 

25617 Mvar ( 2 ) -1 

25618 CALL Routine ( 22 , 239) ! TURN ON data files 

25619! 

25620 Temp$ [1]-"R" 

25621 SELECT I 

25622 CASE 1 TO 4 

25623 Temp$ [ 2 ]-"A" 

25624 CASE 5 TO 7 

25625 Temp$[2]-"B" 

2‘5626 CASE 8 TO 99 

25627 Temp$[2]-"C" 

25628 END SELECT 

25629 SELECT INT(0bs_rec(2,This_obs) ) 

25630 CASE 0 TO 9 

25631 Temp$[3 ,4]-"0"&VAL$(INT(Obs_rec(2 ,This_obs) ) ) !RUN 

25632 CASE 10 TO 99 

25633 Temp$ [3,4] -VAL$ ( INT(Obs_rec (2 , This_obs ) ) ) ! RUN 

25634 END SELECT 

25635 SELECT INT(0bs_rec(3 ,This_obs) ) 

25636 CASE 0 TO 9 

25637 Temp$ [ 5 , 6 ]-"0"&VAL$(INT(Obs_rec(3 , This_obs) ) ) ! POINT 
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25638 

25639 

25640 

25641 

25642 

25643 

25644 

25645 

25646 

25647 

25648 

25649 

25650 ! ! ! 

25651 ! ! ! ! 

$(I) 

25652 ! ! ! ! 

25653 

25654 

25655 ! 

25656 

25657 

25658 

25659 

25660 

25661 

25662 

25663 ! 

25664 

25665 

25666 

25667 

25668 ! 

25669 ! 

25670 ! 

25671 Enable link: ! SEND A PACKET TO THE OTHER CPU WHICH CONTAINS THE "GO" WORD 


25672 ! 

25673 ! 

25674 

Go word$-"GO" 


25675 

OUTPUT Pkc sc; Go word$ 


25676 

RETURN 


25677 ! 

25678 ! 


- 

25679 ! 




25680 Disable_link: ! SEND A PACKET TO THE OTHER CPU WHICH CONTAINS THE "STOP" WO 
RD 

25$81 ! 

25682 Go_vord$-"STOP" 

25683 OUTPUT Pkt_sc ; Go_word$ 

25684 RETURN 

25685 ! 

25686 ! 

25687 ! 

25688 ! 

25689 ! 

25690 ! 

25691 Get_packet: ! 

25692 ENTER Pkt_sc;Pkt$(*) 

25693 ! 

25694 ! ! ! PRINTER IS PRT 

25695 CALL Pline(0 , "RECEIVED PACKET FOR SAMPLE "&VAL$(I)&" ") 


CASE 10 TO 99 

Temp$ [ 5 , 6 ] -VAL$ ( INT(Obs_rec ( 3 , This_obs ) ) ) ! POINT 
END SELECT 
SELECT I 
CASE 1 TO 9 

Temp$[7,8]-"0"&VAL$(I) 

CASE 10 TO 99 

Temp$ [7,8] -VAL$ ( I ) 

END SELECT 
F i lename$ ( I ) -Temp$ 

Mivar(14)-I 

CALL Routine(21 , 232) ! DISPLAY THE FILENAME 
PRINTER IS PRT 

PRINT "FOR OBS # ".Recnum;", THE FILENAME IS: Filename 

PRINTER IS CRT 
NEXT I 

GOSUB Compute_vars 

GOSUB Log_vars 
Ending_obs-Nuo_obs_recd 

ELSE 

CALL Pline(0, " LOG FILE FULL ") 

WAIT 1 

CALL Pline(0, " ") 

END IF 

CALL Pline(0, "OBSERVATION LOGGING COMPLETE ") 

WAIT 1 

CALL Pline(0, " ") 

GOTO End 4108 
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25696 

25697 

25698 

25699 

25700 

25701 

25702 

25703 

25704 

25705 

25706 

25707 

25708 

25709 

25710 

25711 

25712 

25713 

25714 

25715 

25716 

25717 

25718 

25719 

25720 


! ! ! FOR Vars_ctr-1 TO RCVD_vars 
M! PRINT Pkt$(Vars_ccr+l) 

! ! NEXT Vars_ctr 
! ! ! PRINTER IS CRT 
RETURN 

I 

I 

i 

| 

Compute_vars : ! compute values associated with mean logged_vars 
GOSUB Tun_vars 

FOR Probe-Initial_probe TO Ending_probe 
SELECT Probe 
CASE 1 

GOSUB Wall_vars 
CASE 2 

GOSUB Floor_vars 
CASE 3 

! GOSUB OTHER_VARS 
END SELECT 
NEXT Probe 
RETURN 

! 

Tun vars: ! 


25721 ! Construct variables for the tunnel conditions 

25722 Strut$-"TUNNEL" 

25723 GOSUB Compute_u_rho_m! Compute velocity, density, massflow 

25724 Obs_rec(46 ,This_obs)-Local_velocity 

25725 Obs_rec(47 ,This_obs)-Local_density 

25726 RETURN 

25727 ! 

25728 Wall vars: ! 


25729 ! Construct variables for the wall strut computed variables logfile (A) 

25730 A_rec(l , This_obs)-Obs_rec(l ,This_obs) ! TEST 

25731 A_rec(2 , This_obs)-0bs_rec(2 ,This_obs) !RUN 

25732 A_rec(3 ,This_obs)-Obs_rec(3 ,This_obs) ! POINT 

25733 A_rec (4 , This_obs )-0bs_rec ( 11 , This_obs ) ! P_T0TAL 

25734 A_rec(5 ,This_obs)-Obs_rec(12 ,This_obs) ! P_STATIC 

25735 A_rec(6 ,This_obs)-Obs_rec(8 ,This_obs)+460 !T_T0TAL (local) !Us 
e TUNNEL total 


25736 ! 

25737 

25738 

25739 

25740 

25741 

25742 

25743 

25744 

25745 

25746 

25747 

25748 

25749 

25750 

25751 

25752 

25753 

25754 


Strut$-"WALL" 

GOSUB Compute_u_rho_m! Compute velocity, density, mass flow 

A_rec ( 7 , This_obs ) -Local_ve loc i ty 

A_rec(8 ,This_obs)-Local_density 

A_rec ( 9 , This_obs )-Ts 

IF A_rec(7 , This_obs)>0 THEN 

A_rec(10,This_obs)-LGT(A_rec(7 ,This_obs) ) !Log(U) 
END IF 

IF A_rec(8 ,This_obs)>0 THEN 

A_rec ( 11 , This_obs ) -LGT ( A_rec ( 8 , This_obs ) ) !Log(Rho) 

END IF 

IF A_rec(6 ,This_obs)>0 THEN 

A_rec(12 ,This_obs)-LGT(A_rec(6,This_obs) ) !Log(T0) 
END IF 


IF Local_mass_flow>0 THEN 

A_rec(13 , This_obs)-LGT(Local_mass_flow) ! Log(Rho_inf*U) 
END IF 


A_rec(14, This_obs)-A_rec(ll ,This_obs)*A_rec(12 , This_obs) ! Log(Rh 
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A_rec( 15 ,This_obs)-A_rec(10, This_obs)*A_rec( 11 ,This_obs) ! Log(U) 
A_rec(16 ,This_obs)-A_rec(10 , This_obs)*A_rec(12 , This_obs) ! Log(U) 
A_rec(17 ,This_obs)-A_rec(15,This_obs)*A_rec(12 ,This_obs) ! Log(U) 


END IF 

SELECT Obs_rec(36 ,This_obs) 

CASE I TO 16 

A_rec (41+Wire , This_obs ) -Gain_code_14 ( Obs_rec ( 36 , This_ob 
CASE ELSE 

A_rec (4i+Wire , This_obs ) -0 ! Default gain-0 
END SELECT 
NEXT Wire 
RETURN 


o ) *Log(T0) 

25755 
*Log(Rho) 

25756 
*Log(T0) 

25757 

*Log(Rho)*Log(T0) 

25758 FOR Wire-1 TO 4 

25759 A_rec(17+Wire ,This_obs)-0bs_rec(19+Wire , This_obs) 

25760 IF A_rec(17+Wire ,This_obs)>0 THEN 

25761 A_rec(21+Wire , This_obs)-LGT(A_rec(17+Wire , This_obs ) ) ! Lo 
g(E) 

25762 

25763 

25764 

25765 
s)+l) 

25766 

25767 

25768 

25769 

25770 

25771 ! 

25772 Floor_vars: ! 

25773 ! Construct variables for the floor strut computed variables logfile (B) 

25774 B_rec ( 1 , This_obs ) -Obs_rec ( 1 , This_obs ) ! TEST 

25775 B_rec(2,This_obs)-Obs_rec(2,This_obs) !RUN 

25776 B_rec ( 3 , This_obs )-Obs_rec ( 3 , This_obs ) .'POINT 

25777 B_rec ( 4 , This_obs )-0bs_rec ( 14 , This_obs ) ! P_T0TAL (local) 

25778 B_rec ( 5 , This_obs)-Obs_rec (15 , This_obs ) ! P_STATIC ( local) 

25779 B_rec(6 , This_obs)-0bs_rec(8 ,This_obs)+460 !T_T0TAL (local) 

Use TUNNEL total 

25780 ! 

Strut$-" FLOOR" 

GOSUB Compute_u_rho_m! Compute velocity, density .mass flow 
B_rec ( 7 , This_obs ) -Local_ve loc i ty 
B_rec ( 8 , This_obs ) -Local_dens i ty 
B_rec(9 ,This_obs)-Ts 
IF B_rec ( 7 , This_obs ) >0 THEN 

B_rec(10 ,This_obs)-LGT(B_rec(7 , This_obs) ) 

END IF 

IF B_rec(8 ,This_obs)>0 THEN 

B_rec(ll,This_obs)-LGT(B_rec(8 ,This_obs) ) 

END IF 

IF B_rec ( 6 , This_obs )>0 THEN 

B_rec(12 ,This_obs)-LGT(B_rec(6 ,This_obs) ) 

END IF 

IF Local_mass_flow>0 THEN 

B_rec(13 ,This_obs)-LGT(Local_mass_flow) 

END IF 

B_rec ( 14 , This_obs) -B_rec ( 11 , This_obs ) *B_rec ( 12 , This_obs ) ! Log (Rh 


!Log(U) 

!Log(Rho) 


25781 

25782 

25783 

25784 

25785 

25786 

25787 

25788 

25789 

25790 

25791 

25792 

25793 

25794 

25795 

25796 

25797 

25798 

o)*Log(T0) 

25799 
*Log(Rho) 

25800 
*Log(T0) 

25801 
*Log(Rho)*Log(T0) 

25802 FOR Wire-1 TO 3 

25803 B_rec ( 17+Wire , This_obs )-Obs_rec ( 23+Wire , This_obs ) 


! Log(TO) 
!Log(Rho_inf*U) 


B_rec ( 15 , This_obs ) -B_rec ( 10 , This_obs ) *B_rec ( 11 , This_obs ) ! Log (U) 
B_rec ( 16 , This_obs )-B_rec ( 10 , This_obs ) *B_r ec (12 , This_obs ) ! Log (U) 
B_rec ( 17 , This_obs ) -B_rec ( 15 , This_obs ) *B_rec ( 12 , This_obs ) ! Log(U) 
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25804 

25805 

25806 

25807 

25808 

25809 

is_obs)+l) 


IF B_rec ( 17+Wire , This_obs)>0 THEN 

B_rec ( 21+Wire , This_obs )-LGT(B_rec( 17+Wire , This_obs ) ) 

END IF 

SELECT Obs_rec( 36+Wire , This_obs) 

CASE L TO 7 

B_rec(41+Wire , This_obs)-Gain_code_57 (Obs_rec(36+Wire ,Th 


25810 

25811 

25812 

25813 

25814 

25815 ! 

25816 Other vars: 


CASE ELSE 

B_rec(41+Wire ,This_obs)-0 '.Default gain-0 
END SELECT 
NEXT Wire 
RETURN 

t 


25817 ! 

25818 

25819 

25820 

25821 

25822 

25823 

25824 ! 

25825 

25826 

25827 

25828 

25829 

25830 

25831 

25832 

25833 

25834 

25835 

25836 

25837 

25838 

25839 

25840 

25841 

25842 

o)*Log(T0) 

25843 
*Log(Rho) 

25844 
*Log(T0) 

25845 

*Log(Rho)*Log(T0) 

25846 

25847 

25848 

25849 

25850 

25851 

25852 

25853 

25854 

25855 

25856 

25857 

25858 


Construct variables for the 'other' strut computed variables logfile (C) 
C_rec ( 1 , This_obs ) -Obs_rec ( 1 , This_obs ) ! TEST 
C_rec(2 ,This_obs)-0bs_rec(2 ,This_obs) !RUN 
C_rec(3 ,This_obs)-Obs_rec(3 ,This_obs) ! POINT 
C_rec(4,This_obs)-Obs_rec(17 ,This_obs) !P_T0TAL (local) 

C_rec(5 ,This_obs)-0bs_rec(18 ,This_obs) ! P_STATIC( local) 

C_rec(6 ,This_obs)-Obs_rec(19 ,This_obs) !T_TOTAL (local) 

Strut$-”0THER" 

GOSUB Compute_u_rho_m ! Compute velocity .density .mass flow 

C_rec ( 7 , This_obs ) -Local_ve loc i ty 

C_rec(8 , This_obs)— Local_density 

C_rec(9 ,This_obs)-Ts 

IF C_rec(7 ,This_obs)>0 THEN 

C_rec(10,This_obs)-LGT(C_rec(7,This_obs) ) ! Log(U) 

END IF 

IF C_rec ( 8 , This_obs )>0 THEN 

C_rec(ll ,This_obs)-LGT(C_rec(8 , This_obs) ) ! Log(Rho) 

END IF 

IF C_rec(6 ,This_obs)>0 THEN 

C_r e c ( 1 2 , Thi s_ob s ) -LGT ( C_r e c ( 6 , Th i s_ob s ) +4 6 0 ) !Log(T0) 

END IF 

IF Local_mass_flow>0 THEN 

C_rec(13 ,This_obs)-LGT(Local_mass_flow) !Log(Rho_inf*U) 
END IF 

C_rec ( 14 , This_obs )-C_rec ( 11 , This_obs ) *C_rec ( 12 , This_obs ) ! Log (Rh 
C_rec(15 ,This_obs)-C_rec(10 ,This_obs)*C_rec(ll ,This_obs) !Log(U) 
C_rec ( 16 , This_obs )-C_rec ( 10 , This_obs ) *C_rec ( 12 , This_obs ) ! Log (U) 
C_rec ( 17 , This_obs )-C_rec ( 15 , This_obs ) *C_rec ( 12 , This _obs ) ! Log(U) 
FOR Wire-1 TO 3 

C_rec (17+Wire , This_obs ) -Obs_rec ( 2 3+Wire , This_obs ) 

IF C_rec( 17+Wire, This_obs)>0 THEN 

C_rec ( 21+Wire , This_obs ) -LGT ( C_rec ( 17+Wire , This_obs ) ) 

END IF 
NEXT Wire 
RETURN 


Compute_u_rho_m : 
SELECT 


! COMPUTE 
Strut$ 


LOCAL VELOCITY, DENSITY, MASS FLOW 
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CASE "TUNNEL" 

Pt-0bs_rec(7 , This_obs) 
Ps-0bs_rec(6 ,This_obs) 
Tt-0bs_rec(8 ,This_obs)+460 
CASE "WALL" 

Pt-A_rec (4 , This_obs ) 
Ps-A_rec ( 5 , This_obs ) 
Tt-A_rec(6 ,This_obs) 

CASE "FLOOR" 

Pt-B_rec ( 4 , This_obs ) 
Ps-B_rec ( 5 , This_obs ) 
Tt-B_rec ( 6 , This_obs ) 

CASE "OTHER" 

Pt-C_rec (4 , This_obs) 
Ps-C_rec(5 ,This_obs) 
Tt-C_rec ( 6 , This_obs ) 

END SELECT 


IF 


25859 

25860 

25861 

25862 

25863 

25864 

25865 

25866 

25867 

25868 

25869 

25870 

25871 

25872 

25873 

25874 

25875 

25876 ! 

25877 

25878 

25879 

25880 

25881 

25882 

25883 

25884 

25885 

25886 

25887 

25888 

25889 

25890 

25891 

25892 

25893 

25894 

25895 

25896 

25897 

25898 

25899 

25900 

25901 ! 

25902 ! 

25903 ! 

25904 

25905 ! 

25906 ! 

25907 

25908 ! 

25909 

name$(*) .Numsubfile , Subfile_names$(*) , Subf ile_chartst(*) 


Ps-0 THEN 
Prac-0 
GOTO 25883 
END IF 
P_rat-Ps/Pt 

IF P_rat>l THEN P_rat-1 
IF P_rat<-0 THEN 
Ts-MAX(.l.Tt) 

Local_ve loci ty-0 . 

Local_dens i ty-Ps/ (53. 3*Ts ) 

Local_mass_f low*-0 . 

ELSE 

Ts-(Tt)/(P_rat A ( - . 285714285) ) 

IF P_rat-1 THEN 

Local_velocity-0 . 

ELSE 

Sound-SQR(2402. 764*Ts) 

Local_mach-SQR(5*(P_rat A ( - . 285714285) - 1) ) 
Local_velocity-Local_mach*Sound! U 
END IF 

Local_density—Ps/(53 . 3*Ts) ! Rho 

Local_mass_f low-Local_ve loci ty*Local_dens icy! M 
END IF 
RETURN 


Log_vars: ! LOG THE CURRENT OBSERVED VARIABLES 


ASSIGN @Disk TO Log_fn$ 

OUTPUT @Disk, 1 ; Data_set_title$ ,Max_obs_rec , Max_vars , Logged_var_ 


25910 

25911 

25912 

25913 

25914 

25915 

25916 

) , Numsubfile, A_subfile_names$(*) , A_sub_chartst(*) 


ENTER @Disk, 2 
OUTPUT @Disk; Obs_rec (*) 

ASSIGN @Disk TO * 

ASSIGN (ZjDiska TO Log_fn$&"A" 

OUTPUT @Diska, 1 ;A_set_title$ ,Max_obs_rec ,Max_vars , A_var_name$(* 
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ASSIGN 0DISKC TO LOG FN$&"C" 


RETURN 


"SAMPLES TO AVG" 


25917 ENTER @Diska,2 

25918 OUTPUT @Dlska;A_rec(*) 

25919 ASSIGN @Dlska TO * 

25920 ! 

25921 ASSIGN @Dlskb TO Log_fn$&"B" 

25922 ! 

25923 OUTPUT @Dlskb , 1 ; B_set_tltle$ .Maxobsrec ,Max 
) , Numsubf lie , B_subf lle_names$ (*) , B_sub_chartst(*) 

25924 ENTER @Dlskb,2 

25925 OUTPUT @Dlskb;B_rec(*) 

25926 ASSIGN @Dlskb TO * 

25927 ! 

25928 ! ! ! ! ASSIGN @DISKC TO LOG_FN$&"C" 

25929 I 

25930 !!!! OUTPUT @Dlskc , 1 ; C_set_title$ , Max_obs 

_name$(*) , Numsubf lie , C_subfile_naraes$(*) , C_sub_chartst(*) 

25931 ! ! ! ! ENTER @Dlskc,2 

25932 !!!! OUTPUT @Dlskc;C_rec(*) 

25933 !!! ! ASSIGN @Dlskc TO * 

25934 ! 

25935 ! 

25936 RETURN 

25937 ! 

25938 ! 

25939 ! 

25940 End_4108 : ! 

25941 END SELECT 

25942 CASE 4110 ! "SAMPLES TO AVG" 

25943 SELECT Routinl 

25944 CASE 21 

25945 Mvar(2)-Num_avgs 

25946 Mvar(4)-1 

25947 Mvar(5)-1. 

25948 Mvar(6)-300. 

25949 Mlvar(9)-0 

25950 CASE 22 

25951 Num_avgs-Mvar ( 2 ) 

25952 END SELECT 

25953 CASE 4112 ! "INITIAL OBS" 

25954 SELECT Routinl 

25955 CASE 21 

25956 Mvar(2)-Initial_obs 

25957 Mvar(4)-1 

25958 Mvar(5)-1 . 

25959 Mvar(6)-300. 

25960 Mivar(9)-0 

25961 CASE 22 

25962 Inltlal_obs-Mvar(2) 

25963 END SELECT 

25964 CASE 4113 ! "ENDING OBS" 

25965 SELECT Routinl 

25966 CASE 21 

25967 Mvar(2)-Endlng_obs 

25968 Mvar(4)-1 

25969 Mvar(5)-1 . 

25970 Mvar (6)-300 . 

25971 Mivar(9)-0 

25972 CASE 22 

25973 Ending_obs-Mvar(2) 

'25974 END SELECT 


"INITIAL OBS" 


vars.B var name$(* 


rec.Max vars.C var 
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25975 ! 

25976 ! 

25977 ! 

25978 

25979 

25980 

25981 

25982 

25983 

25984 

25985 

25986 

25987 

25988 

25989 

25990 

25991 

25992 

25993 

25994 

25995 

25996 

25997 

25998 

25999 

26000 ! 
26001 ! 
26002 ! 

26003 

26004 

26005 

26006 

26007 

26008 

26009 

26010 
26011 
rs) 
26012 

26013 

26014 

26015 

26016 

26017 

26018 

26019 

26020 
26021 
26022 
rs) 

26023 

26024 

26025 

26026 

26027 

26028 

26029 

26030 

26031 

26032 


CASE 4130 ! "INITIAL PROBE" 

SELECT Routlnl 
CASE 21 

Mvar(2)-Inltial_probe 

Mvar (4)-l 

Mvar(5)-1. 

Mvar(6)-2 . 

Mivar (9)-0 
CASE 22 

Initial_probe-Mvar(2) 

END SELECT 

CASE 4131 ! "ENDING PROBE" 

SELECT Routlnl 
CASE 21 

Mvar(2)-Ending_probe 

Mvar(4)-1 

Mvar(5)-1 . 

Mvar(6)-2 . 

Mivar(9)-0 
CASE 22 

Endlng_probe-Mvar ( 2 ) 

END SELECT 


CASE 4127 ! "LOGFILE TO PC" 

SELECT Routlnl 
CASE 31 

Pc_device-8 ! Send data out thru HP- IB bus 8 

! bus 8 is NOT system controller 
FOR Probe-Initial_probe TO Ending_probe 
SELECT Probe 
CASE 1 

OUTPUT Pc_device;VAL$(Ending_obs-Initial_obs+l) ,VAL$(Max_va 
FOR 1-1 TO Max_vars 

DISP "NAME( "&VAL$(I)&" ) IS: "&A_var_name$(I) 

OUTPUT Pcdevice ; A_var_name$( I ) [ 1 , 10 ] 

NEXT I 

FOR Thls_obs-Initial_obs TO Ending_obs 
FOR 1-1 TO Max_vars 

OUTPUT Pc_device ; A_rec ( I , This_obs ) 

NEXT I 

NEXT Thlsobs 
CASE 2 

OUTPUT Pc_devlce ; VAL$(Endtng_obs- Initial_obs+l) , VAL$(Max_va 
FOR 1-1 TO Max_vars 

DISP ’* NAME ( ” &VAL$ ( I ) &" ) IS: ”&B_var_name$(I) 

OUTPUT Pc device ; B_var_name$ ( I ) [ 1 , 10 J 
NEXT I 

FOR This_obs-Initial_obs TO Ending_obs 
FOR 1-1 TO Max_vars 

OUTPUT Pc_device ; B_rec(I ,This_obs) 

NEXT I 

NEXT Thlsobs 
END SELECT 


26033 NEXT Probe 

26034 END SELECT 

26035 END SELECT 

26036 END SELECT 

26037 SUBEND 

26038 ! 

26039 ! 

26040 ! 

26041 ! 
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APPENDIX B. Function definitions 

This appendix contains the DDAS functions added to ACQUIRE by 
the user to provide the necessary functionality to acquire, 
process, display, store and transmit the hot wire data. 

The function definition sheets provide essential definition 
data for each function. Full documentation for each function 
is contained within the program source code listings. 
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CALC VEL etc 


Function Name : 
Function- Number : 

Module : 


CALC VEL etc 
4104 

M0DUSR2* 


This function computes velocity, density and temperature 
fluctuations as ratios of the fluctuating quantitys to the 
mean quantitys: 

u'/IJ p'/P t'/T 

The log file for each observation in the range of INITIAL OBS 
to ENDING OBS is processed for each probe in the range of 
INITIAL PROBE to ENDING PROBE. 

For each observation, and each probe, the three fluctuating 
data files related to the Run and Point (variables 2 and 3 in 
the observation record) are loaded into traces (and memories) 
1,2, and 3 . 

The computations retrieve mean values, sensitivities, and 
gains related to the data in channels 1, 2 and 3. 

For each simultaneous sample in each of the three traces (the 
beginning and ending samples are defined by the cursors on 
trace 1) the instantaneous value is divided by the equivalent 
mean value and the gain, then the three ratios are matrix 
multiplied by the nmatrix inversion of the sensitivities. This 
process essentially solves a set of simultaneous equations for 
three unknowns: u'/U, p'/P, and t'/T. 

The solutions are placed in traces (memorys) 4, 5 and 6. The 

mean value is removed, and the RMS (root mean square) value of 
each of the answers is stored in the logfile for the 
appropriate observation and probe. Traces 4, 5, and 6 are then 

stored in separate disk files. 


[ CALC VEL etc ] 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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CAT GROUP 


Func Cion- Name : 

Function Number: 
Module : 


CAT GROUP 
4021 

MODUSR1* 


This function catalogs all selected files. The selection criteria is 
defined in detail in the HP BASIC language refernce manual for HP 
function CAT (SELECT). 


[ CAT GROUP ] 

I I 

I I 

| - - [ - ] - - [ name ] - - | 


Item 

| DescriD t ion 

- 1 

Ranee 

name 

| string expression 
1 
1 
1 

1 

1 

1 

1 

any valid 
char ac ter s 
that are allowed 
in a file name 


*. This routine - MODUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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CODE TO GAINS 


Func t i on Nam? : 
Function Number: 

Module : 


CODE TO GAINS 
4125 

MODUSR2* 


This function computes actual voltage gains fom gain codes. 

The fluctuating data is gained to provide adaquate voltage for 
filtering an d digitization, and the gain codes are sent in a 
data packet from the static data computer during the logging 
of a data point. For the probes in the range INITIAL PLROBE to 
ENDING PROBE, and for observations in the range INITIAL OBS to 
ENDING OBS, the gain codes are retrieved from the observation 
file, and thru a table lookup algorithm, the actual gains are 
retrieved, and stored back into the appropriate place in the 
logfile for ht eprobe and observation being processed. 


[ CODE TO GAINS ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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COEF FILENAME 


Function Name: COEF FILENAME 

Function- Number: 4101 

Module: M0DUSR2* 


This function selects the file name relaven to the 
hotwire coeficient file being processed. See functions LOAD 
COEFS and STORE COEFS . 


[ COEF FILENAME ] 


*. This routine - MODUSR2 - was written by Steven J. C.ukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB ider 
contract NAS1-17919, Task 36. This work was done begin; .g in 
October, 1986 and continues thru April, 1988. 
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COMPUTE COEFS 


Function Name: COMPUTE COEFS 

Function- Number: 4123 

Module: M0DUSR2 


This function computes hotwire coefficients for the probe 
previously selected - by function INITIAL PROBE - utilizing a 
custom multiple linear regression rouutine that generates up 
to 10 coefficients based on calibration data already in the 
logfile. These coefficients are then stored - by function 
STORE COEFS - in a coefficient file whose name has been 
previously defined - by function COEF FILENAME. Function 
COMPUTE SENS utilizes these coefficients to generate 
senstivities necessary to compute velocity, density and 
temperature turbulence ratios -by function CALC VAL etc. 


[ COMP COEFS ] 


This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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COMPUTE R etc 


Function Name : 
Function- Number : 
Module : 


COMPUTE R etc 

4132 

M0DUSR2 


This function computes the correlations between velocity, 
density and temperature fluctuations and then computes mass 
flow fluctuation [ m ' /M ( rms ) ] and pressure fluctuation 
[p'/P( rras )] using the correlation between velocity, density, 
and temperature fluctuations. 

The log file for each observation in the range of INITIAL OBS 
to ENDING OBS is processed for each probe in the range of 
INITIAL PROBE to ENDING PROBE. 

For each observation, and each probe, the three fluctuating 
data files related to the Run and Point (variables 2 and 3 in 
the observation record) are loaded into traces (and memories) 
2 , 3 , and 4 . 

As each pair of traces is multiplied together, the resulting 
trace ends up in trace 1. The correlation of the two traces 
multiplied is the ratio of the rms value of the first trace 
multiplied by the rms value of the second trace to the mean 
value of trace 1. The rms values of the first and second 
traces have been previously calculated by function CALC VEL 
etc, and were called u'/ u > p'/P» and t ' /T . 

The correlation between these three components of turbulence 
are stored in the appropriate observation record for the 
observation and probe being processed as R(RhoU), R(UTO), and 
R(RhoTO) . 

Massflow and pressure fluctuations are then computed from the 
correlations just computed, and these are also stored in the 
appropriate logfile as M ' /M and P' /P. 


[ COMPUTE R etc ] | 


This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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COMPUTE SENS 


Function Name: COMPUTE SENS 


Function Number: 
Module : 


4111 

MODUSR2 


This function computes hotwire sensitivities for the 
probes in the range INITIAL PROBE to ENDING PROBE for all 
observations in the range INITIAL OBS to ENDING 
OBS . The hotwire coefficients previously defined - see 
functions GET COEFS, COEF FILENAME, LOAD COEFS , ENTER COEFS , 
EDIT COEFS, STORE COEFS, AND COMPUTE COEFS. 

The computed sensitivities are stored in the appropriate probe 
file for each appropriate observation. 

These senstivities are read from the logfile - by function 
CALC VAL etc - to compute velocity, density and temperature 
turbulence ratios. 


[ COMPUTE SENS ] 


This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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EDIT COEFS 


Function Name : 
Function Number: 

Module : 


EDIT COEFS 
4133 

M0DUSR2* 


This function allows the operator to manually edit hot wire 
calibration coefficients (up to 10) for three wires. These 
coefficients have been previously generated. The operator 
views a list of the entered coefficients, and, by 
responding to prompts, select the coefficient to be 
edited. After each coefficient is edited, the operator may 
accept the coefficients, or be prompted to select another 
coefficient to be edited. The operator should then invoke 
function STORE COEFS. 


[ EDIT COEFS ] 


NOTE: This function is highly interactive, and is no t 

recommended for inclusion in a SEQUENCE PROGRAM. 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 


C-3- 
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ENDING OBS 


Function Name : 
Function Number : 

Module : 


ENDING OBS 
4113 

M0DUSR2* 


This function declares the ending observation to be processed 
by other M0DUSR2 functions - see function INITIAL OBS. 
Functions utilizing this feature to define the range of 
observations to be processed are: COMPUTE SENS, LOG DATA 

POINT, CODE TO GAINS, CALC VEL etc, COMPUTE R etc, Remake 
probe and LOGFILE TO PC. 


[ ENDING OBS ] 


I - - [ - ] - - [value ] - - | 


Item I Description I Range 

value | numeric integer | 1 to 300, the 

| | max number of 

j j observations 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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ENTER GOEFS 


Function Name : 
Function Number: 

Module : 


ENTER COEFS 
4103 

MODUSR2* 


This function allows the operator to manually enter hot wire 
calibration coefficients (up to 10) for three wires. These 
coefficients have been previously generated elsewhere, and are 
not available for entry via disc (see function LOAD COEFS). 

The operator is prompted for each coefficient by wire, number, 
and name. Once all 30 coefficients are entered (unused 
coefficients should be set to 0.0), the operator views a list 
of the entered coefficients, and chooses to accept or reject 
the coefficients. If they aree accepted, the function is 
complete. The operator should then invoke function STORE 
COEFS. If the coefficients are not accepted - because they 
are not correct, this function automatically enters the EDIT 
COEFS function. 


[ ENTER COEFS 


NOTE: This function is highly interactive, and is not 

recommended for inclusion in a SEQUENCE PROGRAM. 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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FILE COPY 


Func t ion Name : 
Function Number: 

Module : 


FILE COPY 
4025 

MODUSR1* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ FILE COPY ] 


*. This routine * MODUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 


Page 100 


FILE GROUP 


Func t ion Name : 
Function Number: 

Module : 


FILE GROUP 
4028 

M0DUSR1* 


This function defines the files to be selected for copying - 
function COPY FILES - or moving - function MOVE FILES. the files 
selected will begin with, or be equal to the character(s) defined 
by this function. For example, if this function defines "ABC' 1 , then 
all of the files in the TO DISK device that begin with M ABC ,r would 
be selected for copying or moving. Note that this function only 
defines the char ac t e r ( s ) : no selection is done until COPY FILES or 
MOVE FILES is invoked. 


[ FILE GROUP ] 

I I 

I I 

| - - [ = ] - - [ name ] - - | 


Item 


1 Description 


Range 


name 


| string expression 


any valid 
characters 
that are allowed 
in a file name 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 



Page 101 


FILE TRANSFERS 


Function Name: 
Function Nunrber : 

Module : 


FILE TRANSFERS 
4126 

M0DUSR2* 


This function serves merely as a label for the path for ac- 
cessing other related functions via the softkeys. 


[ FILE TRANSFERS ] 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NASl-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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FILE UTIL I T Y S 


Func t ion Name : 
Func t i on Numb e r : 

Mo dul e : 


FILE UTILITYS 
4020 

M0DUSR1 * 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ FILE UTILITYS 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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FROM DISK 


I 

I 

i 

i 

Function Name: FROM DISK 

Function ' Number : 4026 

Module: MODUSR1* 


This function defines the mass storage device from which the 
files will be copied - function COPY FILES - or moved - 
function MOVE FILES. 


[ FROM DISK 






] -- [ MSI].-- | 


Item 

| Descr I d t ion 

| Ranee 

value 

| MSI 

| any valid HP 


1 

| storage 


1 

1 device 


*. This routine - MODUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April,1988. 
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GET GO EF 


Function Name : 
Function Num-ber: 

Module : 


GET COEF 
4122 

MODUSR2 * 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ GET COEF ] 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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HOTWIRE MENU 


Func t ion Name : 
Function Number: 

Module : 


HOTWIRE MENU 
4100 

M0DUSR2* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ HOTWIRE MENU ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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HOTWIRE CALC 


Func t i on Name : 
Function- Number : 

Module : 


HOTWIRE CALC 
4100 

MODUSR2* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ HOTWIRE CALC ] 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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INITIAL OBS 


Func t ion Name : 
Function -Number : 

Module : 


INITIAL OBS 
4112 

M0DUSR2* 


This function declares the be ginning observation to be 
processed by other M0DUSR2 functions - see function ENDING 
OBS. Functions utilizing this feature to define the range of 
observations to be processed are: COMPUTE SENS, LOG DATA 
POINT, CODE TO GAINS, CALC VEL etc, COMPUTE R etc, Remake 
probe and LOGFILE TO PC. 


[ INITIAL OBS ] 

I I 

I I 

| - - [ - ] - - [value ] - - | 


Item 

1 DescriDtion 

1 Ranee 

value 

| numeric integer 

| 1 to 300, the 


1 

| max number of 


1 

I observations 




*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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INITIAL PROBE 


Func t i on Name : 
Function Number : 

Mo dul e : 


INITIAL PROBE 
4130 

M0DUSR2* 


This function declares the beginning probe to be processed by 
other M0DUSR2 functions - see function ENDING PROBE. Functions 
utilizing this feature to define the range of probes to be 
processed are: COMPUTE SENS, LOG DATA POINT, CODE TO GAINS, 
CALC VEL etc, COMPUTE R etc, Remake Probe and LOGFILE TO PC. 


[ INITIAL PROBE ] 


I 

I I 

I I 

| - - [ - ] - - [value ] - - | 


Item | Description I Range 

value | numeric integer | 1 to 3 , the 

j j max number of 

I | probes 


*. This routine - M0DUSR2 - was written by Stevn J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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LOAD COEFS 


Func Cion Name : 
Function- Number : 

Module : 


LOAD COEFS 
4102 

M0DUSR2* 


This function loads coefficients previously stored by 
function STORE COEFS. 

The coefficient filename must have been previously defined - 
see function COEF FILENAME. 


[ LOAD COEFS ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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LOAD LOGFILE 


Function Name: LOAD LOGFILE 

Function- Number : 4106 

Module: MODUSR2* 


This function loads the files 
observation data. The functi 
name, and function DISC DEVIC 
Identifier. See LOG FILENAME 
loaded into memory. 


used to log hotwire calibration 
on LOG FILENAME declares the file 
E declares the Mass Storage 
for a description of the files 


[ LOAD LOGFILE 


NOTE: 

This function mus t be performed before the 
function LOG DATA POINT can be invoked. 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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LOAD LOGFILE 


Func t ion Name : 
Function Number: 

Module : 


LOAD LOGFILE 
4106 

M0DUSR2* 


This function loads the files used to log hotwire calibration 
observation data. The function LOG FILENAME declares the file 
name, and function DISC DEVICE declares the Mass Storage 
Identifier. See LOG FILENAME for a description of the files 
loaded into memory. 


[ LOAD LOGFILE ] 


NOTE: 

This function must be performed before the 
function LOG DATA POINT can be invoked. 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April,1988. 
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LOG DATA 


Func t ion Name : 
Function Number : 

Module : 


LOG DATA 
4109 

MODUSR2* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ LOG DATA ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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LOG DATA POINT 


Function Name: 
Function - Number : 

Module : 


LOG DATA POINT 
4108 

M0DUSR2* 


This function logs the hotwire calibration observation data 
into the logfile(s). The function LOG FILENAME declares the 
file name, and function DISC DEVICE declares the Mass Storage 
Identifier. See LOG FILENAME for a description of the files 
which receive the data. 

The hotwire calibration data is received from a static data 
acquisition system - MODCOMP, other HP, etc - thru a GPIB 
interface in the form of an ASCII data packet. The static 
data system is regularly sending about 1 packet per second, 
and it contains the necessary data in engineering units, 
including test conditions and test identification. 

This function also calculates various data items for inclusion 
in some of the data files. 

This function also generates fluctuating data file names base 
on the RUN, POINT and channel of the data. These names relate 
to the data being logged, and allow storing of the fluctuating 
data in appropriately named files. See functions STORE or 
TF/STORE ALL. 

Example: "RA031701" 

R -Realtime digitization 

A -probe A 

03-RUN 

17-POINT 

01-data channel 1 


[ LOG DATA POINT ] 


*. This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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LOGFILE TO PC 


Function Name : 
Function Number: 

Module : 


LOGFILE TO PC 
4127 

MODUS R2 * 


This function transfers logfile data to another computer via a 
GPIB. The information transmitted is all in ASCII to assure 
compatabili ty between systems. 

For each probe in the range INITIAL PROBE to ENDING PROBE and 
for observations in the range INITIAL OBS to ENDING OBS, the 
appropriate observations records are transmitted. Prior to 
transmitting the set of each probes observations, the names of 
the variables contained in the observation record are 
transmitted . 


[ LOGFILE TO PC ] 


t 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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MOVE FILES 


* 


Func t i on Name : 
Function Numb e r : 

Module : 


MOVE FILES 
4030 

MODUSRl* 


This function defines the files to be selected for moving, and then 
copy the selected files from a disk device to a disk device using 
the HP COPY command, when the files are selected, a report is 
generated on the printer which declares the number of files 
selected, the from device, the to device, and the files selected for 
copying. As each file is actually copied, a message is displayed on 
the CRT, and printed on the printer. Once all selected files have 
been copied, all successfully copied files are purged from the from 
device, completing the "move". Functions FROM DISK, TO DISK, AND 
FILE GROUP all effect the results of this function. Selection of 
files by indicating a group name with this function overrides the 
group name selection previously made by function FILE GROUP. 


[ MOVE FILES ] 


| - - [ - ] - - [ name ] - - | 



Item 

| Description 

| Ranee 


name 

| string expression 

1 

1 

1 

| any valid 
j characters 
j that are allowed 
| in a file name 

4 


1 

1 

1 

1 

| see FILE GROUP for 
j details 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 



Page 116 


PICTURE LOG E 


Function Name : 
Function Number : 

Module : 


PICTURE LOG E 
4006 

M0DUSR1* 


This function generates x-y plots on the CRT which represent 
the Log(RhoU) vs. Log(E) for each hotwire. The axes, titles, 
etc are internally generated using a plot file called "RhoU" - 
see function PLOT NAME. (It should be noted that this 'RhoU' 
picture is actually generated for dynamic channels 8 and 9 - 
which are assumed to be 1 point long, - and offscale as well.) 

This function represents data from all observations in the 
range INITIAL OBS to ENDING OBS for all probes in the range 
INITIAL PROBE TO ENDING PROBE. 


[ PICTURE LOG E ] 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PLOT EJECT 


Function Name: 
Function Number : 

Module : 


PLOT EJECT 
4002 

M0DUSR1* 


This function sends a "PG" command to the plot device (if the 
plot device is no t the CRT device. 

This function also resets the 'number of observations plotted' 
pointer, the number of observations printed' pointer, and the 
'ending observations' pointer, which affects the operation of 
functions ENDING OBS , PRNT LOGFILE, PLOT L0G_E and PICTURE 
LOG E . 


[ PLOT EJECT ] 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PLOT LOG E 


Function Name: PLOT LOG E 

Func t ion Numb e r : 4005 

Module: MODUSR1* 


This function generates x-y plots on the Plotter which 
represent the Log(RhoU) vs. Log(E) for each hotwire. The 
axes, titles, etc are internally generated using a plot file 
called "RhoU" - see function PLOT NAME. (It should be noted 
that this 'RhoU' picture is actually generated for dynamic 
channels 8 and 9 - which are assumed to be 1 point long, - and 
offscale as well.) 

This function represents data from all observations in the 
range INITIAL OBS to ENDING OBS for all probes in the range 
INITIAL PROBE TO ENDING PROBE. 


[ PLOT LOG E ] 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PLOT UTILITYS 


% 


Function Name: 
Function' Number: 

Module : 


PLOT UTILITYS 
4001 

M0DUSR1 * 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ PLOT UTILITYS ] 




This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PRNT LOGFILE 


Function Name : 
Function Nunfber : 

Module : 


PRNT LOGFILE 
4107 

MODUS R2* 


9 


This function prints the data hotwire calibration observation 
data currently in the logfile(s). The function LOG FILENAME 
declares the file name, and function DISC DEVICE declares the 
Mass Storage Identifier. See LOG FILENAME for a description 
of the files which contain the data printed. The format of 
the printout is customized to best demonstrate the hotwire 
calibration data. This function prints all observations 
beginning with INITIAL OBS, and ending with ENDING OBS unless 
previously printed. When the program first starts, and when 
the PLOT EJECT function is invoked, the number of observations 
printed is reset to zero, causing all observations already 
logged to be printed when PRNT LOGFILE is invoked. This 
technique allows the easy implementation of a sequence program 
loop including both function LOG DATA POINT and PRNT LOGFILE 
without having to manipulate observation pointers. 


[ PRNT LOGFILE 


ft 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PURGE 


Function Name: PURGE 

Func tion ' Number : 4022 

Module: MODUSRl* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ PURGE ] 


*. This routine - MODUSRl - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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PURGE FILE 


Function Name : 
Function ' Number : 

Module : 


PURGE FILE 
4024 

MODUSR1* 


9 


This function defines the file to be selected for pur ging ( de le t ing ) 
from the disk, and then purges the selected file. 


[ PURGE FILE ] 

I I 

I I 

| - - [ - ] - - [ name ] - - | 


Item 

1 DescriDtion 

I Ranee 

name 

| string expression 
1 
1 
1 
1 

| any valid 
j characters 
| that are allowed 
j in a file name 

i 


*. This routine - MODUSRl - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 


! 
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PURGE GROUP 


Func t ion 

Name : 

PURGE GROUP 

Func t ion 

Number : 

4023 

Module : 


M0DUSR1* 


This function purges all selected files. The selection criteria is 
defined in detail in the HP BASIC language refernce manual for HP 
function PURGE (SELECT). 


[ PURGE GROUP ] 

I I 

I I 

I * - [ “ ] - ' [ name ] - - | 


Item 


Description 


Range 


name 


string expression 


any valid 
characters 
that are allowed 
in a file name 


i 

! 

*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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Remake Probe 


Func t ion 

Name : 

Remake Probe 

Func t ion 

Numb e r : 

4128 

Module : 


M0DUSR2* 


This function reproduces the computations normally 
accomplished during the execution of the function LOG DATA 
POINT, but without actually acquiring any new data. The 
purpose is to recompute data should modifications to the 
computations become necessary. 

This function performs these calculations for all probes in 
the range INITIAL PROBE to ENDING PROBE, and for all 
observations in the range INITIAL OBS to ENDING OBS . 

NOTE: 

This function permanently overwrites previous data in the 
logf iles . 


[ LOG DATA POINT ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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SAMPLES TO AVG 


Function Name: SAMPLES TO AVG 

Function Number: 4110 

Module: MODUSR2* 


This function declares the number of data samples to be 
included in an average of the hotwire calibration data. See 
LOG DATA POINT for a description of the actual data 
acquisition process. 


[ SAMPLES TO AVG ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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STORE COEFS 


Function Name: STORE COEFS 

Function Number: 4124 

Module: M0DUSR2 


This function stores coefficients previously defined by 
functions ENTER COEFS, LOAD COEFS, COMPUTE COEFS, etc. 

The coefficient filename must have been previously defined - 
see function COEF FILENAME. 


[ STORE COEFS ] 


This routine - MODUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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TAG PICTURE 


Function Name: TAG PICTURE 

Function Number: 4004 

Module: MODUSR1* 


This function "tags'* the picture with the first few 
parameter names and values from the current observation. These 
few values are intended to identify the environment from 
which the "picture” was taken. This function is therefore 
intended to be invoked just after function REDRAW PICTURE or 
PICTURE LOG E. 


[ TAG PICTURE ] 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NASl-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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SELECTOR 


Func t ion Name : 
Function Number: 

Module : 


SELECTOR 

4129 

M0DUSR2 * 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ SELECTOR ] 


*. This routine - M0DUSR2 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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TAG PLOT 


Function Name: TAG PLOT 

Function Number: 4003 

Module: M0DUSR1* 


This function "tags" the plot with the first few parameter 
names and values from the current observation. These few 
values are intended to identify the environment from which the 
plot was taken. This function is therefore intended to be 
invoked just after function REDRAW PLOT or PLOT LOG E. 


[ TAG PLOT ] | 

i 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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TO DISK 


Function Name : 
Func t ion Number : 

Module : 


TO DISK 
4027 

MODUSR1* 


This function defines the mass storage device to which the 
files will be copied - function COPY FILES - or moved - 
function MOVE FILES. 


[ TO DISK ] 

1 I 

I I 

|--[ - ]-*[ MSI ]--| 

I tern 1 Description 

value | MSI | any valid HP 


[...Range. 


| storage 
| device 


*. This routine - M0DUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NASl-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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UTILITY 


J 


Func t ion 
Function 

Module : 


Name : 
Numb e r : 


UTILITY 

4000 

M0DUSR1* 


This function serves merely as a label for the path for 
accessing other related functions via the softkeys. 


[ UTILITY ] 


* 


*. This routine - MODUSR1 - was written by Steven J. Clukey of 
Vigyan Research Associates, Inc. for NASA LaRC TAD/FDB under 
contract NAS1-17919, Task 36. This work was done beginning in 
October, 1986 and continues thru April, 1988. 
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c SYS UAPb NAME-AUTOUARS 
6 LOAD SYS UARS 
^ LOG FILENAME-974 
8 LOAD LOGFILE 
10 PLOT EJECT 
15 GRAPH OFF 

20 TRACE ACTIUE-1 .2 ,3 .4,5 ,6 ,7 

30 C-ALL 

40 C TO WHOLE 

50 X CALIB TYPE-TIME 

60 Y CALIB TYPE-UOLTS 

65 BIN SW NAME-8FTSW 

66 LOAD BIN SW 

^0 SEO PPOG NAME - 3FTRSEQ 

71 DATA DISC DEU - :,1400 

72 F ILENAMEC 1 ) -T9 3 4 0100 

7? F I LENAME ( 2 )-T934 0200 

74 FI LENAME (3 )-T934 0300 

75 F I LENAME C 4T-T934 0400 

76 F ILENAMEC 5 )-T934 0500 

77 F I LENAME ( 6 ) -T934 0600 

78 F I LENAME ( 7 ) -T934 0700 

90 DATA FILE MAPCD-YES 

91 DATA FILE MAPC2)-YES 

92 DATA FILE MAPC3)«YES 

93 DATA FILE MAPC4)-YES 

94 DATA FILE MAPC5)-YES 

95 DATA FILE MAPC6)«YES 

96 DATA FILE MAPC7)«YES 

97 OATA FILE MAP(8)-NO 
"98 DATA FILE MAPC9)-NO 

99 DATA FILE MAPC10)-NO 

100 DATA FILE MAPC11)- NO 

101 DATA FILE MAPC12)-NO 

102 DATA FILE MAPC13)-NO 

103 DATA FILE MAPC14)-N0 

110 MEM LENGTH- 64K 

111 MEM LENGTH C 1 )-64K 

112 MEM START C 1 )«0 
119 MEM START C 0 ) - 179 IK 

121 TF MAP CHANCD-l 

122 TF MAP CHANC2 )-2 

123 TF MAP CHAN ( 3 ) -3 

124 TF MAP CHAN(4>-4 

125 TF MAP CHANC5 )-5 

126 TF MAP LriAN ( 6 ) -6 

127 TF MAP CHANC7W 

128 TF MAP CHAN (8 >-0 

129 TF MAP CHAN (9 )-0 

130 TF MAP CHANC 10 )-0 

131 TF MAP CHANC 11 )-0 

132 TF MAP CHANC 12 )-0 

133 TF MAP CHANC 13 )-0 

134 TF MAP CHANC 14)-0 
200 LOAD SEO PROG 
210 RUN SEO PROG 


APPENDIX C.l. Initialization Sequence Program - "8FTINI" 


7 PLOT EJb.LT 

4 INITIAL OBS-49 

5 PLOT LOG_E 

6 PRNT LOGFILE 

7 SEQUENCE MENU 

8 LET 1-0 

10 PRINT "PRESS ENTER WHEN READY TO LOG DATA" 
20 WAIT ? 

70 ARM 

80 LOG DATA POINT 

65 PRNT LOGFILE 

90 PRINT "POINT IN PROGRESS" 

95 PLOT LOG_E 

96 IF 1-0 THEN TAG PLOT 

97 1-1 

100 WAIT * 

101 PRINT "TRANSFERRING DATA TO MEMORY" 

105 TF FROM REC 

106 PRINT "RECORDING DATA" 

107 STORE DATA 

110 PRINT "POINT COMPLETE " 

800 GOTO 1C 


APPENDIX C.2. Run Sequence Program 


"8FTRSEQ W 
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LIST 
1 LIST 

10 ' PROGRAM XFR.HP S. J. Clukey, Vigyan Research Associates 

20 
30 

40 'Initialization from " EXAMPLE . 3AS " of the HP-IB Command Library 
50 'Copyright Hewlett-Packard 1984, 1985 
50 

70 ' Set up program for MS-DOS HP-IB I/O Library 

80 ' For use independent of the PC instrument bus system 

90 CLS 
100 

110 DEF SEG 
120 CLEAR , &HFE00 
130 I-&HFE00 
140 

150 ' PCIB.DIR$ represents the directory where the library files 

160 ' are located 

170 ' PCIB is an environment variable which should be set from MS-DOS 

180 ' i.e. A:> SET PCIB-A:\LIB 

190 ' 

200 ' If there is insufficient environment space a direct assignment 

210 ' can be made here, i.e 

220 ' PCIB.DIR? - "A:\LIB" 

230 ' Using the environment variable is the preferred method 

240 ' 

250 PCIB . DIR$ - ENVIRON$("PCIB") 

260 1$ - PCIB . DIR$ + '’\PCIBILC . BLD" 

270 BLOAT) IS.&HFE00 

280 CALL I (PCIB . DIR$ , 1%, J%) 

290 PCIB. SEG - 1% 

300 IF J%-0 THEN GOTO 370 
310 PRINT "Unable to load."; 

320 PRINT " (Error 
330 STOP 
340 ' 

350 ' Define entry points for setup routines 

360 ' 

370 DEF SEG - PCIB. SEG 
380 0.S - 5 

390 C.S - 10 
400 I.V - 15 
410 I.C - 20 
420 L.P - 25 
430 LD. FILE - 30 
440 GET. MEM - 35 
450 L.S - 40 
460 PANELS - 45 
470 ' 

480 ' Establish error variables and ON ERROR branching 

490 ' 

500 DEF. ERR - 50 

510 PCIB . ERR$ - STRING$(64, 32) 

520 PCIB. NAMES - STRINGS (16 , 32) 

530 CALL DEF. ERR( PCIB. ERR, PCIB. ERRS, PCIB. NAMES, PCIB. GLBERR) 

540 PCIB . BAS ERR - 255 
550 ON ERROR GOTO 370 
560 ' 

570 J--1 

580 I$-?CIB . DIR$+"\HPIB . SYN" 
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590 CALL O.S(I$) 

500 IF PCIB.ZRRoO THEN ERROR ?CI3 . 3ASERR 
510 ' 

620 ' Determine entry points for HP- IB Library routines 

530 ' 

640 1-0 

650 CALL I . V ( I , IOABORT , I0CLEAR , I0C0NTR0L , IOENTER) 

660 IF PCIB.ERRO0 THEN ERROR PCIB . BAS ERR 
670 CALL I . V(I , IOENTERA, IOENTERS , IOEOI , IOEOL) 

680 IF PCIB . ERRO0 THEN ERROR PCIB . BAS ERR 

690 CALL I . V( I , IOGETTERM , I0LL0CK0UT , IOLOCAL , IOMATCH) 

700 IF PCIB . ERROO THEN ERROR PCIB . BAS ERR 

710 CALL I . V ( I , IOOUTPUT , IOOUTPUTA , IOOUTPUTS , IOPPOLL) 

720 IF PCIB . ERRoO THEN ERROR PCIB . BAS ERR 

730 CALL I . V(I , IOPPOLLC, IOPPOLLU, IOREMOTE, IORESET) 

740 IF PCIB. ERROO THEN ERROR PCIB . BAS ERR 

750 CALL I . V(I , IOSEND, IOSPOLL, IOSTATUS , IOTIMEOUT) 

760 IF PCIB. ERROO THEN ERROR PCIB . BAS ERR 
770 CALL I . V( I , IOTRIGGER , J , J , J ) 

780 IF PCIB. ERROO THEN ERROR PCIB . BAS ERR 
790 CALL C.S 

800 I$— PCIB . DIR$-t-"\HPIB . PLD" 

310 CALL L. P(I$) 

820 IF PCIB. ERROO THEN ERROR PCIB . BAS ERR 
330 GOTO 1000 
340 ' 

850 ' Error handling routine 

360 '■ 

370 IF ERR-PCIB . BAS ERR THEN GOTO 900 

880 PRINT "BASIC error #" ; ERR; " occurred in line " ; ERL 
390 STOP 

900 TMPERR - PCIB. ERR 

910 IF TMPERR - 0 THEN TMPERR - PCIB. GLB ERR 

920 PRINT "PC Instrument error #";TMPERR;" detected at line ";ERL 
930 PRINT "Error: PCIB. ERRS 

940 STOP 
950 ' 

960 ' COMMON declarations are needed if your program is going to chain 

970 ' to other programs. When chaining, be sure to call DEF.ERR as 

980 ' well upon entering the chained- to program 

990 ' 

1000 COMMON PCIB. DIR$, PCIB. SEG 

1010 COMMON LD. FILE, GET. MEM, PANELS, DEF.ERR 

1020 COMMON PCIB. BASERR, PCIB. ERR, PCIB. ERRS, PCIB. NAMES, PCIB. GLBERR 

1030 COMMON IOABORT, IOCLEAR, I0C0NTR0L, IOENTER, IOENTERA, IOENTERS , IOEOI , IOEOL, IOG 

ETTERM , I0LL0CK0UT , IOLOCAL , IOMATCH , IOOUTPUT , IOOUTPUTA , IOOUTPUTS , IOPPOLL , IOPPOLLC , 

IOPPOL 

1040 ' 

1050 FALSE - 0 

1060 TRUE - NOT FALSE 

1070 NOERR - 0 

1080 EUNKNOWN - 100001 ! 

1090 ESEL - 100002! 

1100 ERANGE - 100003! 

1110 ETIME - 100004! 

1120 ECTRL - 100005! 

1130 EPASS - 100006! 

1140 ENUM - 100007! 

1150 EADDR - 100008! 

1160 COMMON FALSE, TRUE, NOERR, EUNKNOWN, ESEL, ERANGE, ETIME, ECTRL, EPASS, EN 
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UM, EADDR 
1170 ' 

1180 ' End Program Sec-up 

1190 ' User program can begin anywhere past "his poinc 
1200 ' Program for a system Co receive data from Che Dynamic Data 
1210 ' Acquisition System. 

1220 ' 

1230 ' 

1240 OPTION BASE 1 

1250 MAX. VARIABLES- 50 

1260 DIM NAMES $ (MAX. VARIABLES) 

1270 DIM X(3) 

1280 ACT. VARIABLES- 0 
1290 NAMES $ - SPACES (50) 

1300 ' 

1310 ' Set up HP- IB addressing and initialize system 

1320 ' 

1330 ISC-7 
1340 DEV-1 

1350 DEV - ISC * 100 + DEV 
1360 CALL IORESET (ISC) 

1370 IF PCIB . ERR O NOERR THEN ERROR PCIB . 3ASERR 
1380 TIMEOUT - 5 

1390 CALL IOTIMEOUT (ISC, TIMEOUT) 

1400 IF PCIB. ERR O NOERR THEN ERROR PCIB . BAS ERR 
1410 CALL IOCLEAR (ISC) 

1420 IF PCIB. ERR O NOERR THEN ERROR PCIB . BAS ERR 
1430 ' 

1440 ' 

1450 ' 

1460 CALL IOEOI (ISC, FALSE) 

1470 IF PCIB. ERR O NOERR THEN ERROR PCIB . 3ASERR 
1480 ' 

1490 ' 

1491 PRINT "DO YOU WISH TO RECEIVE A FILE FROM THE HP COMPUTER? (Y or N) " 

1492 INPUT ANS$ 

1493 IF ANS$— "N" THEN GOTO 1730 

1494 IF ANS$0"Y" THEN GOTO 1491 

1500 PRINT "ENTER THE NAME OF THE FILE TO RECEIVE TRANSFERED DATA: " 

1510 INPUT FIL$ 

1520 OPEN FIL$ FOR OUTPUT AS #1 
1530 FOR J-l TO 2 
1540 MM— 6 : AA-0 

1550 T$(J)— SPACES (MM) 

1560 CALL IOENTERS (DEV,T$(J) ,MM,AA) 

1570 IF PCIB. ERR O NOERR THEN ERROR PCIB. BAS ERR 
1580 N(J)— VAL(LEFT$(T$(J) ,AA-2) ) 

1590 PRINT N(J) 

1600 NEXT J 

1610 NOBS— N(l)+1 : NVARS-N(2) 

1620 FOR J-l TO NOBS 
1630 FOR 1-1 TO NVARS 
1640 M-20 : A— 0 

1650 TEMPS-S PACES (M) 

1660 CALL IOENTERS (DEV, TEMPS, M, A) 

1670 IF PCIB. ERR O NOERR THEN ERROR PCIB . BAS ERR 
1680 NAMES$(I)-LEFT$(TEMP$,A-2) 

1690 ' PRINT "Record ",J.", item ",(1),"- " NAMES$(I) 

1701 IF J-l THEN PRINT *1,NAMES$(I) ELSE PRINT *L , VAL(NAMES$ ( I) ) 

1710 NEXT I 
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1720 NEXT J 

1721 CLOSE =1 

1722 GOTO 1490 
1730 END 

Ok 
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